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A.  PURPOSE 

The  aeries  of  Investigations  reported  herein  were  aimed 
at  gaining  physical  understanding  of  certain  transport  processes 
of  importance  in  electrochemical  systems,  and  particularly  in 
energy  conversion  and  storage  devices.  Emphajls  was  placed 
on  the  elucidation,  of  the  effect  of  hydrodynamic  conditions  on 
the  maximum  feasible  rates  of  anodic  and  cathodic  processes ^ 
and  on  the  clarification  of  the  role  of  diffusion  and  migration 
of  reacting  species  on  the  distribution  of  current  along  elec¬ 
trode  surfaces.  Further,  the  relationships  among  the  various 
transport  properties  were  investigated,  with  the  purpose  of 
developing  techniques  for  the  application  of  available  trans¬ 
port  data  to  electrochemical  systems  of  practical  significance. 

Several  distinct  projects  were  pursued  simultaneously, 
each  of  them  representing  a  Master  of  Science  or  Doctor  of 
Philosophy  thesis  subject. 

Projects  completed: 

I.  Current  distribution  at  a  gas -electrolyte -electrode 
interface . 

II.  Ionic  mass  transport  in  channels  by  combined  free 
and  forced  convection. 

III.  The  role  of  ionic  migration  in  electrolytic  mass 
transport. 

IV.  Mass  tran.-  in  concentrated  binary  electrolytes: 
a  method  for  the  prediction  of  ionic  mobilities. 

V.  The  solubility  and  dlffuslvlty  of  oxygen  in  aqueous 


KOH. 
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Projects  Initiated  under  this  contract,  and  continued 
a:%cr  the  termination  date  under  new  sponsorship:* 

VI.  Ionic  mass  transport  and  current  distribution  in 
channels  in  high  velocity  flow. 

VII.  Generalized  method  for  the  evaluation  of  transport 
properties  in  concentrated  electrolytes. 


Under  the  sponsorship  of  the  Inorganic  Materials  Research 
Division  of  the  Lawrence  Radiation  Laboratory,  Berkeley. 
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B.  ABSTRACT 

Project  I . ;  Current  Distribution  at  a  Qas-Electrode-Electro- 
lyte  In';erface. 

The  variation  of  current  density  along  a  cylindrical, 
partially  immersed  electrode  relative  to  the  position  of  the 
apparent  or  Intrinsic  meniscus  has  been  measured  as  a  function 
of  total  applied  current.  Measurements  of  current  distribution 
are  reported  for  the  cathodic  reduction  of  Og  gas  In  aqueous 
KOH  electrolyte  on  sectioned  Ni  and  Ag  electrodes,  A  thin  film 
of  electrolyte  was  observed  to  exist  above  the  Intrinsic  menis¬ 
cus.  It  was  found  that  the  charge  transfer  which  takes  place 
at  the  electrode  electrolyte  surface  occurs  almost  totally 
above  the  bottom  of  the  intrinsic  meniscus,  and  that  a  large 
portion  of  this  chn.rge  transfer  occurs  in  the  thin  film  as 
high  as  1  to  5  millimeters  above  the  top  of  the  intrinsic 
meniscus.  However,  for  high  total  currents  the  relative  por¬ 
tion  of  current  at  large  distances  above  the  intrinsic  menis¬ 
cus  decreases  while  near  the  intrinsic  meniscus  it  increases. 

On  nickel  the  current  was  more  uniformly  distributed  into  the 
upper  areas  of  the  thin  film  than  on  silver. 

A  mathematical  model  is  presented  which  takes  into  account 
diffusion  and  migration  of  relevant  substances  and  the  solu¬ 
bility  of  oxygen  in  the  electrolyte.  The  model  accounts  for 
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the  concentration  dependence  of  transport  properties.  A  com¬ 
parison  of  the  theoretical  and  experimental  results  indicates 
that  the  current  density  distribution  is  controlled  by  a 
balance  between  charge  transfer  overpotential  and  ohmic  resis¬ 
tance  drop  In  the  electrolytic  film. 

Project  II. :  Ionic  Mass  Transport  by  Combined  Free  and  Forced 

* 

Convection. 

The  effects  of  buoyancy  forces  on  Ionic  mass  transfer  at 
horizontal  planar  electrodes  were  studied  experimentally  wi^h 
both  laminar  and  turbulent  bulk  flow.  Copper  cathodes  and 
an  aclfldled  cupric  sulfate  electrolyte  were  used  to  obtain 
mass  transfer  data  by  the  limiting  current  technique  in  the 
ranges  of  parauneters: 

10  <  Gr  (Grashof  Number)  <  10 

75  <  Re  (Reynolds  N'imber)  <  7,000 

600  <  Sc  (Schmidt  Nujnber)  <  12,000 

0.05  <  ^  (duct  dlam/electrode  length)  <  20 

40  <  Nu.Sh  (Nusselt  Sherwood  Number)  <  1,500 

With  turbulent  bulk  flow,  buoyancy  forces  were  found 
to  have  a  negligible  effect  on  the  mass  transfer  process 
in  the  range  of  variables  covered. 

With  laminar  bulk  flow,  buoyanc’^  forces  can  induce  free 
convective  secondary  flows  near  the  electrode  interface.  If 
log  (Gr  >  2.96  +  log  (Re),  the  process  Is  dominated  by  free 
convection,  and  the  equation  of  Fenech  and  Tobias, 

*  This  project  was  sponsored  by  the  Lawrence  Radiation  Labora¬ 
tory  up  to  March  1961;  by  NOL-Corona  (Contract  No.  N123 
(6738)17178A)  to  April  1962. 
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Nj'  -  0.19  (Gr.Sc)V-‘, 

may  be  used  to  predict  average  mass  transfer  coefficient. 
Otherwise,  forced  convection  dominates,  and  the  equation  of 
Norris  and  Streld, 

Nu'  =  1.85  (Re.Sc.^)^/^, 

yields  acceptable  prediction  of  the  experimental  data.  The 
average  deviation  of  all  the  data  from  the  appropriate  corre¬ 
lations  was  less  than  ±7^. 

The  variation  of  local  mass  transfer  coefficients  paral¬ 
lel  and  normal  to  the  direction  of  flow  may  be  explained  by 
considering  the  Interaction  of  the  two  modes  of  convection. 

Project  III. :  The  Role  of  Ionic  Migration  In  Electrolytic 

Mass  Transport;  Diffusivities  of  [Fe{CN)g] 
and  [Pe(CN)g]^“  in  KOH  and  NaOH  Solutions. 

The  role  of  migration  In  Ionic  mass  transport  Is  dis¬ 
cussed,  and  a  rigorous  method  Is  presented  for  the  evaluation 
of  the  contribution  of  migration  when  a  constant  potential 
gradient  can  be  assumed  to  exist  In  the  mass  transfer  boundary 
layer. 

The  influence  of  supporting  electrolytes  Is  shown  to 
depend  on  the  relative  mobilities  of  all  ionic  species  present 
In  the  solution. 

Theorf  :5al  results  are  applied  to  the  K2[Fe(CN)gJ  ~ 
K^[Fe(CN)g]  -  KOH  or  NaOH  systems.  Limiting  currents  were 
measured  by  a  rotating  disk  electrode,  and  the  apparent  diffu¬ 
sivities  calculated  from  these  measurements  were  corrected 
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taking  into  account  the  influence  of  the  presence  of  noiireacting 


ionic  species.  It  is  shown  that  the  corrected  product,  Dii, 
is  ccnptant  over  a  large  range  of  concentrations 

r  [Fe(CN)g]^",  2.40  x  10"^^ ( crn^/scc )  poise/’K 

for  [Fe(CN)g]^“,  2.02  x  10“^^ ( cm^/sec )  poise/®K. 


Project  IV. :  Mass  Transfer  in  Concentrated  Binary  Electrolytes. 

For  a  multicomponent  system  the  diffusion  law 
can  be  expressed  in  the  form 


^T®iJ 


where  is  the  electrochemical  potential  and  v.  the  average 

a.  TX 

velocity  of  species  1  and  is  a  diffusion  coeificient 
accounting  for  the  interaction  between  species  1  and  J.  This 
alffuslon  law  is  self-consistent  and  can  be  used  to  develoo 
equations  in  a  form  useful  for  electrochemical  transport 
problems . 

It  IS  suggested  that  the  concentration  and  temperature 
dependence  for  ion-solvent  interactions  is  described  by 

=  kT/A, 

and  for  ion-ion  interactions  is 

where  p,  is  the  viscosity  and  A  and  G  are  approximately  constant. 


Project  V. :  The  Solubility  and  Diffusion  Coefficient  of  Oxy¬ 
gen  in  KOH. 

The  solubility  of  oxygen  in  aqueous  KOH  solutions  has 
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been  measured  by  a  Van  Slyke  apparatus  and  by  an  absorption 
technique  developed  by  Hildebrand.  In  the  range  of  concentra¬ 
tion  of  KOH  between  0  and  12  N,  at  25®C,  the  two  methods 
yielded  identical  results: 

log  S  =  .1746  C  +  log  1.26  x  lO"^ 
where  .1746  is  the  solubility  coefficient,  S  =  solubility  of 
oxygen,  gram-moles  per  liter,  and  C  =  concentration  of  KOH, 
gram-moles  per  liter.  The  solubility  was  found  to  decrease 
with  Increasing  temperature,  in  dilute  solutions,  and  to  in¬ 
crease  with  temperature  in  concentrations  of  KOH  higher  than 
8.25  normal. 

Diffusion  coefficients  of  oxygen  in  aqueous  KOH  were 
evaluated  from  the  limiting  current  of  oxygen  on  a  rotating 
disk  electrode,  and  also  by  a  stagnant  tube  technique  similar 
to  that  used  by  von  Stackelberg.  The  diffusivity  drops  sharply 
with  increasing  KOH  concentration,  and  increases  with  tempera¬ 
ture.  At  25®C  and  concentrations  higher  than  2N  KOH,  the 
product  of  the  diffusivity  and  the  viscosity  Is  constant : 

_7 

Dp.  =  1.30  x  10  gram-cm 
2 

where  D  is  the  diffusivity,  cm  /sec,  and  \i  is  the  viscosity 
in  poise.  At  60® C  the  value  of  this  product  is 

Dp.  =  1.88  X  10”^  gram-cm 

and  is  constant  for  concentrations  higher  than  1  mole /liter. 
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Project  VI t  Ionic  Maas  Transport  and  Current  Distribution 
in  Channels  at  High  Velocity  "^low. 

As  a  tool  for  extended  and  refined  investigations  into  the 
nature  of  superimposed  free  and  forced  convective  mass  transfer 
at  electrodes,  a  combination  of  a  high-capacity  liquid  circula¬ 
tion  loop,  a  controlled  current  source  and  an  electric  measure¬ 
ment  system  has  been  designed  and  constructed. 

The  central  element,  the  electrode  flow  cell,  consists  of 

a  rectangular  channel  of  15  by  2  cm  cross-section  containing 

2 

40  cm  long  plane  parallel  electrrdes  of  600  cm  area,  each, 

separated  by  a  porous  diaphragm.  Alternatively,  a  flow  cell  of 

2 

5  by  2  cm  cross-section  containing  two  electrodes  of  180  cm  area, 
each,  can  be  used.  The  cathodes  are  divided  into  Insulated, 
coplanar  sections  which  allows  measurements  in  axial  or  lateral 
directions.  The  electrolyte  circulation  system  can  provide  1 
linear  velocities  along  the  electrodes  of  up  to  450  cm/sec,  cor¬ 
responding  to  Re  =  70,000.  Special  attention  has  been  given  to 
the  design  of  the  inlet  section  of  the  flow  cell. 

A  regulated  d.c.  power  combired  with  a  linear  current  ramp 
generator  has  been  designed  and  constructed.  This  source  can 
provide  any  current  maximum  up  to  100  amp  in  a  preset  time  range, 
between  1  and  20  minutes.  The  section  current  measuring  device 
has  been  modified  to  ensure  equipotentlal  conditions  at  the 
cathode,  within  ±1  mV. 

The  electric  measurement  system  allows  simultaneous  record¬ 
ing  of  total  current,  sectional  current  densities,  local 
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reference  electrode  potentials  and  electrolyte  flow  rate. 

The  entire  cell  system,  power  supply  and  measuring  cir¬ 
cuitry  is  presently  subjected  to  extensive  calibration  and 
evaluation  procedures. 

Project  VII t  Transport  Properties  of  Concentrated  Electrolytic 
Solutions . 

Available  data  on  transport  properties  of  binary  electro¬ 
lytes  have  been  used  to  calculate  coefficients  as  defined 
by  ^,ewrllan,  Bennlon,  and  Tobias.  The  concentration  dependence 
of  conductivities  of  various  systems  corrected  for  viscosity 
Indicates  that  the  development  of  a  more  generally  useful  method 
for* the  prediction  of  transport  properties  in  binary  and  mul¬ 
ticomponent  systems  requires  more  reliable  and  more  extensive 
data  on  diffusion  coefficients. 

Of  the  possible  methods  of  measuring  dlffuslvltles,  the 
rotating  disk  electrode  and  a.  new  type  of  restricted  diffusion 
cell  which  uses  interferometry  for  the  detection  of  concentra¬ 
tion  gradients  are  shown  to  be  most  promising  for  obtaining 
precise  data  because  it  is  possible  to  Interpret  the  measure¬ 
ments  with  rigorous  matxhematlcal  analysis.  Therefore  an  im¬ 
proved  rotating  disk  apparatus  (see  Project  III)  and  a  modi¬ 
fied  restricted  diffusion  cell  with  optical  measurements  which 
applicable  to  concentration  solutions  are  presently  under 
development . 
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1 .  Lectures : 
Tobias,  C.  W. 


Tobias,  C.  W. 


Newman,  J.  S. 


Bennlon,  D.  N 


HI ckman ,  R .  G 


LECTURES,  CONFERENCES  AND  REPORTS 

"The  Solubility  and  Diffusion  Coefficient  of 
Oxygen  In  KOH" j  contribution  to  the  discussion 
of  oxygen  electrodes,  on  the  Sixth  Advisory 
Group  for  Aeronautical  Research  and  Development 
Combustion  and  Propulsion  Colloquium,  Cannes, 
France,  March  15-18,  1S64. 

"Convective  Mass  Transport  In  Electrolytes"; 
Symposium  on  Mass  Transport  In  Electrochemical 
Processes,  a  keynote  address  at  the  125th 
National  Meeting  of  The  Electrochemical  Society, 
Toronto,  Canada,  May  3-7,  1964. 

"The  Plow  of  Charge  and  Mass  Through  an  Agi¬ 
tated  Electrolyte";  Symposium  on  Mass  Trans¬ 
port  In  Electrochemical  Processes  at  the  125th 
National  Meeting  of  The  Electrochemical  Society, 
Toronto,  Canada,  May  3-7,  1964. 

.  "Current  Distribution  at  a  Gas -Electrode - 

Electrolyte  Interface";  Symposium  on  Fuel  Cells 
at  the  126th  Meeting  of  The  Electrochemical 
Society,  Washington,  D.C.,  October  11-14,  1964. 

.  "The  Effect  of  Buoyancy  Forces  on  Forced  Con¬ 
vection  Ionic  Mass  Transfer  at  Horizontal 


Eler  erodes" ;  Symposium  on  Mechanism  of  Electro- 
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deposition  at  IPSth  ^'t-etlng  of  The  Electro¬ 
chemical  Society^  Washington,  D.  C.,  October 
11-14,  19G4. 


2 .  Conferences ; 

June  23,  1962  -  Port  Monmouth,  New  Jersey. 

Organization  represented:  U.S.  Army  Electronics  Re¬ 
search  and  Development  Laboratories;  Department  of 
Chemical  Engineering,  University  of  California,  Berkeley. 
The  meeting  was  held  to  review  the  status  and  future 
plans  for  this  research  program. 

April  17,  1963  -  Fort  Monmouth,  New  Jersey. 

Organizations  represented:  U.S.  Army  Electronics  Re- 
se.irch  and  Development  Laboratories;  Department  of 
Chemical  Engineering  (Charles  W.  Tobias),  University  of 
Calliornia,  Berkeley.  The  meeting  was  held  'to  review 
the  status  and  future  plans  for  this  research  program. 
August  20-21,  1963  -  University  of  California,  BerkeDey. 
Organizations  participating:  U.S.  Array  Electronics  Re¬ 
search  and  Development  Laboratories  (Mr.  John  Murphy); 
Department  of  Chemical  Engineering,  University  of  Cali¬ 
fornia,  Berkeley. 

August  6,  1964  -  University  of  California,  Berkeley. 

Organizations  participating:  U.S.  Army  Electronics  Re¬ 
search  and  Development  Laboratories  (Mr.  John  Murphy); 
Department  of  Chemical  Engineering-  University  of  Cali¬ 
fornia,  Berkeley. 
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3.  Reports ; 

Quarterly  letter  reports: 

Report  No.  1-30  June  1962 
Report  No.  2  -  30  September  1962 
Report  No.  3-31  December  1962 
Report  No.  4-31  March  1963 
Report  No.  5-30  June  1963 
Report  No.  6-30  September  1963 
Report  No.  7-31  December  1963 
Report  No.  8-31  March  1964 
Report  No.  9-30  June  1964 
Report  No.  10  -  30  September  1964 
Report  No.  11  -  31  December  1964 
Report  No.  12  -  31  March  1965. 

Technical  Reports: 

Report  No.  1.  "Investigation  of  Ionic  Diffu¬ 
sion  and  Migration  by  a  Rotating  Disk  Electrode",  by 
Stanley  L.  Gordon  (M.S.  Thesis),  Issued  May  15,  1963. 

Report  No.  2.  "The  Effect  of  Buoyancy  Forces 
on  Forced  Conveciilon  Ionic  Mass  Transfer  at  Horizontal 
Planar  Electrodes",  by  Robert  G.  Hickman  (Ph.D.  Thesis), 

Issued  December  1,  1963. 

Report  No.  3.  "Phenomena  at  a  Gas-Electrode- 

Electrolyte  Interface",  by  douglas  N.  Bennlon  (Ph.D.  Thesis) 

Issued  June  1,  1964. 
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Publications ; 

The  results  of  investigations  under  this  contract,  have  not 
as  yet  been  published  in  the  open  scientific  literature. 

Papers  with  contents  closely  similar  to  the  descrlpt^-ons  of 
projects,  will  be  submitted  to  Journals  following  the 
issuance  of  this  report. 
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PROJECT  I. 

Current  Distribution  at  a  Gas  - 
Electrode -Electrolyte  Interface 


a)  Experimental  Treatment 


Introduction 

Many  of  the  possible  fuels  and  oxidants  proposed  for  use 
in  fuel  cells  are  gases.  To  gain  an  understanding  of  elec¬ 
trodes  In'^olvlng  gase<,as  reactants,  analysis  of  transport 
phenomena  at  gas-electrcde-electrolyte  interfaces  Is  necessary. 
Various  stuoies  of  porous  media  involving  only  electrode- 
electrolyte  Interfaces  have  been  presented  recently  . 

The  presence  of  a  gas  phase  provides  additional  problems  which 
are  to  be  considered  here. 

Porous  electrodes  are  used  to  Increase  the  reaction  area 
per  unit  volume  and,  in  the  case  of  gaseous  reactants,  to  pro¬ 
vide  separation  between  gas  and  electrolyte.  In  such  elec¬ 
trodes,  many  simultaneous  and  crnsecutive  processes  occur, 
and  the  nature  and  extent  of  the  actual  reaction  area  within 
the  pores  is  not  well  defjned  or  understood.  A  plane  elec¬ 
trode  partially  immersed  in  an  electrolytic  solution  offers 
a  simpler  system  for  studying  the  behavior  of  ga;.  electrode - 
electrolyte  interfaces. 

In  1957  Wagner^  suggested  that  a  thin  film  of  electrolyte 
might  exist  above  the  visible  meniscus  on  the  surface  of  a 
partially  immersed  electrode.  From  his  calculations  it 
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appealed  as  if  the  electrode  surface  contacting  this  thin 
film  might  contribute  appreciably  to  the  charge  transfer  reac¬ 
tion. 

Weber,  Meissner,  and  Sama^  performed  experiments  using 
an  oxygen  cathode  to  study  the  reaction  rate  distribution. 

They  did  not  consider  the  possibility  of  a  thin  film  above  the 
Intrinsic  or  readily  visible  meniscus,  and  copper  oxide 

formation  may  have  Influenced  their  experimental  results. 

7  8 

will  ’  has  perforraec  experiments  with  a  h^^drogen,  sul¬ 
furic  acid,  and  platinum  system.  For  a  give  applied  over¬ 
potential,  he  observed  the  variation  in  current  with  the 
length  of  platinum  extending  above  the  liquid  level  of  the 
electrolytic  solution.  Will  concluded  that  a  thin  film  did 
exist  above  the  intrinsic  meniscus  and  that  most  of  the 
current  was  transferred  into  the  film  in  a  narrow  region 
adjacent  to  the  upper  edge  of  the  meniscus.  He  also  shov/ed 
that  surface  migration  of  hydrogen  along  the  platinum  was 
not  a  significant  mode  of  mass  transfer.  Diffusion  of  hydro¬ 
gen  through  the  film  was  determined  to  be  the  rate  controlling 
step.  This  was  the  first  really  definitive  experlmens,  to  show 
the  role  of  thin  electrolyte  films.  However,  as  the  platinum 
electrode  was  extended  further  and  further  out  of  the  electro¬ 
lyte,  new  area  on  which  an  electrolyte  film  could  exist  was 
produced.  Thus,  the  film  geometry  and  resulting  current 
distribution  clianged  during  the  course  of  an  experiment. 

Our  work  has  been  undertaken  to  obtain  a  more  quantita¬ 
tive  descilptlcn  of  the  reaction  rate  distribution  and  to 
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investigate  the  role  of  the  thin  film  during  oxygen  reduction. 
For  the  experimental  work,  oxygen  gas,  aqueous  3.5  molar  KOK, 
and  either  nickel  or  silver  electrodes  were  chosen. 
Experimental  Equipment  and  Procedures 

The  objective  was  to  measure  the  variation  of  current 
density  along  a  partially  Immersed  electrode  relative  to  the 
bulk  electrolyte  surface.  In  order  to  accomplish  this, 

9 

sectioned  electrodes  were  employed.  The  sections  must  be 
insulated  electrically  from  one  another,  and  the  insulating 
separators  should  no':  significantly  disturb  the  geometry  of 
the  thin  liquid  film  on  the  surface  of  the  electrode.  Also, 
the  sections  of  the  metal  electrode  must  be  at  the  same 
p  -ential. 

In  order  to  satisfy  these  requirements  as  closely  as 
possible,  only  two  sections  and  a  single  separation  were 
used  (see  Figure  l) .  With  this  approach,  the  variation  in 
the  current  from  each  section  could  be  measured  as  a  function 
of  the  distance,  AL,  of  the  bulk  electrolyte  level  below  the 
separation  (see  Figure  1) .  The  sum  of  the  two  currents 
equals  '-he  total;  therefore,  only  ore  of  the  currents  is 
lndepenue.it.  In  order  to  compare  the  results  of  different 
total  currents  more  easily,  the  percent  curr-ent  out  of  the 
top  section,  defined  as  ^.O.T.  =  (Current  out  top  section) 
(10C)/Total  Currenc),  was  used  as  the  Independent  variable. 

The  experimental  procedure  Involved  determination  of 
the  ^’C.O.T,  a?  a  function  of  AL  at  various  total  currents 
on  nickel  and  silver  electrodes.  To  determine  the  actual 
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current  density  variation  one  plots  the  ^.O.T.  vs.  AL  and 
takes  the  slope  of  this  curve.  The  slope  at  a  given  value  of 
AL,  in  cm~^,  times  the  applied  current,  I  in  mA/cid,  gives  the 
transfer  current  density  at  a  height  AL  above  the  electrolyte 
level.  The  readings  were  consistent  to  within  3^  for  a 
single  run  except  at  the  lowest  current  levels  when  varia¬ 
tions  up  to  10^  of  the  measured  current  were  Inherent  in 
the  circuit  design. 

To  avoid  edge  effects  a  cylirdrlcal  electrode  was  used. 

The  working  length  of  the  electrode  was  approximately  1?  cm; 
tnl  was  found  to  be  long  enough  to  avoid  appreciable  end 
effects  at  the  top  and  bottom.  The  electrode  was  constructed 
using  a  cylindrical  plexiglass  spindle  over  which  two  cylin¬ 
ders  of  nickel*  (or  silver**)  were  slipped  (see  Figure  2), 

These  were  separated  hv  a  thin  "gasket*’  made  from  a  sheet  of 
oriented  polystyrene  25  microns  thick.  This  thickness  was 
necessary  to  insure  that  no  metallic  electrical  connection 
occurred  between  the  two  sections. 

In  order  to  avoid  the  possibility  of  systematic  errors 
and  to  confirm  the  reproducibility  of  the  data,  2  different 
silver  cells  and  5  separate  nickel  cells  were  built.  The 
procedures  in  constructing  the  separations  in  the  working 
electrodes  (such  as  rounding  the  corners  of  the  electrode 
where  it  contacts  the  insulating  jacket)  were  modified 
*  . 

The  nickel  vjas  99.5%  pure.  The  principal  impurities 

were  Mn  and  Fe . 

** 


The  silver  was  99.95%  pure. 
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slightly  to  see  if  the  results  were  affected.  No  effects 
on  the  results  were  noted  from  these  changes.  For  one  of 
the  silver  electrodes,  runs  were  also  made  with  the  25  micron 
thick  polystyrene  gasket  replaced  by  a  separator  made  from 
filter  paper  approximately  70  microns  thick.  This  was  done 
to  check  the  results  for  any  dependence  on  the  nature  of  the 
material  used  as  the  separator. 

The  cell  consisted  of  two  cylindrical  compartments  as 
shown  In  Figure  2.  The  cylindrical  electrode  spindle  des¬ 
cribed  above  was  screwed  Into  the  left-hand  compartment. 

In  the  right-hand  compartment,  nickel  screens  served  as  the 
counter  electroc^e.  The  reaction  at  the  working  electrode  was 

Og  +  SHgO  +  4e~  =  40H". 

The  reaction  at  the  counter  electrode  was  just  the  reverse 

of  the  above  reaction.  The  reactions  were  driven  externally 

* 

by  a  constant  current  power  supply. 

The  two  compartments  mentioned  above  were  connected  to 
a  storage  reservoir,  not  shown  In  F"*  ?ure  2.  By  carefully 
adjusting  the  height  of  the  storage  reservoir,  the  liquid 
level  In  the  cell  compartments  couid  be  accurately  adjusted. 
The  value  of  AL  was  determined  by  measuring  the  height  of 
the  bulk  liquid  and  the  height  of  the  separation,  to  ±0.05  .run, 
using  a  cathoraeter. 

The  electrolyte  entered  the  working  electrode  compart¬ 
ment  through  radial  holes  drilled  in  the  lower  part  of  the 

■K 

Electronics  Measurements,  Inc.,  Constant  Current  Power 
Supply,  Model  C612. 
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plexiglass  spindle.  This  provided  for  an  axially  symmetric 
electric  field.  To  align  the  electric  field  more  nearly 
parallel  to  the  electrode  and  to  reduce  convection  effecfs, 
a  precision-bore  glass  tube  was  placed  over  the  electrode 
as  shown  In  Figure  2 .  The  bulk  electrolyte  was  confined 
In  an  annular  space  about  1  mm  wide .  The  bottom  of  the 
meniscus  was  used  as  the  reference  level  for  measuring  AL. 

The  distance  from  the  bottom  of  the  meniscus  to  the  topi  of  the 
apparent  or  ''Intrinsic”  meniscus  was  In  the  range  of  0.6  - 
1.1  mm.  The  electrolyte  In  the  cell  outside  the  jacket  did 
not  participate  In  the  reaction  and  need  not  be  considered. 

Convection  currents  In  the  electrolyte  below  the  menis¬ 
cus  would  tend  to  give  abnormally  high  currents  from  the 
bottom  section.  The  glass  jacket  v/hlch  confined  the  elec¬ 
trolyte  next  to  the  electrode  In  an  annular  region  reduced 
this  effect.  Smaller  spacing  would  have  reduced  the  proba¬ 
bility  of  convection  further^^,  but  this  would  have  caused 
variations  In  the  annular  space .  The  value  chosen  was  con¬ 
sidered  to  be  a  suitable  compromise. 

The  readings  were  sensitive  to  the  surface  conditions 
of  the  electrodes.  If  any  oxide  formed  on  either  the  nickel 
or  silver  electrodes,  the  results  were  not  reproducible. 

Care  was  always  taken  thatt,  following  polishing*,  the  elec¬ 
trode  was  immediately  placed  In  the  cell  and  cathodlcally 

Of  the  various  polishing  compour.  Is,  "Flow-five"  aluminum 
oxide  proved  to  be  the  most  satisfactory. 

Although  the  roughness  of  surfaces  was  different  depending 
on  the  polishing  compounds  used,  no  significant  differences 
In  the  pattern  of  current  distribution  resulted. 
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polarized  to  prevent  oxide  formation.  The  current  was  In¬ 
creased  to  approximately  200  ma  yielding  vigorous  hydrogen 
evolution.  During  the  24-48  hour  cleaning  period,  the 
electrolyte  was  changed  each  2-5  hours.  The  cleanliness 
of  the  surface  was  considered  satisfactory  when,  upon  lowering 
the  electrode  level,  the  surface  remained  uniformly  wetted. 

Electrolytic  contact  was  made  to  a  Hg,  HgO  reference 
electrode  through  a  hole  drilled  lengthwise  down  the  spindle 
and  out  through  a  small  hole  at  the  bottom  of  the  metal  elec¬ 
trode^^.  Electrical  connection  to  the  bottom  section  of 
the  electrode  was  made  with  a  nickel  rod  down  a  hole  in  the 
spindle.  A  stainless  stell  screw  went  through  the  electrode 
and  was  threaded  Into  the  nickel  rod.  The  electrical  connec¬ 
tion  to  the  top  section  was  made  by  a  nickel  plate  In  the 
bottom  of  the  cap.  When  the  cap  was  screwed  down,  the  nlL,kel 
plate  pressed  down  on  the  upper  electrode  section  and  completed 
the  contact.  The  resistances  due  to  these  connections  were 
less  than  0.1  ohm. 

Figure  3  shows  a  schematic  drawing  of  the  electrical 
circuit.  Before  the  measurements  were  started,  the  variable 
resistors  were  set  at  zero;  the  two  electrode  sections  were 
then  shorted  togather  and  thus  were  at  the  same  potential. 

When  measurements  were  to  be  made,  the  resistances  were  in- 
creased  until  a  measurable  voltage  drop  was  obtained.  When 

properly  adjusted,  the  potential  drop  across  each  resistor 

* 

General  Radio  type  1432 -U  or  1432 -K  decade  resistance 
boxes  were  used. 

This  voltage  drop  was  measured  using  a  Keithly  Model 
149  milli-mlcrovoltmeter . 
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was  equal.  As  small  a  voltage  drop  as  practical  was  used 
so  that  errors  In  balancing  wou‘'.d  be  negligible.  These  poten¬ 
tial  differences  varied  between  10  and  1000  microvolts  for 
the  lowest  and  highest  total  applied  currents.  From  the 
voltage  drop  and  the  known  resistance,  the  current  from  each 
section  could  be  calculated.  The  poteritlal  of  the  working 

electrode  with  respect  to  the  Hg,  HgO  reference  cell  was 

* 

measured  with  an  electrometer  . 

The  start-up  procedure  was  to  refill  the  cell  with  3.5 
molar  KOH,  begin  bleeding  In  Og  saturated  with  water  with 
respect  to  3.5  molar  KOH,  and  to  set  the  total  current  at  a 
fixed  level.  The  liquid  level  was  then  lowered  to  below  the 
lowest  value  of  AL  desired.  The  cell  was  left  for  about  15 
hours  to  allow  the  liquid  film  to  drain  and  reach  an  equili¬ 
brium  thickness.  The  level  was  next  raised  a  small  distance 
and  a  set  of  readings  taken.  After  each  set  of  readings, 
the  level  was  raised  again  and  In  this  manner  a  series  of 
readings  for  different  values  of  AL  were  obtained.  Observa¬ 
tions  were  made  at  thirty-minute  Intervals  which  was  suffi¬ 
cient  to  allow  steady  state  conditions  to  prevail  during 
readings.  At  each  electrolyte  level  the  value  of  ^.O.T. 
could  easily  be  calculated. 

Experimental  Results 

12 

Exploratory  experiments  showed  that  srortly  after  the 

electrolyte  level  was  lowered,  a  thin  electrolyte  film  was  left 
* 

Kelthly  Model  610R  electrometer  was  used  here.  Its 
Internal  resistance  Is  rated  at  1014:  ohms  or  greater. 
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behind  on  the  exposed  electrode,  the  fraction  of  current  from 
the  upper  section  was  very  high  but  that  It  dropped  off  with 
time.  This  decrease  with  time  continued  for  up  to  90  hours, 
after  which  a  steady  state  value  remained  as  long  as  the  run 
was  continued,  for  several  more  days  In  some  cases.  This 
behavior  was  attributed  to  the  draining  of  the  electrolytic 
film  which  was  Initially  thick  but  became  thinner  as  draining 
progressed.  Thicker  films  have  a  small  resistance  and  thus 
tend  to  allow  penetration  of  the  current  higher  into  the  film 

Simple  optical  observations  were  made  in  an  attempt  to 
observe  the  film  and  determine  Its  thickness  by  light  inter¬ 
ference  techniques.  The  existence  of  the  film  was  confirmed 
by  the  observation  of  v/eak  interference  colors  when  white 
Incident  light  was  used  and  by  the  observation  of  a  "wake" 
when  a  small  corner  of  tissue  paper  was  drawn  over  the  sur¬ 
face.  Quantitative  measurements  of  the  film  thickness  were 

“  15 

beyond  the  scope  of  this  work.  Muller  has  made  a  detailed 
Investigation  of  the  optical  piopertles  of  thin  electrolytic 
films  on  nickel  electrodes  and  found  that  In  the  range  of 
2  to  5  mm  above  the  meniscus  the  film  thickness  Is  In  the 
order  of  1  micron  for  3.5  molar  KOH. 

Plots  showing  ^C.O.T.  as  a  function  of  AL  are  shown  In 
Figures  4,  5,  and  6  for  4.84,  15.4,  and  48.4  jiA/cm  applied 
current,  respectively.  Figure  7  shows  the  observed  depen¬ 
dence  of  the  electrode  potential  on  the  applied  current. 

A  feature  to  be  noted  Is  that  at  higher  currents  the 
reaction  tends  to  be  concentrated  nearer  the  top  of  the 
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intrinsic  meniscus.  A  comparison  of  the  results  for  equal 
appf^ied  currents  shows  that  the  reaction  on  nickel  spreads 
more  into  the  film  than  on  silver. 

One  of  the  requirements  of  the  sectioned  electrode  was 
that  the  separation  should  not  disrupt  the  film.  This  becomes 
especially  important  when  it  is  realized  that  the  film  is 
about  1  micron  thick  and  the  separation  is  25  microns  wide. 
Such  a  wide  gasket  obviously  causes  some  distrotlon.  The 
distortion  can  be  broken  down  into  three  categories  as  follows 

1)  Assume  that  the  insulating  gasket  is  perfectly  flush 
with  the  metal  surface  and  that  the  film  thickness  does  not 
vary.  The  only  error  from  such  an  ideal  situation  results 
from  no  electrode  reaction  occurring  along  the  insulating 
section.  This  Inactive  length  of  the  film  will  act  as  an 
"unnatural”  resistance.  If  this  occurs  near  the  bottom  of  the 
film  where  the  current  might  reach  100  |iA/cm  as  a  maximum 

and  if  the  film  is  one  micron  th.’’ tae  "unnatural"  resis¬ 
tance  will  cause  a  potential  drop  of  4.2  millivolts.  This 
compares  to  a  total  potential  drop  in  the  film  of  about  140 
millivolts . 

2)  Referring  to  case  l),  the  net  effect  of  the  separa¬ 
tion  might  be  such  thut  the  film  appears  to  be  thicker  as 

it  passes  the  separation.  For  this  situation,  the  "unnatural" 
resistance  will  be  lowered,  and  the  error  will  be  less  than 
for  case  l). 

3)  The  situation  might  also  be  an  effective  thinning  of 
the  film  as  it  passes  the  separation.  At  worst,  this  thinning 
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will  occur  along  the  Tull  length  of  the  separation.  A  thinning 
of  this  type,  for  example,  to  one -fourth  of  the  normal  thick¬ 
ness,  will  Increase  the  "unnatural"  resistance  so  that, 
using  the  example  of  case  l),  the  added  potential  drop  would 
be  19  millivolts.  Such  an  added  potential  drop  would  cause 
a  shift  In  the  results,  such  that  less  current  would  pass  to 
the  upper  section  than  on  an  uninterrupted  surf ace . 

The  reproducibility  of  the  data  was  confirmed  by  using 
different  cells  of  slightly  modified  design.  Similarly,  no 
change  In  the  results  was  detected  by  replacing  the  25  micron- 
thick  polystyrene  separator  with  a  70  micron  thick  separator 
made  from  filter  paper  (see  Figures  and  6). 

Prom  these  observations  It  is  concluded  that  the  film  on 
the  separator  was  not  significantly  thinner  than  on  the  ad¬ 
jacent  metal  surfaces,  and  therefore,  the  maximum  ^rror  In 
the  current  out  the  tope,  due  to  the  presence  of  the  separa¬ 
tor,  is  estimated  to  be 
Discussion 

These  experiments  have  confirmed  that  a  stable,  thin 
film  of  electrolyte  adheres  to  the  electrode  surface  above 
the  top  of  the  intrinsic  meniscus.  A  large  part  of  the  reac¬ 
tion  takes  place  on  the  surface  of  the  electrode  exposed  to 
this  film.  The  reaction  rate  drops  off  very  rapidly  wlthln 
the  visible  meniscus  and  Is  very  small  along  the  electrode 
are?,  exposed  to  the  bulk  electrolyte.  The  obvious  reason 
for  the  decrease  In  reaction  rate  In  the  bulk  electrolyte  Is 
mass  transport  limitation  of  the  molecular  oxygen.  However, 
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the  controlling  feature  In  the  film  region  can  arise  from  two 
sources . 

7  8 

will  •’  concluded  that  the  rate  limiting  step  was  the 
mass  transfer  of  the  gaseous  reactant j  hydrogen  in  this  case, 
through  the  film  to  the  electrode  surface.  He  also  found 
that  the  reaction  rate  dropped  off  very  rapidly  above  the 
"top  of  the  intrinsic  meniscus."  The  observations  with  tne 
oxygen  -  KOH  system  show  that  the  reaction  rate  remains  appre¬ 
ciable  for  several  millimeters  above  the  "top  of  the  intrinsic 
meniscus."  In  additiom,  it  is  observed  that  the  reaction 
extends  further  up  the  film  for  nickel  than  for  silver. 

It  is  also  possible  that  the  limiting  step  can  be  charge 

transfer  at  the  electrode  surface.  It  is  known  that  silver 

is  a  better  catalyst  for  the  reduction  of  oxygen  than  is 
14 

nickel  .  Thus  one  would  expect  that  if  the  dependence  of 
charge  transfer  "overpotential"  on  current  density  were  higher 
for  nickel  than  for  silver,  the  reaction  would  tend  to  "spread 
out"  or  be  distributed  over  a  wider  area,  thus  further  up 
into  the  film. 

Both  the  charge  transfer  and  mass  transfer  overpotential 
tend  to  cause  the  reaction  to  spread  out  farther.  It  is  the 
resistance  drop  in  the  film  which  limits  the  extention  of 
the  reaction  up  into  the  film.  From  either  standpoint,  if 
the  film  were  thicker,  the  reaction  would  tend  to  extend  fur¬ 
ther  into  the  film.  Observations  made  during  the  initial 

♦ 

stages  when  the  film  was  draining  confirm  this  argument. 

♦ 

Measurements  of  film  thickness  on  silver  were  not  avail¬ 
able  at  the  conclusion  of  this  work. 
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* 

Cone lus Iona 

1)  The  existence  of  a  thin  film  of  electrolyte  above 
the  intrinsic  meniscus  on  clean  nickel  and  silver  has  been 
confirmed.  The  thickness  of  this  film  is  established  as 
being  in  the  order  of  magnitude  of  the  wave  length  of  visible 
light. 

2)  For  the  cathodic  oxygen  electrode,  current  density 
distribution  measurements  on  cylindrical  nickel  and  silver 
electrodes  established  that  a  large  portion  of  the  current 
arises  in  the  thin  electrolytic  film  above  the  intrinsic 
meniscus.  The  fraction  of  current  transferred  as  much  as 
10-20  mm  above  the  bottom  of  the  intrinsic  meniscus  is  not 
negligible . 

3)  The  thicker  uhe  film  and  the  smaller  the  applied 
current,  the  further  the  current  spreads  over  the  electrode 
surface  covered  by  the  film.  At  a  given  value  of  applied 
current,  the  reaction  is  concentrated  nearer  the  intrinsic 
meniscus  on  silver  than  on  nickel  surfaces. 


After  the  completion  of  this  manuscript,  an  account  of  the 
work  of  Yu  A.  Mazltov  e_t  al^^  has  come  to  the  author's  atten¬ 
tion.  Current  distribution  along  a  silver  v/lr-e  in  the  reduc¬ 
tion  of  oxygen  in  10.6  KOH  was  measured  by  the  method 
originally  used  by  Will'»°.  Although  their  experimental  ob¬ 
servations  are  in  qualitative  agreement  with  the  results 
presented  above,  the  reader's  attention  is  called  to  the  use 
of  an  uninterrupted  cathode  and  the  much  shorter  draining 
times  (3-4  minutes)  reported  by  Mazltov  e^  al . 
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Figure  2 .  Cross-secticnai  view  of  cell  used  for  measurement  of 

distribution  of  current  relative  to  position  of  meniscus. 
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Figure  4.  Penetration  of  current  in  the  electrolyte  film.  I  =  4.84  pA 
Upper  curve;  Nickel  cathode,  cells  no.  2  and  3,  runs  6  and 
Lower  curve:  Silver  cathode,  cell  no.  4,  runs  no.  6  and  15. 
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Figure  5.  Penetration  of  current  in  the  electrolyte  film.  I  =  15.4  nA/cm 
Upper  curve:  Nickel  cathode,  cells  no.  2  and  3,  runs  5  and  5'. 
Lower  curve:  Silver  cathode,  cell  no.  4,  runs  no.  5,  11  and  14 
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Figure  6.  Penetration  of  current  In  the  electrolyte  film.  I  =  48.4  MA/cm. 

Upper  curve:  Nickel  cathode,  cells  2,  3,  ajod  7;  runs  4,  4',  19,  2C 
Lower  curve:  Silver  cathode,  cells  4  and  6,  runs  4,  10,  13,  19.  20 


Figure  7.  Potential  of  the  oxygen  cathode  (not  corrected  for  ohmic  drop) 

relative  to  the  Hg/HgO  reference  electrode.  Upper  curve:  Sll'.er 
cathode,  lower  curve:  Oxygen  cathode. 
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Current  Distribution  at  a  Gas  - 
Ele ctrode  -E.le  ctroly te  Interface 

b)  Theoretical  Treatment 

Introduction 

An  experimental  study ^  In  this  laboratory  has  shown  that 
on  partially  Immersed  gas  electrodes  a  thin  film  of  electro¬ 
lyte  may  exist  on  the  electrode  above  the  bulk  electrolyte, 
and  when  this  film  exists  the  electrochemical  reaction  takes 
place  primarily  on  the  surface  of  the  electrode  which  contacts 
the  thin  film. 

2 

By  considering  several  possible  physical  models  Wagner  , 

In  1957,  demonstrated  the  likelihood  of  the  Importance  of  thin 

electrolyte  films  In  actual ^s  electrodes.  Since  then, 

various  models  have  been  proposed  which  attempt  to  duplicate 

3  4 

more  closely  actual  physical  conditions.  Will  '  was  able 
to  explain  his  experimental  results  obtained  on  partially  im¬ 
mersed,  plane  electrodes  In  terms  of  a  theoretical  model. 

His  model  did  not  Include  consideration  of  charge  transfer 
overpotential,  and  the  only  transport  property  of  the  elec¬ 
trolyte  employed  was  the  conductivity  which  he  assumed  to  be 

5 

Invariant.  Grens  In  a  model  geometrically  similar  to 

the  one  by  Will  Included  an  Erdey-Gruz-Volmer  type  charge 
transfer  overpotential  expression.  Transport  of  Ions  by  both 
diffusion  and  migration  Is  accounted  for,  and  the  movement  of 
water  is  assumed  to  occur  only  in  the  electrolyte  phase. 

Grens  et  al  were  able  to  show  how  systems  can  become  charge 
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transfer  or  mass  transfer  limited  depending  on  the  magnitude 
of  the  physical  parameters  Involved.  Altnough  variations  of 
transport  properties  were  not  considered,  their  model  predicts 
large  concentration  gradients  in  the  thin  film.  It  should  be 
noted  that,  of  these  models,  only  that  of  Will  has  been  com¬ 
pared  directly  to  experiments,  and  that  his  model  Incorporated 
features  characteristic  of  the  hydrogen,  sulfuric  acid,  plati¬ 
num  system. 

The  above  theories  provide  primarily  a  guide  to  under¬ 
standing  the  detailed  behavior  of  the  electrode  reaction 
v;lthln  and  near  a  thin  electrolyte  film  which  may  exist  on 
plane,  metal  electrodes.  Other  models  have  been  proposed 
which  attempt  to  explain  the  behavior  of  porous  electrodes  in 
terms  of  certain  assumed  microstructures  of  the  porous  natrix. 
Some  recent  models  by  Gurevich^,  Rockett"^,  Iczkowskl®,  and 

9 

Burshteln  e^  al^  ,  have  Incorporated  the  possibility  of  thin 
film  existing  on  the  walls  of  the  gas  filled  pores  or  cava- 
ties  within  the  electrode.  Of  these,  the  one  by  Burshteln 
^  ^  is  of  particular  interest,  since  it  includes  a  trans¬ 
formation  of  a  rather  sophisticated  model  into  one  where 
equations  developed  for  two  phase  porous  electrodes  can  be 

applled^^^^^ .  In  Burshteln  et  al*b  model,  a  detailed 

/  \ 

descrlptlop  of  what  is  talking  place  within  the  thin  film  is 
not  under tp.ken,  and  a  simple,  linear  relation  between  the 
local  ovejl’potential  and  the  current  being  transferred  is 

V 

assumed . 

The  theorejbical  model  to  be  presented  here  is  one  in 

\ 


Section  D.,  Project  I,  part  b) 


Page  40 . 


which  the  details  of  the  film  behavior  are  described.  Our 
model  is  to  be  for  the  oxygen,  KOH,  silver  system  and  is  to 
Include  a  realistic  charge  transfer  overpotential  relation. 

It  also  allows  for  variation  of  electrolyte  composition  and 
consequent  changes  in  the  solubility  of  oxygen  and  in  the 
transport  properties  of  oxygen  and  KOH.  Consideration  is 
given  to  whether  water  transport  occurs  in  the  gas  phase  or 
the  liquid  phase.  Direct  comparison  of  the  theoretical 
model  and  experimental  results  is  designed  to  further  our 
understanding  of  the  system  and  its  controlling  features. 

The  mathematical  model 

Geometric  considerations:  The  model  describes  an  oxygen 
half  cell  operating  cathodically  in  aqueous  KOH  electrolyte. 
From  the  analysis  of  experimental  results,  it  has  been  con¬ 
cluded  that  a  large  portion  of  the  reaction  takes  place  in  a 
thin  film  which  extends  above  the  bulk  electrolyte  on  the 
surface  of  the  electrode.  Measurements  of  the  film  thickness 
by  Muller""  show  that  the  thickness  of  this  film  decreases 
with  height  about  the  bulk  liquid.  However,  the  percent 

change  in  thickness  is  small  compared  to  the  observed  percent 

15 

changes  in  the  reaction  rate  over  comparable  distances 
In  our  model,  to  avoid  unnecessary  mathematical  difficulties, 
the  thin  film  is  assumed  to  have  a  constant  thickness,  6. 

Since  the  current  density  was  observed  to  spread  out 
over  a  comparatively  large  area,  the  current  transferred  In 
fnc  vicinity  of  the  curved  meniscus  is  small  compared  to 
the  total  current.  This  suggests  that  the  precise  geometry 
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of  the  meniscus  Is  not  Important"^ ^  ^  ,  Consequently,  In  our 
model  the  thin  film  is  assumed  to  change  its  thickness 
abruptly,  yielding  a  hypothetical  film  (see  Figure  l). 

Mass  transport  considerations:  Mathematically,  the 
model  is  one  dimensional.  No  variations  in  concentration, 
c,  or  current  density,  1,  across  the  width  of  the  film  is 
considered.  Oxygen  is  assumed  to  reach  the  electrode  by 
diffusion  in  the  negative  x  direction  through  the  upper 
film  (see  Figure  l).  In  the  lower  (bulk)  electrolyte  region 
oxygen  is  assumed  to  move  only  in  the  negative  y  direction. 
The  consumption  of  oxygen  through  reaction  at  the  electrode 
is  treated  as  a  homogeneous  sink  term  (see  equation  (10)). 

The  system  is  assumed  to  be  at  steady  state.  The  over¬ 
all  electrode  reaction  is 

The  rate  of  reaction  at  the  electrode  is  proportional  to 
the  transfer  current  density,  J.  J  Is  taken  as  positive 
when  the  reaction  proceeds  from  left  to  right  as  written. 

In  the  thin  film,  0^  migrates  tow  .rd  the  electrode  in  the 
negative  x  direction,  0H~  ions  move  down  the  film,  and  the 
K''"  ions  remain  stationary  ?t  steadv  «tate.  The  water  can 
reach  a  reaction  site  by  two  mechanisms  or  combinations 
thereof.  It  can  diffuse  up  into  the  film  countercurrently 
to  the  0H~  ions  or  condense  from  the  gas  phase  onto  the 
film  surface  and  diffuse  co-current ly  with  the  oxygen. 
Convection  is  assumed  to  be  negligible. 


Zn^O  4e  =  40H  . 
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Equations  for  electrolytic  transport;  Algebraically, 
the  Ionic  fluxes  and  water  flux  can  be  written  as 


N, 


0, 


N 


- 


V*  =  ^  (N  +  N 

''r£  + 


By  electroneutrallt-: 


<=+  =  <=. 


“o)  = 


c . 


(1) 


(2) 


Q  is  equal  to  unity  If  the  water  is  supplied  wholly  by 
diffusion  up  the  film,  and  It  Is  equal  to  zero  If  supplied 
totally  by  condensation  from  the  gas  phase.  It  seems  reason¬ 
able  that  one  mechanism  or  the  other  will  be  controlling 
so  fnat  Q  will  be  assumed  to  equal  either  one  or  zero.  The 
model  Is  not  suitable  for  fractional  values  of  Q. 

The  equations  representing  the  flow  of  current  and  raove- 

+  -  14 

ment  of  the  K  and  OH  Ions  have  been  discussed  elsewhere 

Incorporating  the  concentration  dependence  of  the  transport 

properties,  we  obtain 
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Values  of  A,  B,  and  G  for  KOH  are 

A  =  28.89xlC"®cm,  B  =  8.63x10"®  cm,  G  =  9.63x10"®  cm. 
Equations  for  oxygen  transport:  In  the  thin  film, 
oxygen  diffusion  can  be  represented  by 


PD  (pc°-o  ) 

j  =  — ■  ^  ^  ,  (y  ^  0). 


(6) 

(7) 

(8) 


(9) 


The  transfer  current  density,  J,  is  equal  to  the  oxygen  flux 

times  4F.  It  is  assumed  that  Henry's  law  holds  and  that 

p  is  the  partial  pressure  of  oxygen  in  atmospheres.  is 

the  diffusion  coefficient  of  oxygen  at  the  concentration 

of  KOH  existing  at  the  point  under  consideration  in  the 

film.  The  numerical  calculations  were  easily  adapted  to  allow 

for  the  variation  of  as  a  function  of  the  KOH  uoncentra- 

* 

tion  at  different  positions  in  the  film. 

In  the  lower  or  bulk  electrolyte  region  the  following 
equation  applies  to  steady  state  oxygen  transport. 

d^c 

D  ^  ,  (y  i  0).  (10) 

d  y 


♦ 

For  the  diffusion  coefficient  of  oxygen  in  KOH  solutions, 
see  report  on  Project  V,  pages  173-174. 
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Equations  relating  current  and  overpotentlal* :  A  rela¬ 
tionship  between  local  overpotentlal  and  local  transfer 
current  density,  J,  is  still  needed.  A  review  of  the  litera- 
ture  on  the  kinetics  of  the  oxygen  electrode  *  indicates 

that  a  generally  accepted  kinetic  description  for  the  oxy- 

1 7  18  1 9 

gen  electrode  is  not  yet  available  *  *  .  Based  on  work 

PC)  PI  PP 

by  Krasilshchlkov,  Nefodova,  Bellna,  and  Andreeva  *  *  * 

23,24,25  work  by  Bogotskli  and  Yablokova^^,  a  sequence 

of  elementary  steps  emerges  which  seems  to  be  satisfactory 

for  use  with  Ag,  Au,  and  Hg  electrodes.  Prom  this  reaction 

16 

sequence,  summarized  by  Vetter  ,  the  following  kinetic  ex¬ 
pression  can  be  developed^^ 


J  =  ktO^l^IOirj'lHgO]'^ 


r  4(a-l)F'n 

L 


4aFT) 


(11) 


where  k  is  a  constant,  square  brackets  indicate  the  acti¬ 
vity  of  the  enclosed  species,  and  T)  is  the  overpotentlal 
defined  as  the  potential  of  the  working  electrode  compared 
to  a  reversible  oxygen  electrode  located  just  outside  the 
double  layer.  For  a  silver  electrode  in  alkaline  solution 

a  0.625,  X  =  0.5,  y  =  0.5.  (l2) 

Work  on  carbon  electrodes'^  *  ’  *  also  seems  to  sup¬ 

port  the  above  mentioned  reaction  sequence.  However,  infor¬ 
mation  regarding  the  kinetics  of  oxygen  reduction  on  nickel 
is  quite  contradictorj' .  The  surface  oxides  of  nickel  form 
at  about  the  same  potential  at  which  oxygen  reduction  takes 
place  and  a  generally  applicable  reaction  sequence  or  kine- 

tic  rate  expression  does  not  appear  available  * 

*  A  more  detailed  discussion  and  references  are  given  in 
an  Appendix  (page  56  ) . 
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Because  of  the  lack  of  availability  of  a  suitable  kinetic 
expression  for  the  reduction  of  oxygen  on  nickel,  no  attempt 
is  made  to  apply  the  present  mathematical  model  to  the  nickel 
electrode . 

Two  relations  are  needed  in  addition  to  the  equations 
already  mentioned.  Since  charge  must  be  conserved,  the 
transfer  current  density,  J,  and  the  current  density  in  the 
solution,  1,  whlc’*.  is  assumed  to  vary  only  in  the  y  direc¬ 
tion,  are  related  by 

J  =  -  B  (y  S  0);  J  =  -  A  (y  s  0).  (13) 

It  is  convenient  to  replace  the  potential  of  the  reversible 
oxygen  electrode  by  that  of  a  Hg/HgO  reference  electrode  with 
respect  to  the  working  electrode. 

Activities  used  in  kinetic  equation;  The  activities 
of  the  reacting  species  are  assuaied  to  be  represented  by 
the  following: 

[Og]  is  assumed  to  be  the  concentration  of  oxygen  with 
unit  activity  being  the  concentration  of  oxygen  in  equili¬ 
brium  with  gaseous  oxygen  at  a  partial  pressure  of  one  at¬ 
mosphere.  In  the  numerical  work  it  was  possible  to  allow 
the  solubility  of  the  oxygen  to  depend  on  the  concentration 
of  KOH.* 

As  reported  under  Project  V  (page  ),  the  solubility  of 
oxygen  at  25® C  in  aqueous  KOH  as  determined  in  this  labora¬ 
tory  may  be  represented  by 

c^  =  (1.18  X  10'"^mol/liter)  exp (-0.405  c), 
where  c  =  concentration  of  KOH,  moles/liter. 
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[0H"j  is  assumed  to  be  the  same  as  the  mean  ionic  acti¬ 
vity,  a,  of  the  KOH.  The  activity  coefficients  were  taken 

2  9 

from  Robinson  and  Stokes  .  The  activity  at  a  KOH  concen¬ 
tration  of  3.866  molar  is  taken  to  be  unity. 

[c^]  is  assumed  to  be  the  concentration  of  free  water. 
This  is  determined  by  subtracting  from  the  total  water  con¬ 
centration  the  concentration  of  water  that  is  associated 

15 

with  the  KOH  ions  .  Unit  activity  is  assumed  to  be  the 
activity  of  water  in  3.866  molar  KuH  solution. 

Numerical  Solution 

There  are  5  dependent  variables,  c,  c^,  J,  Tj,  and  i. 
These  are  related  by  equations  (3),  (4),  (9)  or  (10),  (11), 
and  (13).  Equation  (9)  is  used  in  the  upper  film  and  equa¬ 
tion  (10)  in  the  lower  region.  Equation  (14)  allows  one  to 
convert  from  fi  to  when  desirable.  Since  the  equations  are 
slightly  different  for  the  upper  and  lower  sections,  the 
two  sections  are  solved  for  separately  and  their  boundary 
conditions  matched,  l.e.,  the  potential  and  current  are  con¬ 
tinuous  in  the  elecT^rolyte . 

The  boundary  conditions  at  the  top  of  the  upper  film 

are : 

y  — ►  +  00,  T1  — ►-  0  and  1  — 0 . 

The  boundary  conditions  at  the  bottom  of  electrode  are: 

at  y  =  bottom  of  electrode,  T]  =  Tl  or  i  =  i  . 

0  0 

The  equations  are  nonlinear,  Ist  and  End  order  differ¬ 
ential  equations.  They  were  transformed  into  finite  dif¬ 
ference  form  and  solved  numerically  using  an  IBM  7090 
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digital  computer.  In  a  numerical  sense,  it  was  possible 
to  solve  the  equations  for  the  thin  film  region  ezactly. 

That  is,  no  Iterations  were  necessary.  This  via.z  accom¬ 
plished  through  an  appropriate  change  of  variables.  In  the 

lower  section,  it  was  necessary  to  linearize  the  equations, 

15 

put  them  in  matrix  form,  and  solve  the  matrix  .  The  non- 
linearities  were  mild  and  only  three  or  four  Iterations 
were  necessary.  The  method  appeared  to  be  stable  and  to 
converge  rapidly  to  the  correct  answer. 

Theoretical  Results  and  Comparison  with  Experimental 
Observations . 

Explanation  of  figures;  The  experimental  results  are 
summarized  in  Figures  2  through  6.  A  comparison  between  the 
calculated  variation  of  current  above  a  fixed  point  on  the 
electrode,  AL,  and  the  experimental  observations  corres¬ 
ponding  to  similar  conditions  is  given  in  Figure  2.  Ex¬ 
perimental  observations  at  48.4,  15.4,  and  1.48  M.A/cm  were 
presented  in  Part  a,  Flg.4-G.  The  agreement  of  the  theoreti¬ 
cal  results  with  the  experimental  results  for  the  silver 
electrode  is  equally  as  good  as  the  best  curves  drawn  in  by 
hand  to  represent  the  data.  No  experimental  data  are  avail¬ 
able  for  the  voltage  variation  in  the  film. 

Figure  3  shows  the  theoret’'  ..,1  variation  of  potential 
and  the  percent  current  above  a  fixed  point  on  the  electrode 
as  a  function  of  applied  current.  Figures  4  and  5  present 
the  current  density  variation  along  the  electrode.  The 
curves  in  Figure  3  are  the  Integrals  of  the  curves  in 
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Figures  4  and  5  for  corresponding  total  applied  currents . 

Figure  6  shows  examples  of  concentration  variation  in  the  film 
as  a  function  of  electrode  height.  A  comparison  Is  made 
between  the  concentration  variations  for  both  modes  of  water 
transport  to  the  reaction  sites. 

Current  distribution  and  the  square  meniscus  assumption: 

AL  Is  the  Independent  variable  used  in  the  experimental 
work  which  Indicates  position  on  the  electrode  relative  to  the 
bottom  of  the  Intrinsic  or  visible  meniscus.  The  variable  y 
Is  used  In  the  theoretical  work  to  In^ -cate  position  along 
the  electrode  relative  to  the  sudden,  right  angle  widening  of 
the  thin  film.  In  order  to  make  direct  comparison  between  the 
two  approaches,  y  and  AL  must  bv  properly  relate 

The  distance  from  the  top  of  the  Intrinsic  meniscus  over 
which  the  current  density  remains  appreciable  is  to  be  nr  ^d. 
Even  at  the  lowest  currents  observed  (1.45  M«A/cm),  the  area 
lying  under  the  level  of  the  bottom  of  the  Intrinsic  meniscus 
contributes  less  than  a  third  to  the  total  current.  At  the 
highest  currents  observed  (154  pA/cm)  this  fraction  falls  to 
belov/  a  hundreth.  The  discontinuity  of  the  current  density 
distribution  curves  at  AL  of  1  ram  (see  Figures  4  and  5)  re¬ 
sults  from  the  square  shaped  meniscus  used  in  the  model.  In 
reality,  the  transition  from  the  charge  transfer  and  conductance 
coi,*" rolled  region  In  the  film  to  the  oxygen  mass  transport 
controlled  area  below  the  top  of  the  Intrinsic  meniscus  and 
then  drops  off  rapidly  as  the  film  begins  to  thicken  from  ap¬ 
proximately  1  micron  to  a  thickness  In  the  order  of  1  mm  at 
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the  bottom  of  the  intrinsic  meniscus-  The  refion  of  hdgh  cur¬ 
rent  density  at  the  bottom  of  the  thin  film  contributes  a  large 
portion  of  the  total  current.  However,  the  reaction  remains 
appreciable  for  several  millimeters  above  the  intrinsic  meniscus. 
Due  to  the  large  area  available  in  the  upper  film,  the  contri¬ 
bution  to  the  total  current  from  that  region  is  appreciable. 

In  matching  the  theoretical  and  experimental  current 
distribution  curves,  it  was  found  that  y  =  0  corresponded  to 
a  AL  of  1  mm.  This  implies  that  at  one  millimeter  above  the 
bottom  of  the  intrinsic  meniscus  the  film  thickness  is  suffi¬ 
cient  to  cause  oxygen  transport  to  become  controlling. 

A  theoretical  treatment  using  a  curved  meniscus  and  a 
tapered  film  is  possible.  However,  it  appears  this  would 
greetly  Increase  the  numerical  complexity  due  to  a  second  geo¬ 
metric  dimension.  Such  a  modification  is  not  expected  to  yield 
significant  further  information  or  insight.  It  is  established 
by  experimental  evidence^  and  a  reasonably  realistic,  theoreti¬ 
cal  model  that  an  overwnelming  proportion  of  the  current  ori¬ 
ginates  well  above  the  bottom  of  the  intrinsic  meniscus.  The 
changes  in  the  pattern  of  the  current  density  distribution 
that  can  be  expected  if  one  substitutes  a  smoothly  curving 
meniscus  for  the  square  meniscus  are  not  expected  to  be 
appreciable . 

Effect  of  exchange  current  density,  1^:  Of  the  parameters 
which  appear  in  the  numerical  calculations,  many  can  be  calcu¬ 
lated  from  Independent  experimental  data.  Where  values  are 
not  available,  reasonable  estimates  have  to  be  made. 
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The  parameter  known  with  the  least  certainty  is  k  In 
equation  (ll).  When  multiplied  by  the  concentration  term 
raised  to  the  appropriate  power  in  equation  (ll),  k  can  be 
identified  vfith  an  exchange  current  density,  1^.  A  value  for 
i^  of  1.76  X  10”^  A/cra^  was  found  to  "fit"  the  data  the  best. 
However,  this  could  be  varied  by  plus  or  minus  three  orders 
of  magnitude  with  little  effect  on  the  position  or  shape  of 
the  theoretical  curve.  Although  literature  reports  on  measure¬ 
ments  of  i^  for  the  oxygen  electrode  vary  widely,  a  comparl- 

■P5 

son  to  the  Tafel  lines  shown  by  Bellna  and  Krasilshchlkov 
indicate  that  the  value  used  is  not  unreasonable. 

Effect  of  the  film  thickness  parameter,  6;  The  theoretical 
results  were  very  sensitive  to  changes  in  the  film  thickness 
parameter,  6.  It  should  be  remembered  that  in  reality  the 
film  tapers  while  the  model  assumed  5  to  be  constant.  Since 
the  current  extends  further  up  the  film  at  low  currents,  it 
is  to  be  expected  that  the  effective  average  film  thickness 
will  decrease  with  decreasing  total  current.*  The  following 
values  of  6  were  used  for  the  indicated  total  applied  currents. 


I 

1.48  iiA/cm 

4.84  MA/cm 

15 .4iiA/cra 

48.4  pA/cm 

154  iiA/cm 

6 

0.33  micron 

0.43  micron 

0.75  micron 

1.50  micron 

2  .25  micron 

The  values  of  6  used  in  the  numerical  solutions  are  in  reason¬ 
able  agreement  with  the  average  film  thicknesses  obtalmd  from 
Miiller’s^-^"^^  data. 


* 

It  should  be  remembered  that  at  low  total  current  the  rela¬ 
tive  current  distribution  is  more  spread  out.  Thus,  the 
averaging  of  the  film  thickness  includes  more  of  the  thinner, 
upper  parts  of  film. 


. . . . . . . . . . . 
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As  confirmed  by  our  models,  experimental  observations 
have  shovm  tnat  the  current  density  dlstrlbutjon  varies  a 
rreat  deal  while  the  film  is  draining,  l.e.,  while  the  film 
thickness  is  dlmlshlng  with  time.  In  tlilnnlng  films,  the 
overall  overpotential  Increases  and  the  reaction  spreads  less 
and  less  into  the  fixm. 

Influence  of  variable  transport  paramete rs  and  water  trans¬ 
port  :  In  the  present  work,  current  distributions  were 
calculated  both  with  constant  and  variable  transport  properties. 
Although  the  assumption  of  constant  (average)  properties  did 
not  yield  significantly  different  distributions,  this  finding 
should  be  regarded  as  valid  only  under  the  specific  experi¬ 
mental  conditions  employed.  By  providing  water  to  the  reac¬ 
tion  sites  through  the  film  from  the  vapor  phase,  the  steady 
state  concentration  gradient  in  the  film  is  quite  mild  compared 
to  the  case  when  water  transport  occurs  by  diffusion  from  the 
bulk  electrolyte  in  the  liquid  phase  (see  Figure  6). 

The  high  concentration  gradients  which  are  predicted  when 
water  is  assum.€;d  to  be  transported  only  within  the  film  cause 
the  KOH  concentration  to  exceed  the  saturation  limit.  It  is 
suggested  that  in  actual  porous  gas  electrodes  water  is  supplied 
to  the  reaction  site,  at  least  Ir  part,  by  condensation  from 
the  vapor  phase . 

Charge  trai-sfer  kinetics;  The  shape  and  position  of  the 
theoretical  current  distribution  curves  also  depends  on  the 
reaction  orders,  a,  X,  and  y.  These  appear  as  the  exponents 

and  Og  activities,  respectively,  in  equation 


of  the  OH  ,  HgO 
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(ll).  The  values  of  a,  X,  y  given  in  equation  (12)  were  based 
on  the  reaction  sequence,  rate  controlling  step,  +  e  = 

and  value  of  the  transfer  coefficient  recommended  by 
Kraoilshchikov'^® .  There  are,  of  course,  a  very  large  number 
of  reaction  sequences  and  ra-e  controlling  steps  which  might  be 
chosen.  Calculations  involving  another late  controlling  step 
(Og  t  e”  =  Og)  yielded  a  good  fit  of  the  experimental  data 
only  when  unrealistic  values  of  L  were  employed. 

Oxygen  transport:  Except  for  the  region  below  the  menis¬ 
cus,  oxygen  transport  should  not  be  a  limiting  factor.  The 
numerical  results  show  that  throughout  and  across  the  entire 

length  of  the  film,  the  activity  of  oxygen  is  very  close  to 

** 

unity  ,  i.e.,  nearly  identical  to  the  activity  at  the  electro¬ 
lyte-gas  interface.  However,  once  one  moves  below  the  meniscus 
where  the  diffusion  path  for  the  oxygen  becomes  large,  the 
oxygen  activity  drops  quickly  resulting  in  the  rapid  decrease 
in  current  density  below  the  top  of  the  Intrins'' c  meniscus. 

Rate  controlling  process:  Since  our  model  shows  that  the 
oxygen  activity  does  iCt  vary  noticeably  in  the  thin  film, 
it  follows  Lhat  the  penetration  of  the  reaction  into  the  film 
depends  on  a  balance  between  the  charge  transfer  overpotential 
and  resistance  drop  in  the  electrolyte.  As  stated  in  Part  a, 
the  qualitative  differences  between  behavior  of  the  n-’.ckel  and 
silver  electrodes  also  suggest  this  to  be  the  case. 

*  / 

Note  that  the  parameter  a  is  relaceo.  to  (but  not  the  same 

as)  the  transfer  coef f i  c-i  ent; . 

It  bhoulc’.  oe  remembered  that  although  ohe  activity  remains 
constant,  the  concentration  of  o>^/gen  varies  along  the 
length  of  t.he  film  due  to  the  decreasing  solubility  with 
increasing  KOH  concentration. 
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Will  In  his  work  with  the  hydrogen  electrode  on  plati¬ 
num  found  that  t" e  transport  of  hydrogen  gas  to  the  electrode 
surface  through  the  thin  film  was  controlling,  along  with 
ohmic  resistance  drop,  except  at  very  low  currents.  This  is 
not  necessarily  in  contradiction  with  the  results  found  for 
the  oxygen  electrode. 

The  hydrogen  electrode  on  platinum  is  much  more  reversible 
than  is  oxygen  on  silver.  Consequently,  the  concentration 
of  current  in  the  immediate  vicinity  of  the  intrinsic  meniscus 
is  rendered  less  possible  in  the  case  of  oxygen,  and  charge 
transfer  overpotential  along  with  the  ohmic  drop  in  the  film 
can  be  expected  to  control  the  current  distribution. 

Conclusions 

] .  Oxygen  transport  is  not  limiting  in  the  thin  upper 
film.  The  current  density  distribution  in  this  region  arises 
from  a  balance  between  charge  transfer  overpotential  and 
resistance  drop  in  the  film. 

2.  Below  the  intrins'^c  meniscwis,  the  bulk  electrolyte 
region,  the  oxygen  activity  decreases  rapidly  due  to  the  mass 
transport  limitations  of  dissolved  oxygen.  Thus,  the  current 
density  is  small  in  this  region. 

3.  The  high  concentration  gradients  which  result  when 
water  is  not  supplied  from  the  gas  phase  point  to  the  desir¬ 
ability  of  supplying  water  as  vapor  in  the  oxygen  gas. 
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Nomenc latrre 

”  subsidiary  transport  properties  defined  so  as  to  reduce 
the  concentration  dependence  (a). 

-  mean  molar  activity  of  KOH  =  f^c. 

•Z 

-  c_^/v_^  =  c_/v_  -  concentration  of  KOH  (mole/cm  ). 

-  conc3ntration  of  oxygen  (mole/cm  ). 

-  saturation  concen oration  of  oxygen  under  a  partial  pressure 
of  1  atm.  (mole/cm^-atm) . 

-  concentration  of  K"^  (mole/cm^). 

-  concentration  o^*  OH”  (mole/cm  ). 

-  concentration  of  water  (mole/cn  ). 

«  c_j_  +  c_  +  c^  -  total  concentration  (Including  solvent). 

-  diffusion  coefficient  of  oxygen  in  aqueous  KOH  (cm  /sec). 

-  molecular  diffuslor  coefficient  of  KOH  based  on  activity 
gradients  (cra^/secy. 

-  electronic  charge  (coulombs). 

-  symbol  for  an  electron. 

-  mean  molar  activity  coefficient.  (Values  from  reference 
33  were  used.) 

-  Faraday’s  constant  ( coulombs/equl valent ) . 

-  current  density  in  the  solution  (amp/cm^). 

p 

-  exchange  current  density  (amp/cm  ). 

2 

-  transfer  current  density  (amp /cm  ). 

-  kinetic  parameter  related  to  exchange  current  density. 

-  flux  of  k  (me" ^/cm  -sec ) . 


. . . . . . . 
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N_  -  flux  of  OH"  (mole/cm^-sec) . 

p 

-  flux  of  HgO  (mole/cm  -sec), 
p  ~  partial  pressure  of  oxygen  (atm) . 

Q  -  one  when  H2O  is  transferred  in  liquid  phase,  zero  when 
HgO  transferred  in  the  gas  phase 

R  -  universal  gas  constant  ( joule /nole -®K) . 

t^  -  transference  numoer  of  species  i  with  respect  to  molar 
average  velocity. 

T  -  absolute  temperature  (®K) . 

v^  -  molar  average  vdoclty  (  cm/sec). 

X  -  coordinate  perpendicular  to  the  plane  of  the  electrode 
(cm) . 

y  -  coordinate  parallel  to  the  electrode  in  the  vertical 
direction,  the  independent  variable  (cm). 

z.  -  valence  or  charge  number  of  species  i  (  negative  for 
^  anions). 

a  -  kinetic  parameter  related  to  reaction  order  of  oxygen 
and  to  the  transfer  coefficient. 

y  -  kinetic  param'=’^‘'  indicating  reaction  order  of  H2O. 

6  -  film  thicknes  crons). 

A  -  width  of  annular  el.  ctrolyte  region  below  meniscus  (cm). 

T\  -  overpotential  referred  to  a  reversible  oxygen  electrode 
(volts) . 

K  -  conductivity  (mho/cm). 

X  -  kinetic  parameter  Indicating  reaction  order  of  OH", 
p.  -  chemical  potential  of  KOH  (  JoUle/mole).  (RT  Ik  f+  c) . 

p,  -  viscosity  (poise). 

-  number  of  cations  per  molecule  of  electrolyte. 

V  -  number  of  anions  per  molecule  of  electrolyte. 


$  -  overpotential  referred  to  a  reversible  lig/HgO  electrode 

(volts) . 
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APPENDIX 

Notes  on  the  Kinetic  of  the  Oxi’’ger  Electrode 

There  has  been  accumulated  a  body  of  Information  v;hlch 
suggests  that  hydrogen  peroxide  is  an  Intermediate  In  oxygen 
reductlon^^”^^^  34-38^  Although  the  evidence  Is  fairly  conclu¬ 
sive  for  such  electrodes  as  Ag,  Hg,  C,  and  Au,  It  appears  that 
there  may  be  competing  sequences  of  steps  on  such  electrodes  as 
Pt  and  this  point  only  the  former  category 

of  electrodes  appears  amenable  to  a  simple  kinetic  analysis. 

A  review  of  the  above  mentioned  references  indicates  that 
the  assumption  of  the  following  sequence  of  elementary  kinetic 
steps  is  reasonable  on  Ag,  Hg,  c,  and  Au  electrodes. 

(1)  Og  +  M  =  MOg 

(2)  MOg  +  e~  =  MO" 

(3)  MOg  +  HgO  =  MHOg  +  OH" 

(4)  MHOg  +  e~  =  MHO“ 

(5)  MHO"  +  HgO  =  MHgOg  +  OH" 

(6)  MH^Og  +  e"  -  MOH  OH" 

(7)  MOH  i-  e-  =  M  OH" 

The  M  represents  an  adsorption  site,  and  when  associated  with  a 

molecular  symbol  Implies  that  tho  molecule  Is  adsorbed.  It 

Is  possible  for  any  adsorbed  molecule  on  the  right  side  of  an 

equation  to  desorb  before  reacting  In  the  next  step.  The  Og 

and  HOp  are  so  reactive  that  this  Is  unlikely.  The  desorption 

of  HOp  or  HpOg  (depending  on  the  pH)  could  occur. 

20  25 

However,  there  Is  evidence  *  which  indicates  that  on  sliver 
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and  In  alkaline  solution  the  rate  of  desorption  is  smalx  com¬ 
pared  to  reaction  by  steps  (6)  and  (7).  It  should  be  noted  that 
because  the  O'  and  HO.-,  radicals  are  very  reactive,  they  have  not 
been  positive  identified. 

It  appears  that  the  adsorption  step  is  not  rate  controlling 
and  that  either  step  (2)  or  (4)  is  rate  determining.  The  choice 
appears  to  depend  on  pH  and  on  th(,  electrode  material.  For 
silver  in  alkaline  solution  step  (4)  appears  to  be  rate  con¬ 
trolling.  With  these  assumptions,  the  kinetic  expression  can 
be  written  as 

I  (a^-2)Fn  (a4-t2)FTl 

-  e 


<^4+2 


a, 


J  =  k[02l  ^  [HgO]  '^[0H"1 


(!-«/) 


I 


RT 


RT 


where  J  is  the  local  transfer  current  density,  is  the  trans¬ 
fer  coefficient  for  step  (4),  Tj  is  the  local  overpotential, 
and  k  is  a  kinetic  rate  constant.  Brackets  signify  the  acti¬ 
vity  of  the  enclosed  species. 
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Comparison  of  experimentally  observed  current  distribution  vrith  the 
distribution  calculated  using  the  theoretical  model.  Nickel  cathode 
I  =  154  iiA/cm.  Cell  no.  4,  runs  7,  9  and  12.  Numerical  solution 
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Figure  5.  Theoretical  current  density  distribution  relative  to  the  position  of  the 
meniscus.  High  current  range. 
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Figure  6.  Calculated  variation  of  KOH  concentration  in  the  electrolyte  film. 
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PROJECT  II. 

Ionic  Mass  Transport  by 
Combined  Free  and  Forced  Convection 

The  problem  of  hydrodynamic  flow  caused  by  density 
gradients  existing  in  the  vicinity  of  working  electrodes  has 
been  the  subject  of  numerous  experimental  and  theoretical 
studies  (10,17,18,31,32).  The  complex  nature  of  the  role 
of  solution  properties,  geometric  factors  and  current  den¬ 
sities  in  the  case  of  laminar  or  turbulent  free  convection 
at  vertical  electrodes  is  reasonably  well  understood,  and 
satisfactory  quantitative  relations  describing  the  functional 
relationships  between  pertinent  process  variables  are  avail¬ 
able.  Mass  transport  by  free  convection  at  horizontal 
electrodes  has  also  been  treated,  and  although  the  convec- 
tional  mode  is  not  well  understood  in  this  case,  correla¬ 
tion  of  variables  in  dimensionless  form  proposed  by  Fenech 
and  Tobias  (lO)  quite  satisfactorily  predicts  experimental 
behavior.  Free  convection,  because  of  the  high  velocity 
gradients  generated  in  the  boundary  layer  adjacent  to  the 
electrode,  provides  a  very  effective  mode  of  mass  transport. 
In  so-called  unstirred  solution  "equivalent"  mass  transfer 
film  thicknesses  (29)  in  the  order  of  a  few  tenths  of 
millimeters  are  commonly  encountered. 

Examination  of  typical  laboratory  and  industrial  elec¬ 
trolysis  conditions  reveals  that  stirring  (or  bulk  flow) 
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rates  are  often  not  adequate  to  overshadow  the  effects  of 
free  convection. 

The  scientific  literature  appears  to  be  void  of  studies 
aimed  at  the  clarification  of  the  effects  of  combined  free 
and  forced  convection  on  mass  transfer  rates.  A  general 
theoretical  description  of  the  convectional  modes  in  the 
presence  of  free  and  forced  convection  represents  a  formidable 
task  (2,12,24,30).  Rudimentary  experimentation  on  the  analo¬ 
gous  heat  transfer  case  in  the  1930 's  (3,14,26)  did  not  lead 
to  quantitative  description  of  the  phenomena  observed.  At¬ 
tempts  at  correlating  experimental  results  have  been  unsuccess¬ 
ful  (1,20). 

For  the  case  of  pure  laminar  forced  convection  experimental 
and  theoretical  investigations  have  provided  a  reasonable 
basis  for  prediction  of  ionic  mass  transport  rates.  Lin  et  al 
(22)  found  that  in  the  cell  of  annular  geometry,  limiting 
currents  on  the  outer  electrode  are  well  represented  by 

Nu^  =  l-62(Re.Sc.d/L)^/^  (l) 

This  equation  is  based  on  the  Leveque  approximation  and  is 
applicable  when  the  velocity  profile  is  fully  developed  in 
the  electrolyte  approaching  the  leading  edge  of  the  electrode. 

It  should  be  noted  that  the  equation  applies  only  when  the 
concentration  profile  is  still  not  fully  developed  at  the 
trailing  edge  of  the  electrode.  Except  for  e..tremely  low 
flow  rates,  this  condition  would  be  met  in  typical  situations 
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Involving,  lectrode  reactions.  On  the  basis  of  Norris 
and  Streid's  work  on  heat  transfer  In  rectangular  flat 
ducts  (25),  one  would  expect  the  numerical  constant  to 
change  from  1*62  to  1*85  for  the  case  of  Ionic  mass  trans¬ 
port  in  laminar  flow  to  planar  electrodes  in  rectangular 
channels : 

Nu^  =  l*85(Re.Sc.d/L)^/^  (2) 

The  present  investigation  was  undertaken  to  gain  some 
insight  into  the  nature  of  mass  transport  under  combined 
free  and  forced  con\ectlon  conditions.  As  an  experimental 
model,  a  horizontal  channel  type  electrolytic  cell  was  chosen. 
In  such  a  cell,  when  the  electrodes  at  which  ions  are  con¬ 
sumed  (e.g.,  a  cathode  at  which  a  metal  is  deposited)  faces 
upward,  in  the  absence  of  bulk  flow  the  electrolyte  layer 
adjacent  to  the  electrode  becomes  unstable,  and  except  for 
the  case  of  extremely  dilute  solutions .  turbulent  edaies  are 
formed.  The  distribution  of  eddies  over  the  electrode  sur¬ 
face  is  statistically  uniform,  and  as  i  consequence  the 
effective  mass  transfer  boundary  layer  thickness  over  ex¬ 
tended  horizontal  electrode  surfaces  does  not  vary  with  posi¬ 
tion  (ll).  In  a  typical  case,  the  turbulent  eddy  motion 
extends  only  a  few  m'.llimeters  from  the  surface;  above  this 
distance  any  opposing  llmiiting  surface  has  no  noticeable 
effect  on  the  transport  rates  (ll). 

Llmliing  currents  in  pure  free  convection  extended  over 
horizontal  electrode  surfaces  were  successfully  correlated 
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by  Penech  and  Tobias  (lO),  who  found  that 

Nu|^  -  0*19(Gr  3Sc)V2  ^2^ 

L 

0 

represents  the  experimental  values  In  the  range  of  10  <  Gr. 

Sc  <  1*4  X  10^^  with  a  correlation  coefficient  of  0*094. 

If  one  Introduces  a  steady  state  bulk  flow  Into  a  hori¬ 
zontal  electrolyzer,  the  viscous  shear  near  the  walls  Inter- 
i'eres  with  the  development  of  the  turbulent  free  convection 
motion  In  the  boundary  layer.  At  the  leading  electrode  edge 
the  mass  transfer  boundary  layer  thickness  Is  zero.,  and  only 
at  some  distance  downstream  will  appreciable  boundary  layer 
thickness  develop.  Eventually  this  thickness  exceeds  the 
stability  limit,  and  motion  due  to  the  gravitational  field 
will  begin.  Proceeding  further  downstream,  the  bormdary  layer 
reflecting  the  contributions  of  both  forced  and  free  convec¬ 
tion  reaches  a  steady  state  conditions.  This  qualitative 
Image  of  the  transport  process  Indicates  that,  unlike  In  the 
case  of  pure  free  convection,  current  densities  will  change 
In  the  direction  of  flow  -  probably  from  a  high  Initial  value 
at  the  leading  edge  to  some  steady  state  value  downstream. 

EXPERIMENTS 

Limiting  currents  were  measured  In  a  cell  that  was  de¬ 
signed  to  yield  values  of  average  and  local  mass  transfer 
rates  over  a  broad  range  of  the  pertinent  variables.  The 
cell  had  separate  parallel  flow  channels  for  the  anolyte  and 
catholyte  as  shown  In  Isometric  section  In  Figure  1.  The 
channels  were  separated  by  a  luclte  plate,  except  between  the 
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anode  and  cathode,  where  a  porous  ceramic  separator  was  used. 
This  separator  permitted  free  passage  of  ions  but  Insured 
that  there  would  be  no  disturbance  at  the  cathode  by  the  hy¬ 
drodynamic  effects  from  the  dissolution  process  occurring 
at  the  anode.  A  saw-toothed  profile  v/as  cut  into  the  anode 
to  minimize  the  possibility  of  reaching  the  saturation  con¬ 
ditions  at  the  anode  during  experiments  with  high  concentra¬ 
tions  of  copper.  This  possibility  was  further  decreased  by 
the  dual  set  of  channels  which  made  possible  a  higher  flow- 
rate  past  the  anode  than  past  the  cathode.  The  60  cm  entrance 
length  necessary  to  always  Insure  a  fully  developed  velocity 
profile  was  calculated  by  the  methods  proposed  by  Han  (15). 

The  flow  channel  dimensions  could  be  altered  by  placing 
luclte  Inserts  Into  the  channel. 

The  cathode  fit  Into  a  holder  that  was  clamped  to  the 
slae  of  the  cathode  flow  channel,  so  that  the  electrode  sur¬ 
face  was  coplanar  with  the  rest  of  the  duct  wall.  This 
holder  could  accommodate  two  types  of  electrodes:  continuous 
electrodes  of  various  lengths  for  use  In  experiments  designed 
to  yield  average  mass  transfer  coefficients  and  a  sectioned 
electrode  (11,23)  for  use  In  obtaining  local  mass  transfer 
rates . 

In  the  latter  case,  the  general  technique  of  cutting 
the  electrode  Into  segments,  insulating  the  segments  from  one 
another,  and  then  meas' _ ig  the  current  flowing  to  each  seg¬ 

ment.  Within  the  cathode  holder  were  the  electrode  contacts. 
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the  low  resistance  shunts  for  measuring  segments  currents, 
and  a  switching  circuit  for  calibrating  the  shunt  resistances. 
Starting  from  the  lesdirg  edge  the  segment  lengths  were  0*5, 
0*5,  1,  2,  3,  3,  5,  5,  8,  and  8  cm  for  a  total  electrode 
length  of  36  cm.  The  segments  were  insulated  from  one  another 
by  0*004  Inc.  polyethylene  sheeting  and  electrically  contacted 
from  the  rear  by  tv/lst  screws.  The  entire  assembly  was 
fitted  with  la  .^e  wheels  and  mounted  on  rollers  so  that  the 
cell  could  be  oriented  with  the  cathode  facing  either  upward 
or  downward.  This  feature  was  found  to  be  very  useful,  as 
will  be  shown  later. 

For  the  segmented  e],ectrode  to  truly  represent  the  cor¬ 
responding  continuous  electrode,  it  was  Important  that  the 
segments  all  be  coplanar  and  equlpotential.  The  first  condi¬ 
tion  was  met  by  polishing  the  segments  to  a  coplanar  surface 
in  a  special  Jig  that  held  the  segments  in  the  same  position 
as  they  assumed  in  the  cathode  holder.  The  second  requirement 
could  not  be  exactly  met.  Ing  the  electrodiposltion  of 
copper  fi’om  the  acid  solution,  the  overpotential  plateau  at 
limiting  current  reaches  approxl.mately  half  a  volt  before 
hydrogen  evolution  begins.  Limiting  current  measurements 
are  made  on  the  upper  end  of  the  plateau,  so  the  electrode 
overpotential  is  usually  in  the  order  of  400  to  500  mv.  It 
was  arbitrarily  decided  that  a  satisfactory  approximation 
to  an  equlpotential  electrode  would  be  one  in  which  the  poten¬ 
tial  of  any  two  segments  did  not  differ  by  more  than  two 
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millivolts.  With  this  approximation,  measurements  of  the 
individual  section  currents  was  straightforward.  Very  low 
resistance  shunts  were  introduced  between  a  busbar  and  the 
electrode  segments.  The  potential  drops  across  them  were 
selected  sequentially  amplified  vilth  a  low  noise  D.C.  ampl_- 
fler.  aud  then  read  out  on  a  multipoint  printing  recorder. 

This  circuit  is  shown  schematically  in  Figure  2 .  Contact 
resistances  between  the  busbar  and  the  electrode  segments 
v/ere  kept  to  a  minimum  by  silverplating  all  the  mating  parts. 
Limiting  current  curves  v;ere  plotted  continuously  on  an 
X  -  y  recorder  while  current  was  passed  through  the  cell  by 
a  programmed  current  source  which  provided  a  preset  linear 
increase  of  current  with  time.  This  device  is  described  in 
detail  elsewhere  (16).  Electrolyte  was  continuously  pumped 
through  both  flow  channels  and  its  flowrate  was  measured 
with  calibrated  rotameters. 

PROCEDURE 

The  following  consecutive  steps  were  made  during  the 
performace  of  each  experiment: 

1.  The  cathode  was  polished  with  number  600-A  silicon 
carbide  paper  until  all  trace  of  the  previous  deposit  disap¬ 
peared.  Surface  was  dusted  with  lint-free  paper  and  degreased 
with  acetone. 

2 .  Cathode  v/as  mounted  In  the  cathode  holder  and  the 


holder  clamped  to  tlie  body  of  the  cell.  Electrodes  were 
connected  to  the  current  source  and  measuring  circuitry. 
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3.  The  pump  was  turned  on  and  the  flow  channels  were 
filled  with  electrolyte,  carefully  expelling  all  air  bubbles. 
Then  the  desired  flowrate  was  set,  current  rate  selected, 
and  the  measuring  circuits  put  into  operation. 

4.  The  current  was  turned  on.  When  limiting  current 
was  reached,  the  following  data  were  recorded:  limiting 
current,  flowrate,  electrode  length,  channel  cross  section, 
electrolyte  used,  measuring  circuits  used,  and  electrolyte 
temperature . 

5.  Following  eacn  run,  the  cell  was  disassembled  and 
preparations  were  made  for  the  next  experiment. 

6.  The  bulk  concentrations  of  species  in  the  electro- 
-.yte  were  determined  by  standard  techniques  of  analysis  (8). 

7.  The  scheme  developed  by  Fenech  and  Tobias  (10)  was 
used  throughout  to  calculate  the  physical  properties  and  the 
dimensionless  groups. 

Table  1.  Range  of  variables  covered. 

The  ranges  of  variables  covered  in  this  investigation  were 
as  follows: 

Electrode  length (cm) 

Channel  cross  section  (cm) 

Viscosity  (cps) 

Density  (g/cm^) 

Velocity  (cm/sec) 

Concentration  of  Cu(m) 

Schmidt  number 
Reynolds  number 
Grashof  number 
Electrode  Length/Duct  Dlam. 


0*2,  0-5,  1,  2,  5,  10,  20, 

50 ,  36 

2  X  10,  1  X  10,  1x8,  1x5 
1-12 
1-1-1-3 
0-35 

0*01-0*45 
600-12,000 
75-7000 
105-lol2 
0 '05-20 
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QUALITATIVE  OBSERVATIONS 

From  the  appearance  of  the  electrode  deposits  obtained 
on  continuous  cathodes,  it  was  evident  that  there  was  a 
considerable  variation  in  current  density  in  both  the  direc¬ 
tion  of  bulk  flow  and  In  the  transverse  direction.  Variation 
of  current  density  in  the  direction  of  flow  was,  of  course, 
predictable  from  forced  convection  theory,  but  It  appeared 
that  the  magnitude  and  distribution  of  the  nonuniformity 
was  quite  different  from  that  expected.  A  typical  deposit 
Is  shown  In  Figure  3.  To  explain  the  appearance  of  the 
streaked  deposit  It  is  suggested  that  a  secondary  convection 
effect  was  present  which  took  the  form  of  "roll  cells'*  with 
a  velocity  profile  similar  to  that  shown  in  Figure  4  The 
secondary  convection  results  from  an  instability  which  is 
caused  by  the  density  gradient  across  the  mass  transfer 
boundary  layer.  Such  an  explanation  finds  support  In  the 
fluid  dynamics  and  heat  transfer  literature  (3,4,13,30). 

It  should  be  noted  that  the  scale  of  periodicity  of  the 
electrode  deposit  was  In  the  order  of  half  a  millimeter 
after  the  roll  cells  have  grown  to  a  fully  developed  array 
across  the  electrode  surface. 

On  the  basis  of  this  observation,  the  significant  length 
in  the  Rayleigh  number  (equal  to  the  product  of  the  Grashof 
and  Schmidt  numbers)  which  determines  the  stability  of  the 
system  should  be  the  distance  across  which  the  density 
difference  is  manifest.  In  the  system  under  consideration 
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here .  ^hls  distance  Is  the  mass  transfer  boundary  layer 
talckiie  i  6.  For  flow  between  flat  plates  with  a  fully 
developed  velocity  profile  at  the  leading  edge  of  the  trans¬ 
fer  section,  this  is  given  by  (2l) 

6  =  A(d^ .x/Re.Sc)^/^  (4) 

Experiments  were  performed  on  a  30  cm  electrode  to 
see  how  the  onset  of  this  secondary  convection  pattern  was 
influenced  by  the  bulk  velocity  of  the  electrolyte.  From 
eqn.  (4)  it  can  be  seen  that  if  the  physical  and  geometric 
parameters  of  a  system  are  held  constant,  6  varies  as  (x  V 
(x/V)^/^.  Therefore  one  might  expect  that  the  onset  of 
instability  will  occur  at  positions  on  the  electrode  which 
correspond  to  the  szune  value  of  the  Rayleigh  number.  Since 
the  Rayleigh  number  varies  as  the  cube  cf  6,  at  some  con¬ 
stant  value  of  6  the  position  at  v/hich  the  vorticc'  begin 
should  be  linear  in  the  bulk  velocity.  That  this  is  actually 
the  case  is  shown  in  Figure  5, 

As  might  be  expected,  the  presence  of  roll  cells  has  a 
pronounced  effect  on  the  local  mass  transfer  process.  In 
the  direction  normal  to  the  bulk  flow,  variations  in  deposit 
thickness  were  measured  by  two  optical  methods.  The  first 
were  made  on  a  microscope  equipped  with  a  movable  stage  that 
could  be  located  in  two  directions  with  an  accuracy  of  one 
micron.  The  barrel  could  be  positioned  with  the  same  accu¬ 
racy.  Using  250  x  magnification,  the  depth  of  focus  was 
reduced  to  about  one  micron,  and  a  two  dimensional  relief  of 
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the  surface  normal  to  the  direction  of  flow  was  made  before 
and  after  deposition.  The  thickness  of  the  copper  deposit 
was  then  obtained  with  an  accuracy  of  ±2  microns  by  taking 
the  difference  of  the  two  profiles.  Representative  depo¬ 
sition  profiles  are  shown  in  Figure  6.  Within  the  accuracy 
of  the  optical  technique  employed,  the  measured  average 
thickness  of  the  deposits  agreed  with  that  calculated  from 
Faraday's  law.  The  profiles  indicate  that  the  character 
of  the  deposits  caused  by  the  presence  of  roll  cells  is 
essentially  independent  of  flowrate  and  deposition  rate. 

The  physical  nature  of  the  grooves  in  the  deposit  was 
also  asserted  by  photomicrography.  After  prolonged  (up  to 
one  hour)  electrolysis  at  limiting  current,  electrodes  were 
metallographically  sectioned,  transverse  to  the  direction 
of  flow  and  photographed  at  low-power  magnification.  The 
photomicrographs  gave  similar  profiles  to  the  one  illus¬ 
trated  in  Figure  6. 

Local  mass  trajisfer  rates  in  the  direction  of  flow 
were  measured  with  sectioned  cathodes.  Typical  smoothed 
data  obtained  with  this  technique  are  given  in  Figures  7 
through  11.  On  each  figure,  a  series  of  experiments  with 
the  sane  electrolyte  is  shown,  the  independent  variables 
then  being  limited  to  the  velocity  and  electrode  orienta¬ 
tion.  With  the  cathode  facing  upward,  both  convective  pro¬ 
cesses  are  operative,  whereas  the  inverted  position  causes 
stratification  of  the  electrolyte  during  deposition  and 
only  pure  forced  convection  occurs.  The  natural  convective 
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contribution  to  the  local  mass  flux  can  therefore  be  ob¬ 
tained  as  the  difference  between  the  curves  for  the  two 
cell  orientations. 

The  following  observations  and  comments  can  be  made 
regarding  these  curves: 

1.  At  the  leading  edge  and  for  some  distance  down¬ 
stream,  the  current  density  distribution  is  independent  of 
the  electi  de  orientations. 

2.  At  some  position  downstream  from  the  leading  edge 
a  considerable  increase  of  current  density  above  that  for 
pure  forced  convection  occurs  for  the  orientation  resulting 
in  superimposed  free  convection.  The  position  of  current 
rise  was  always  found  to  correspond  to  the  appearance  of 
the  streaks  in  the  deposit. 

3.  Regardless  of  the  velocity  of  flow,  for  a  given 
electrolyte  with  the  cathode  facing  upward  it  appears  that 
the  value  of  current  density  far  down  the  electrode  reaches 
or  approaches  the  same  steady  state  value.  Cn  many  of  the 
experiments  the  electrode  was  not  sufficiently  long  to  al¬ 
low  this  value  to  be  reached. 

4.  For  system.s  with  relatively  large  Rayleigh  num¬ 
bers  (l.e.,  systems  with  low  velocity  [see  eqn.  (4)],  high 
density  differences  [high  copper  concentrations],  and  low 
viscosities),  the  increase  in  current  density  above  that 
for  the  pure  forced  convection  case  occurs  closer  to  the 
leading  edge. 
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5.  Even  at  high  laminar  Reynolds  numbers,  natural 
convection  can  play  a  very  important  role-  If  the  electrode 
is  sufficiently  long,  it  can  be  expected  to  dominate  the 
mass  transfer  process. 

6.  In  many  of  the  experiments  dominated  by  free  convec¬ 
tion,  the  Initial  rise  in  current  density  is  not  maintained, 
but  gradually  declines  to  some  lower  value  farther  long  the 
electrode.  This  behavior  is  particularly  marked  la  systems 
with  a  high  Rayleigh  number,  vihere  the  character  of  electrode 
deposits  after  the  onset  of  roll  cells  gradually  changes 
from  well-defined  streaks  to  no  streaks  at  all.  This  sug¬ 
gests  a  transition  to  a  turbulent  mode  of  free  convection 
near  the  interface  as  the  mass  transfer  boundary  layer 
Increases  (higher  Rayleigh  number).  In  experiments  where 

the  streaks  ret_ined  their  identity  for  the  full  length  of 
the  electrode,  a  drop  in  current  density  is  still  present, 
suggesting  that  the  initial  vortex  motion  is  in  unsteady 
state,  and  some  tlm.e  is  required  before  the  fully  developed 
convectional  mode  is  established. 

From  these  observations,  a  fair  understanding  of  the 
combined  flow  processes  begins  to  emerge.  As  the  electro¬ 
lyte  flows  past  the  electrode  at  limiting  current,  an  in¬ 
creasing  mass  transfer  boundary  layer  thickness  raises  the 
Rayleigh  number  leading  to  an  instability.  Due  to  the 
inertia  of  the  electrolyte  some  time  passes  before  the  in¬ 
stability  causes  a  signiflcan  roll  cell  motion.  During 
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this  incubation  period  the  pure  forced  convection  mode  of 
mass  transfer  continues  making  the  system  Increasingly  un¬ 
stable.  Finally j  at  seme  position  on  the  ^lectrode  down¬ 
stream  from  the  point  of  Initial  instability,  roll  cell 
motion  begins,  but  on  a  short  electrode,  the  unstable  fluid 
may  be  swepv  off  the  electrode  before  the  buoyancy  forces 
can  Influence  the  fluid  dynamic  conditions.  Prior  to  the 
onset  of  the  natural  convective  motion  the  mass  transfer 
process  appears  to  be  the  same  as  for  pure  forced  convec¬ 
tion,  even  though  buoyancy  forces  are  present. 

Onset  of  roll  cells  was  found  to  be  a  function  of  the 
size  and  number  of  the  Imperfections  on  the  electrode  sur¬ 
face,  Just  as  in  the  case  of  initiating  turbulence  (5,9). 

It  was  observed  that  a  small  dent  in  the  surface  or  a  small 
bubble  .sticking  to  the  surface  could  initiate  a  roll  cell 
2  or  3  cm  before  the  others. 

CORRELATION  OF  DATA 

Several  hundred  experiments  were  performed  wich  con¬ 
tinuous  cathodes  and  laminar  balk  flow  to  develop  a  general 
correlation  which  could  be  used  to  predict  average  mass 
transfer  coefficients. 

It  was  hoped  that  it  would  be  possible  to  define  the 
limits  of  the  flow  regime  in  which  the  superimposed  convec¬ 
tive  processes  both  play  important  roles  in  the  mass  trans¬ 
port,  so  that  all  experimental  systems  outside  this  regime 
could  be  dealt  with  by  using  the  existing  correlations  for 
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either  pure  free  or  pure  forced  convection.  Such  a  scheme 
is  appealing  because  from  physical  reasoning  it  Is  to  be 
expected  tha  ’  at  high  flowrates  and  very  small  density 
differences,  forced  convection  should  dominate  the  mass 
transfer  process.  Conversely,  if  the  bulk  flowrate  is  very 
slow  and  the  density  difference  between  the  buxk  and  inter¬ 
face  is  very  large,  natural  convection  should  dominate. 

Previous  viovk  on  iieat  transfer  from  vertically  orien¬ 
tated  plates  into  laminar  fluid  streams  (6)  indicated  that 
three  flow  regimes  could  be  defined  as  shown  in  Figure  12 . 

As  a  first  attempt  for  the  correlation  of  the  present  data, 
therefore,  the  position  of  these  two  lines  was  sought. 

First,  for  each  value  of  d/L  the  slope  and  position  of  line 
A  was  established  in  such  a  manner  as  to  result  in  the 
least  standard  deviation  of  all  data  points  lying  above  it 
from  eqn.  (3)  for  puia  free  convection  control.  Similarly, 
the  position  of  line  B  v/as  established,  yielding  the  best 
agreement  of  the  data  points  below  it  with  eqn.  (2)  for 
pure  laminar  forced  convection.  In  comparing  the  data  points 
with  the  above  correlations,  the  significant  lengths  may  be 
chosen  as  the  electrode  length,  or  the  hydraulic  diameter. 

It  was  found  that  if  the  electrode  length  was  used  in  the 
Grashof  number,  and  the  hydraulic  diameter  in  the  Reynolds 
number,  lines  A  and  B  coincide  with  a  slope  of  unity.'^ 

^  It  should  be  noted  that  the  length  dimension  in  eqn.  (3) 
cancels  out^  Ihus,  the  choice  of  the  electrode  length  as 
the  slgnlfica^.t  dimension  serves  only  to  determine  which  con- 
vectional  mode  is  dominant.  For  a  given  electrolyte  composi¬ 
tion  and  flow  velocity,  the  longer  the  electrode,  the  more 
influence  free  convectj  n  effects  vflll  have. 
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When  the  mass  transfer  coefiJcient  expressed  from 
eqns.  (2)  and  (3)  are  equated  along  the  separating  line, 
we  obtain 

\1  /3 

^  1-85  (^Re.Sc.  ^  0-19(Gr 

Cubing  both  sides  and  then  taking  logarithms,  we  have 

n 
T  ^ 

log  +  2*964  h  log  (Re. Sc)  =  log  (Gr  ^Sc) 

In  agreement  with  the  empirically  obtained  separating  lines 
log  (Gr  ,Sc)  Is  linear  In  log  (Re. Sc)  with  a  slope  of  unity 
and  the  Intercept  Is  given  by 

t2 

Intercept  -  2*964  +  log 

d"^ 

It  Is  no  longer  necessary,  therefore,  to  Include  the 
Schmidt  number  as  a  parameter,  as  It  appears  on  both  sides 

O 

to  the  same  power.  If,  in  addition,  the  term  log  of 
eqn.  (6)  Is  grouped  with  the  Grashof  number,  a  new  Grashof 
number  may  be  defined, 

gp^^p  2 

fr  a  L 

d  L  M-  PjL 

The  resulting  equation  of  the  separating  line  Is 

2*964  +  log  Re  =  log  Gr  „ 

d'^L 

This  makes  a  separate  plot  for  each  value  of  d/L  unnecces- 
sary  and  all  the  experimental  data  can  be  shown  on  one 
figure  (Fig.  13) .  The  points  Judged  by  this  scheme 
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to  be  dominated  by  free  convection  are  shown  In  Figure  14, 
as  compared  with  eqn.  (3),  while  in  Figure  15  the  points 
b'^low  the  dividing  line  (forced  convection  control)  are 
shown  togecher  with  eqn.  (s).  The  average  deviation  of  all 
the  points  showri  is  seven  percent  If  it  is  assumed  that 
the  points  in  the  two  flow  regions  should  be  exactly  pre¬ 
dicted  by  the  correlating  equations  appropriate  to  the 
respective  modes  of  control,  and  if  it  is  also  assumed  that 
all  errors  in  the  data  are  random,  then  confidence  limiis 
for  the  fit  of  the  data  to  the  equations  may  be  calculated. 
This  has  been  done,  and  the  resulting  confidence  limits 
at  the  95  per  cent  level  are  shown  on  both  Figures  14  and  15. 

The  validity  of  using  a  Grashof  number  with  (d^L)^/" 
as  the  characteristic  length  can  be  partially  rationalized 
from  the  following  considerations:  In  the  Raylelgn  number 
which  determines  the  stability  of  the  system,  the  signifi¬ 
cant  length  should  be  6,  as  given  by  eqn.  (4).  Since  this 
distance  appears  to  the  third  power  in  the  Rayleigh  num¬ 
ber,  at  the  end  of  the  transfer  section  the  Rayleigh  (or 

p 

Grashof)  number  is  proportional  to  d  L.  While  the  actual 
value  of  6  is  not  used,  the  proper  geometric  parameters  do 
appear  and  probably  to  the  correct  relative  powers.  Systems 
Judged  to  be  under  forced  convection  control  are  tho.^e  in 
which  roll  cell  formation  does  not  occur.  Those  under  free 
convection  control  have  both  types  of  flow,  but  the  average 
mass  transfer  coefficient  for  the  entire  electrode  Is 
nearly  the  same  as  for  pure  free  convection. 
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EXPERIMENTS  WITH  TURBULENT  BULK  FLOW 
1  the  case  of  turbulent  bulk  flow,  one  would  not  ex¬ 
pect  as  marked  an  influence  of  the  buoyancy  forces  as  with 
laminar  flow.  The  mass  transfer  bovmdary  layer  is  thinner 
than  for  laminar  flow,  thus  lowering  the  Rayleigh  lumber 
of  the  system.  In  addition,  the  turbulent"  eddies  continu¬ 
ally  penetrating  the  mass  transfer  boundary  layer  might 
be  expected  to  dominate  the  mass  transfer  process,  even  if 
instability  due  to  forces  existed. 

No  theory  or  empirlvcal  correJ^tlon  is  available  for 
describing  the  turbulent  mass  transfer  entry  region  in  a 
flat  duct,  making  it  impossible  to  predict  either  the 
average  or  local  mass  transfer  rates.  Again  two  of  the 
sectioned  cathodes  were  used.  Experiments  were  performed 
under  conditions  which  maximized  the  buoyancy  forces,  l.e., 
low  flowrates  and  high  copper  concentration.  Although  the 
averaged  data  shown  in  Figure  16  Indicate  that  at  Re  -  4000 
and  with  the  highest  Cu"^  concentrations  employed  there  is 
some  contribution  due  to  free  convection,  already  at 
Re  =  5000  such  differences  could  not  be  detected. 

CONCLUSIONS 

1.  In  the  presence  of  buoyancy  forces,  average  rstes 
of  ionic  mass  transport  to  horizontal  electrodes  in  chan.iel 
type  cells  may  be  predicted  with  satisfactory  accuracy  by: 
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a)  The  equation  proposed  by  Fenech  and  Tobias 

Nuj^  =  0*i9(Gr 

L 

v;hen  the  process  is  controlled  by  free  convec¬ 
tion,  or 

b)  The  equation  proposed  by  Norris  and  Streld 

Nu^  =  l*8b(Re.SG.d/X)^/^ 

U  * 

when  superimposed  laminar  flow  is  the  control¬ 
ling  convectlonal  mode. 

2.  I4unlnar  flow  oauaea  a  nonunlfonn  distribution  of 
leoal  transport  rates i  while  the  mass  transfer  ooeffloient 
deoreases  from  the  leading  edge  in  pure  laminar  flow, 
onset  of  secondary  oonveotion  (  probably  roll  eells) 
causes  an  increase  In  local  coefficients  at  a  predictable 
distance  downstream.  At  great  distances  from  the  entry 
region,  the  transport  process  is  dominated  by  free  con¬ 
vection, 

3.  Evidence  of  secondary  convection,  probably  roll 
cell  motion,  is  given  by  stratified  deposits  obtained  in 
locations  where  free  convection  contribution  to  the  trans¬ 
port  process  is  significant. 

4.  When  the  bulk  flov/  is  turbulent  (Re  >  4000),  the 
effect  of  free  convection  on  mass  transport  rates  appears 
to  be  negligible  Within  the  ranges  of  Rayleigh  and 
Schmidt  numbers,  and  d/L  ratios  considered  in  this  study. 
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NOfffiNCLATURE 


A 

a 

b 

C 

d 

n 

AE 


F 


g 

L 

n 

R 

t 

T 

X 

6 

p. 

P 

Ap 

Gr 

Gr 


d^L 


a  constant 
activity  moles/l 

subscript  denotes  tal<en  in  the  bulk 

concentration  moles/l 

equivalent  duct  diameter  cm 

2 

diffusion  coefficient  cm  /sec 
cathodic  overpotential  volts 
’"iraday’s  con^joan-*-  coul/eq 

P 

acceleration  of  gravity  cm^/sec 

2 

limiting  current  density  A/cm 
mass  transfer  coefficient  cm/sec 
electrode  length  cm 

number  of  electrons  per  ion  participating  in  the 
electrode  reaction 

universal  gas  constant  cal/niole 

transference  number  of  cupric  ion 

temperature  ®C 

location  on  electrode  from  the  leading  edge  cm 
mass  transfer  boundary  layer  thickness  cm 
viscosity  gm/cm  sec 

3 

density  gm/cm 

density  difference  between  the  bulk  and  interface 
gm/cm^ 

2  2  ? 

modified  mass  transfer  Grashof  number  p  d  LgAp/ji'p^ 

2  3  2 

rodlfied  mass  transfer  Grashof  number  p  L  gAp/p.  p^ 
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Nul  mass  transfer  Nusselt  number  k^d/D 

mass  transfer  Nu.sselt  number  kj^L/D 

Re  Reynolds  number  dVp/|i 

Sc  Schmjdt  number  M-/pD 


Section  D.,  Project  II. 


Page  88 


1. 

2. 

3. 


4. 


5. 

6 . 

7. 

8. 

9. 

10. 


11. 


12. 


13. 


1^ 

15. 

16. 


17. 

18. 


LITERATURE  CITED 

Altman,  M.,  and  F.  W.  Staub:  CEP  Symposium  Series  5^ 
(1959)  111. 

Casal,  S.  D.,  and  W.  N.  Gill:  J.  AIChE  8  (1962)  570. 

Chandra,  K.:  Roc.  Roy.  Soc.,  London  A  164  (1938)  231 

Chandrasekhar,  S.:  "Hydrod:/ naraic  and  Hydromagnetlc 
Stability”.  Or ford  Unlv.  Press,  London  1961. 

Charters,  A.  C-.:  NACA  TN  891  March  (1943). 

Eckert,  R.  G.,  and  A.  J.  Dlaguila:  Trans  ASME 
(1954)  496. 

Elsenberg,  M.,  C.  W.  Tobias  and  C.  R.  Wilke:  J. 
Electrochem.  Soc.  102  (1955)  415. 

Elsenberg,  M.,  C.  W.  Tobias  and  C.  R.  Wilke:  J. 
Electrochem.  Soc.  103  (1956)  434. 

Emmons,  H.  W. :  J.  Aero.  Sci.  (l95l)  490. 

Fenech,  E.  J.,  and  C.  W.  Tobias:  Electrochlmlca  Acta 
2  (1960)  311. 

Fenech,  E.  J.:  Ph.D.  Thesis,  Univ.  of  Calif.,  Berk., 
Dept,  of  Che  (I960). 

Gill,  W.  N.,  and  E.  D.  Casal:  J.  AIChE  8  (1962)  513. 

Gortler,  H. :  Ingenleur  Archiv.  1^  (1959)  463. 

Graham,  A.:  Philos.  Trans.  A  252  (1933)  291. 

Han,  L.  S.:  J.  App.  Mech.  82  (i960)  403. 

Hickman,  R.  G. :  The  Effect  of  Buoyancy  Forces  on 
Forced  Convection  Ionic  Mass  Transfer  at  Horizontal 
Planar  Electrodes.  Ph.D.  Thesis,  Unlv.  of  Calif., 
Berk.  1953,  Technical  Report  No.  2,  this  contract. 

Ibl,  N.,  and  R.  H.  Muller;  Z.  Electrochem.  ^  (1955) 
671. 

Ibl,  N.,  and  R.  H.  Muller:  J.  Electrochem.  Soc.  105 
(1958)  346. 

Ibl,  N. :  Advances  in  Electrochemistry  and  El'^ctro- 
chemlcal  Engineering,  Vol.  2,  Inters^’lence  F’ olishers. 
New  York  3962. 


19. 


l'&.ge  89» 


Section  D.,  Project  II. 

20.  Jackson,  T.  W.,  J,  M.  Spurlock  and  K.  R,  Purdy: 

J.  AIChE  1_  (1961)  38. 

21.  Levich,  V.  C.:  Physicochemical  Hydrodynamics, 

2nd  Ed.  Prentice-Hall,  Englewood  Cliffs,  N.J.  1962. 

22.  Lin,  C.  S.,  ^  aJ :  Ind.  Eng.  Chem.  43  (1951)  2136. 

23.  Mantzell,  E. :  Z.  Elektrochem.  ^  ^1936)  303. 

24.  Mori,  Y.:  Trans.  ASME  J.  Heat  Transfer  83  (1961)  479. 

25.  Norris,  R.  H.,  and  D.  D.  Streld;  Trans.  ASME  62 

(1940)  525.  — 

26.  Phillips,  C.  T.  :  J.  R.  Met.  Soc.  58  U932)  23. 

27.  Shaw,  P.  y.,  L.  P.  Reiss  and  T.  J.  Hanratty:  J. 

AIChE  9  (1963)  362. 

28.  Sleder,  fi.  N.,  and  G.  E.  Tate:  Ind.  Eng.  Chem.  28 

(1936)  1429.  — 

29.  Tobias,  C.  W.,  M.  Elsenberg  and  C.  R.  Wilke:  J. 
Electrochem.  Soc.  99  (1952)  359C. 

30.  Velte,  W.:  Z.  angew.  .Math,  und  Physlk  (1962)  591. 

31.  Wagner,  C.:  J.  Electrochem.  Soc.  ^  (1949)  161. 

32.  Wilke,  C.  R.,  M.  Elsenberg  and  C.  W.  Tobias:  J. 

Electrochem:.  Soc.  100  (1953)  513. 


Section  D.,  Project  II. 


Page  90 . 


Cathoda  bus  bar-| 
Contacting  post->  t 

Contacting  bushing.  \  \ 


Cothodo  hoidor  ciomp- 


Roforenco  elactrode- 
connoction 

Anodo- 


-Terminol  post 
^Signoi  rsodout 


-Colibroting  circuit 
switch 


Cothode  segnent 
'-Polyathyleno  insuiotor 
Porous  ’oftilc  channel  separator 


Figure  1.  View  of  the  experimental  cell  with  a  sectioned 
electrode  In  place. 
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Figure  2.  Schematic  diagram  of  the  experimental  ceil  and  measuring  circuit 
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PlgTire  4.  Expected  roll  cell  velocity  profile. 
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Figure  5.  The  influence  of  bulk  velocity  on  the  onset  of  roll  cel 
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Figure  6.  Transverse  deposition  profiles  for  three  systems 
as  measured  by  microscope  focusing. 
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Figure  7. 


Smoothed  data  from  the  sectioned  electrode  cell  for 
0.03  M  CuSO.  and  2.7  M  Glycero] . 
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Figure  8. 


Distance  from  leading  edge,  cm 

Smoothed  data  from  the  sectioned  electrode  cell  for 
O.cS  M  CuSO^  and  1.0  M  Glycerol. 
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Figure  9, 


Distance  from  leading  edge,  cm 

Smoothed  data  from  the  sectioned  electrode  cell  for 
0.25  M  CuSOa  and  0.0  M  Glycerol. 
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Figure  10-  Smoothed  data  from  the  sectioned  electrode  cell  for 
0.11  M  CUSO4  and  0.0  M  Glycerol. 
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Figure  11. 


Distance  from  leading  edge,  cm 

Smoothed  data  from  the  sectioned  electr*ode 
0.04  M  CuSO^^  and  0.0  M  Glycerol. 
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Figure  13.  Definition  of  flow  regimes  by  use  of  the  modified 
Qrashof  number. 
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Figure  14.  CoTiparlson  of  experimental  data  on  free  convection  con  rolled  systems  with 
the  equation  of  Fenech  and  Tobias. 
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Figure  15.  Conparison  of  experimental  data  on  laminar  forced  convection 
controlled  systems  with  the  equation  of  I'lorris  and  Streid. 
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Figure  16.  Turbulent  flow  data  from  the  sectioned  electrode  cell 
for  Re  =  4000 . 
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PROJECT  III. 

The  Role  of  Ionic  Migration  in  Electrolytic 
Mass  Transport;  Diffusivi ties  of 
[Pe(CN)g]^"  and  [Fe(CN)gl^"  in  KOH  and  NaOH  Solutions 

a)  Theoretical  Considerations 

Mass  transfer  in  an  electrolyte  can  be  described  by  four 
relationships:  a  flux  equation,  continuity  equation,  and  a 
charge  conservation  equation,  which  are  analogous  to  non- 
electrolytlc  mass  transfer  systems;  and  the  statement  of 
electroneutrality  which  is  unique  to  ionic  mass  transfer.^ 


The  principal  equations  used  are 

=  -  D^Vc^  -  U^ez^c^70  +  vc^  (l) 

7‘N^  =  -  3c^/dt  (2) 

£ 

I  P  ^  (3) 

i=i 

£ 

^  ^1 ^1  ~  ^  ^  ^  ^ 

i^ 


where  c^  is  the  concentration  of  species  i,  is  the  j.onic 
diffusivi ty,  e  is  the  electronic  charge,  P  is  the  Paraday 
constant,  I  is  the  current  density,  Is  the  ionic  flux,  t 
is  the  time,  U^,  is  the  ionic  mobility,  v  is  the  velocity  of 
the  fluid,  z^  is  the  number  of  charges  on  an  ion  and  0  is  the 
potential . 
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V/lthln  the  mass  transfer  boundary  layer  the  concentra¬ 
tions  of  all  ions  change  due  to  migration,  diffusion  and 
convective  fluxes  (equation  l).  Ions  that  migrate  towards 
the  electrode  must  diffuse  away  if  they  are  not  involved 
in  the  electrode  reaction.  The  conr entration  of  each  ionic 
species  departs  from  the  bulk  value,  but  the  solution,  to 
a  good  approximation,  remains  electrically  neutral  at  all 
points  (equation  4) . 

The  potential  gradient,  the  migration  driving  force, 
at  any  point  in  the  solution  is  obtained  by  combining 
equations  (1,2,3)  and  solving  for  V0. 

m 


1  =  1 


V0  = 

I  1=1 

K  K 

(5) 

where 

m 

K  = 

r’  2  2 

P  2^  IJ^ezp^ 

(6) 

1=1 

which  is  commonly  known  as  the  solution  conductivity. 

The  potential  gradient  Is  composed  of  an  ohmic  drop 
and  a  term  which  Includes  the  concentration  gradients  and 
the  conductivity.  The  second  term,  which  appears  in  all 
electrolyte  systems  where  a  concentration  gradient  and  dlffu- 
slvity  differences  exist,  is  the  diffusion  potential. 

To  understand  more  fully  the  effect  of  migration,  con¬ 
sider  the  cathodic  discharge  of  a  metal  ion.  Assume  that 
a  large  amount  of  supporting  electrolyte  relative  to  the 
reacting  ion  is  present,  and  therefore  the  migration  flux 
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is  negligible.  In  the  bulk  of  a  well-stirred  solution,  the 
concentration  remains  practically  constant  becau=-e  the  swiftly 
moving  fluid  rapidly  replenishes  the  loss  due  to  diffusion. 

As  the  electrode  is  approached  the  fluid  velocity  decreases 
and  the  concentration  gradients  Increase  because  the  total 
flux  is  constant  at  steady  state  (equation  2).  Concentra¬ 
tion  gradients  attain  maxima  at  the  electrode  surface  where 
the  fluid  velocity  is  zero. 

The  contribution  of  migration  can  be  increased  by 
decreasing  the  amount  of  supporting  ele  .  ’’yte.  As  in  the 
supporting  electrolyte  case,  the  fluid  velocity  decreases 
In  the  vicinity  of  the  electrode.  However,  a  migration 
velocity  (towards  the  electrode  for  the  metal  deposition 
case  being  considered)  affects  the  entire  concentration  pro¬ 
file.  The  reacting  ion  in  the  faster  moving  lamina  migrates 
to  the  slower  lamina  nearer  the  electrode,  therefore  shifting 
the  concentration  profile  towards  the  electrode  surface  as 
shown  in  Figure  1. 

The  change  in  the  concentration  profile  of  the  reacting 
ion  is  of  Interest  because  the  limiting  current  1^^  is 
directly  related  to  the  concentration  gradient  at  the  inter¬ 
face.  Combining  equation  (l)  and  (5),  one  obtains 


II  =  nF 


(7) 
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fuantltatlve  Calculation  of  the  Limiting  Current  for  a 
pedal  Case 

Until  this  point  a  consideration  of  the  physical  ’Tiodel 
has  not  been  necessary,  but  a  quantitative  analysis  can  be 
performed  only  on  a  specific  flow  systens.  The  rotating 
disk  is  one  of  the  few  physical  models  where  the  velocity 
distribution  is  known  within  the  vicinity  of  the  surface.^ 
The  limiting  current,  including  the  migration  contri¬ 
butions,  to  a  rotating  disk  electrode  can  be  calculated  for 
a  constant  potential  gradient  in  the  mass  transfer  boundary 
layer.  A  constant  potential  gradient  is  possible  if  all 
ionic  dlffuslvities  are  equal  and  the  solution  contains 
enough  supporting  electrolyte  so  that  the  conductivity 
remains  essentially  constant  through  the  boundary  layer. 
Then  the  potential  gradient  from  equation  (5)  will  be  (at 
limiting  current) 

7$  =  -  (8) 

Applying  a  ccr.sjoant  potential  and  eliminating  the 
flux  between  equation  (2)  and  equation  (l),  one  obtains 
at  steady  state 

=  (v  -  (9) 

where  the  velocity  gradient  was  eliminated  by  the  conti¬ 
nuity  condition  and  the  fluid  properties  were  assumed  to 
be  constant. 

For  the  rotating  disk  electrode,  equation  (9)  can 
be  simplified  further.  Levlch  (  )  has  proposed  that 
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the  concentration  is  a  function  of  only  one  dimension,  y, 

which  is  the  perpendicular  distance  from  the  plane  of  the 
disk 


d^c. 


D- 


dy‘ 


p 


=  (y  -  U^ez^  4^)  ^ 
1  1  dy^  dy 


(10) 


In  equation  (10)  v  and  d0/dy  are  known  as  functions  of  y. 


These  functions  do  not  Involve  c-  or  its  derivatives.  There¬ 
fore,  we  have  a  linear  second-order  differential  equation 
which  can  be  readily  solved  f<  and  related  to  the 
limiting  current,  I^.  The  solution  to  equation  (8)  which 
satisfied  the  boundary  conditions: 


dc- 
•L'  W 


is 


(11) 


The  fluid  velocity  component,  v,  in  equation  (9) 

is^ 

v=  (v'.n)^/2(-0.510(^)e2+0.333(^)^/2r-0.1C3(^)^e^+....  (la) 

where  oj  Is  the  angular  frequency  of  the  disk  and  v  is  the 
kinematic  viscosity  of  the  solution. 

The  development  is  simplified  by  letting 


X  =  • 

According  to  Levlch;  6^  --  bj|(Sc)  (diffusion  boundary 
layer),  and  6^^  =  1.805  (— )  ^  (hydrodynamic  boundary  layer), 
where  Sc  =  ^  (Schmidt  number). 
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Perfoiming  tlie  variable  transformation  to  equations  (ll) 
and  (12)  we  have 


=  J  “g-x®+.B85(So)'VV-.103 


U^ez^V06. 


dx  (13) 


The  value  of  the  limiting  current,  from  equation  (6)  and 
equation  (.-'3),  is: 


I.  =  - 

Ll 


nF  D^c° 


(14) 


where 


exp(^-x^+.885(Sc)"^/^x^-.103(Sc}‘^/^x^. . 


/  U,ez,V06,\ 

exp( - ^ - ]  dx  (13) 

\  •^l  / 

The  migration  integral,  M, ,  is  a  function  of  the  para- 
meter  - J  which  can  be  related  to  solution  proper¬ 

ties  by  the  application  of  equations  (7) , (ll) 

^1®^1^^  ^1®"1^L^1  ^  ^1 

°1^1  ^1^1 

where  t^  is  the  transference  number  of  the  reacting  ion  in 
the  bulk. 

t,  ^ 

^  U  ez^oj 

Equation  (15)  can  be  expanded  into  an  infinite  series  of 

gamma  functions.  When  evaluated,  neglecting  cross  products 

♦ 

and  higher  terras,  one  obtains: 

#  - 

Note:  Se^yting  t]^^0  in  equation  (16)  gives  the  relationship 
Mi=.8934-f .2747  Sc"^/^,  which  compares  with  the  empirical  and 
commonly  used  Gregory  and  Rlddiford  correction,^ 
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1  /r,  n  t, 

=  .893  +  ,2  747  -  .4913  (“^y 


n  t-,  /iit- 

'l‘‘l 

Figure  2  gives  the  values  of  ,  for  large  Schmlr't 
numbers,  as  a  function  of  nt^/s^. 


Comparison  of  Migration  Correction  Methods 

A  comparison  can  be  made  between  the  commonly  used 

migration  correction  and  the  present  results.  By  equating 

^l^L 

the  migration  flux  to  the  product,  — — 5-,  we  have 


2^ 


nt,^ 


II  (1-  - 


(17) 


for  an  arbitrary  mass  transport  system.  Applying  the 
effective  film  thickness,  6^,  for  the  rotating  disk  elec- 

•r 

trode  equation  (17)  becomes 


- 


b 

nF  D^c^ 


.89341 


(18) 


(since  7c®  =  ®l/®l  ~  c^/.8934  5^)  whereas  from  equation  (14) 


and  equation  (16),  at  small  t^  and  large  Sc^  we  have 

b 


=  - 


nF  D^Cj 


.8934^ 1  -  .502 
\  z 


^1  ^ 


(19) 


The  transference  number  appears  In  a  rigorous  solutions 
as,  approximately,  the  arithmetic  average  over  the  boundary 
layer  (t^^O  at  the  electrode  surface).  In  equation  (18) 
the  t^  was  considered  to  be  constant. 

To  clarify  further  the  role  of  Ionic  migration  the 
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effect  of  a  changing  potential  gradient  within  the  boundary 
layer  should  be  estimated. 


Changing  Potential  Gradient  Effect  on  Mass  Transfer 

The  Interfaclal  potential  gradient  can  be  estimated 
for  an  arbitrary  number  of  Ions  by  allowing  the  Interfacial 
concentrations  to  be  represented  by 

(20) 


O  b  „  Oe» 
Cl  =  -  Vc^5^ 


where  5^  Is  the  equivalent  dll '‘uslon  In  film  thickness. 

To  satisfy  electroneutrallty  all  6^  must  be  equal  when 
a  linear  representation  such  as  equation  (20)  Is  used  for 


0 

'1* 


c; .  At  I  =  c°  =  0,  hence 


whe  re 


and 


^1  ^  ^'^'l  "  ^1' 


b  ,  b 

7l  »  c^/o^ 


P,  =  7o°/Vo° 


(21) 


By  eliminating  the  potential  gradient  between  the  flux 
equations  for  the  non -discharging  Ions,  applying  the  Nernst- 
Elnsteln  relationship  (Dj^--kTUj^) ,  eliminating  the  Inter¬ 
faclal  concentrations  by  equation  (2l)  and  solving  for 
we  obtain  for: 

Case  1.  Ion  discharge  (deposition)  reaction 
6  = 

1  ^2^2  ^2  ^  ^1  ~  ^2  ^ 


arid  since 


m 


m 
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by  electroneutrality  I3  given  by 


m 


^2 '''2 


0 


and  ana-LOgously  for 

+z,  +s,. 

Caae  2*  Oxidation-Reduction  (Ox  )e  =■  Red.  ^ ) 


1  2  ^3  + 


and 


+ 


m-l 

I 


\  j.  j.  a.  114  —  A 

i«*3 


^2  = 


-Di/Dg  22(73-^3)  +  ^2'Ye^3 

Vi's  P3'"2-"3^ 


^1  “  ZjYj  +  ’ 


The  Interfaclal  potential  gradient  can  be  written  as 

„^o  "^L  [l+I*] 


when 


- - - 

2 

eZiTi 


from  equation  (6)  and  equation  (21) 
and  tn  r. 


1-1  ^ 


A  rigorous  calculation  of  a  proper  average  potential 
gradient  across  the  mass  transfer  boundary  layer  could  be 
a  complex,  and  perhaps,  an  Impossible  task  by  analytical 
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techniques.  For  purposes  of  Illustration  an  arithmetic 
average  gradient  has  been  selected.  Physically  an  arithmetic 
average  would  not  provide  a  proper  weighting  of  the  more 
influential  interfaclal  potential  gradient  where  the  fluid 
velocity  approaches  zero  and  migration  and  diffusion  pre¬ 
vail.  However,  this  analysis  is  only  intended  to  provide 
additional  insight  into  the  effect  of  a  change  in  poten¬ 
tial  gradient  within  the  boundary  layer.  JThls  average 
gradient  is 


The  results  for  remain  equivalent  for  a  constant 
average  potential  gradient  if  an  effective  transference 
number  is  defined.  Applying  the  potential  gradient  of 
equation  (20),  we  arrive  at  the  following  equation  for  the 
effective  transference  number 


t’ _ ± 

^1  “  2 


1  + 


l+l' 

1+a 


■■ 


(23) 


which  embodies  corrections  for  a  diffusion  potential  and 
conductivity  changes  in  the  boundary  layer  that  appear  in 
real  systems. 


Numerical  Illustrations  for  an  Average  Potential  Gradient 
Calculations  were  performed  for  several  typical  elec¬ 
trolytes.  The  parameters  I*,  a  and  t^  were  calculated 
from  limiting  ionic  conductances  at  25,0®C*^  Since 
only  values  of  ionic  properties  are  involved,  the  use  of 
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limiting  conductHnces  is  not  too  unrealistic  and  suitable 
for  illustrative  purposes.  The  results  of  these  calcu¬ 
lations  are  shown  in  Table  I. 

In  each  of  these  tables  the  effective  transference 
number  t^  and  the  transference  number  in  the  bulk,  t^, 
have  been  compared.  If  the  supporting  electrolyte  is 
a  better  conductor  than  the  reacting  ion,  the  diffusion 
potential  is  positive,  thereby  increasing  the  effective 
transference  number.  At  the  other  extreme,  as  shown  in 
Table  I,  the  diffusion  potential  is  negative,  hence 

^1  ^  ^1* 

Conclusior.3 

1.  The  changing  concentrations  in  the  mass  transfer 
boundary  layer  are  important  in  determining  the  contribu¬ 
tion  of  ionic  migration  to  the  limiting  current.  For  the 
special  case  of  a  constant  potential  gradient,  the  limiting 
current,  including  migration  effects,  can  be  calculated 
rigorously. 

2.  The  diffusion  potential  contributes  to  the  elec¬ 
tric  field  which  promotes  migration.  Although  an  approxi¬ 
mate  calculation  was  used  to  Illustrate  the  effect  of  the 
diffusion  potential,  additional  insight  about  the  capacity 
of  various  supporting  electrolytes  for  suppressing  ionic 
migration  vras  obtained  by  evaluating  effective  transference 
numbers  for  some  typical  binary  and  quaternary  electrolyte 
systems.  Highly  conducting  supporting  electrolytes 
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significantly  increase  the  effective  transference  number  of 
the  reacting  ionic  species  decreasing  the  effectiveness 
of  these  salts  for  suppressing  ionic  migrations. 


Dc..ditle3  and  tlscositles  of  Aqueous  Solutions  of 
K3|Fe(CN)g|  -  K^|Fe(CN)gl  -  KOH  (or  NaOK),  and  the 
Diffusion  Coefficients  of  the” b’errous-  and  Ferricyani de 


Ions . 


Density,  viscosity,  and  dlffuslvlty  data  were  obtained 


in  equimolar  ferrous-  ferricyanide  solutions  (0.01  -  0.20 
molal)  in  the  presence  of  NaOH  or  KOH  (0.0  -  3.5  molal). 


Solution  Preparation  and  Analysis 


Solutions  were  prepared  with  distilled  water  and  re¬ 


agent  grade  chemicals.  Ferricyanide  ion  was  titrated  with 

g 

thiosulfate  and  ferrocyanlde  was  determined  by  a  ceric 

7 

titration  .  Hydroxyl  ion  was  titrated  with  standard  HCl 


to  a  bromthymol  blue  endpoint. 


Densities  and  Viscosities 

A.  Experiment. 

Densities  and  viscosities  were  determined  by  stan¬ 
dard  techniques  (pycnometer  and  U-tube  viscometer^  s 
respectively).  All  solutions  were  thermostated  to  25.00±.05®C 
before  measurements  were  taken. 

B.  Results. 

A  second  degree  polynomial  involving 
(average  of  ferri-  and  ferrocyanlde  concentrations)  and 
^KOH  ^NaOH  Pitted  to  the  data  by  a  least  squares 
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technique.  A  total  of  26  points  for  KOH  supported  solutions 
and  28  points  for  NaOH  solutions  were  used  in  the  curve 


plotting. 

The  equations  and  their  associated  standard  errors  of 
estimate  are: 

a.  Potassium  hydroxide  as  a  supporting  electrolyte. 

H[cpl  =  .8937  +  .6178  +  .1064 

+  .00656  +  .0205  (1) 

Standard  Error  of  Estimate  =  l.OjS. 
ptg/cp]  =  .9971  +  .3989  +  .0502 

-  .000140  C™.  -  .0032  0„.„  C„  . 

KOh  KOH  Redox  v*:; 

Standard  Error  of  Estimate  =  0.5^. 

b.  Sodium  hydroxide  as  a  supporting  electrolyte. 

p[cp]  =  .8937  .H  .6969  +  .1539 

+  •'^235  4  .0450  +  .2527  <3) 

Standard  Error  of  Estimate  =  1.0^. 
ptg/op]  =  .9971  +  .3996  +  .0427 

.00113  -  .00657  CNaCH'^Redox  (4) 

Standard  .Error  of  Estimate  =  0.1^. 

Ranges  of  validity:  Cj^g^lox 

^KOH’  ‘^NaOH  • 

The  data  compared  with  literature  values,  through 

equations  (l)  -  (4),  as  siiown  in  Tables  II-IV. 
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Dlffuslvltles 


Experiment 

A  rotating  disk  electrode  was  applied  for  the 

measurement  of  the  Ionic  diffusion  coefficients  of  the 

ferrous  and  ferrlcyanides‘1'  '  .  This  tool,  first 

analyzed  by  Levlch  (S)  for  a  binary  electrolyte,  then 

13 

treated  by  Hogge  and  Kralchman  for  a  completely  flooded 
electrolyte,  h.as  been  extended  to  partially  flooded  systems 
where  migration  effects  can  not  be  neglected. 

The  theory  developed  In  part  a  was  used  to  determine 
the  dlffusivlty  of  the  reacting  Ion  from  limiting  current 
measurements  by  a  rotating  disk  electrode.  The  equation 


used  is: 


r-1.805  I^v  f  A 

’i °  — J 


nF  cjoi-v*-  J 

Migration  effects  are  included  in  the  parameter  M^ 
which  contains  the  effective  transference  number  (see  part  a) 
of  the  reav^tlng  Ion.  has  a  value  of  .893  +  ,2661  Sc^”^/'^ 
when  the  transference  number  of  the  reacting  ion  Is  zero. 

The  electrode  was  a  1  cm  Ft.  disk  embedded  In  epoxy 
resin  and  attached  to  an  aluminum  spindle  rotated  at  pre¬ 
cisely  controlled  speeds  between  600  -  3600  rpm.  Figure  1 
Is  a  schematic  drawing  of  the  apparatus. 

Results 

From  the  limiting  current  measurements  the  diffu- 
slvltles  of  ferrl-  and  ferrocyanlde  ions  were  determined 
by  equation  (5).  The  data  were  correlated  in  the  form  of 
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the  Stokes-Einstein  relation,  Du/T,  as  a  function  of  the 
Ionic  strength  F. 

m 

r  =  1/2 

Preliminary  calculations  performed  on  the  data  of  Nelson  and 
Riddlford  Indicated  that  D,u/T  was  independent  of  M,  and  T 
at  r  =  const. 

Although  temperature  was  not  a  variable  In  the  present 
study,  the  parameter  D|.l/T  should  be  Independent  of  tempera¬ 
ture,  since  the  behavior  of  the  ferrois-  or  ferric  cyanide 
Ion  should  be  similar  In  this  respect  to  that  of  the  large 
trl -Iodide  ion. 

A  least  squares  fit  of  the  data  (values  obtained  v/ith 
KOH  or  NaOH  agreed  v/lthln  range  of  experimental  error, 

±45^)  gives  the  following  empirical  relationships  and  stan¬ 
dard  errors  of  estimate: 

^(fer’rl)  -  (.234  +  .0014  i)  10'^  —  (7) 

±  .005 

^(ferro)  =  (.187  +  .0036  F)  x  10”^  (8) 

1  se  c 

±  .007 

15 

As  a  comparison,  the  data  of  Elsenberg,  ejt  al  gives 

^(ferrl)  =  .250  x  10"^  — 

T  '  '  see  K 


and 


^(ferro)  =  .215 


10-3 


cm 

sec 


poise 


whereas  equations  (7)  and  (8)  give 
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^(ferrl)  =  .240  x  £$13®. 

f  (ferro)  =  .202  x  lO’^  ii-  Eg|2® 

Although  no  conclusion  can  be  reached  about  the  dis' 
crepancies  between  the  data  presented,  the  inherent  errors 
in  the  effusion  capillary  technique  used  by  Eisenbex-g  e^  ^ 
could  be  responsible  for  the  differences.  These  empirical 
relationships  are  considered  adequate  representations  of 
the  data  for  use  in  mass  transfer  studies,  and  are  rcllab'ie 
within  the  error  estimates  presented. 
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NOMENCIATURE 
Lower  Case  Letters 

c  concentration  -  moles/ctn^ 

e  electronic  charge  -  coulombs 

g  parameter  -  dlmenslcnle’us 

j  Index  of  series  *  dimensionless 

k  Toltzr  an  constant  ••  erg/dcg 

m  otal  number  of  spe'les  -  dimensionless 

n  number  of  electrons  reacting 

t  time  sec 

t  transference  number  -  dimensionless 

V  fluid  velocity  -  cm/sec 

V  migration  velocity  -  cm/sec 

X  transformed  variable  -  dimensionless 

y  coordinate  axis  -  cm 
z  charge  on  an  Ion  -  dimensionless 

Upper  Case  Letters 

A  constant  of  integration  -  moles /cm^-cm 
B  constant  of  Integration  -  raoles/cm^ 

C  concentration  -  moles/liter 
D  diffusion  coefficient  -  cm^/sec 
F  Faraday  constant  -  coulombs/equlv 

p 

I  current  density  -  amp /cm 
I  "diffusion  current"  -  dimensionless 
K  solution  conductivity  -  ohm''^(’m~^ 

M  molarity  -  moles/llter 
M  migration  integral  -  dlme.islonler  s 
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N  Ionic  flux  -  moles/cra  -sec 
R  gas  constar  -  -  ergs/deg-mole 
T  absolute  temperature  ~  deg.  absolute 
U  mobility  -  cm-sec~^/dynes 

Greek  letters 

a  fractional  conductivity  change  -  dimensionless 
P  relative  concentration  gradient  -  dimensionless 

y  relative  concentration  -  dimensionless 

6pj  hydrodynamic  boundary  layer  -  cm 
5  diffusion  boundary  layer  -  cm 
6'  equivalent  diffusion  thickness  -  cm 
€  dielectric  constant  -  dimensionless 
X  Ionic  conductance  -  ohm  cm  /equlv 
[I  bulk  vr.scosity  -  poise 

o 

V  kinematic  viscosity  -  cm  /sec 
i  variable  of  Integration  -  cm 

variable  of  Integration  -  cm 

3 

p  density  -  gm/cm 
<t>  electric  potential  -  volt 
0)  angular  frequency  -  radians /sec 
r  Ionic  strength  -  moles/liter 

Subscripts  and  Superscripts 

b  bulk 

1  Index  of  species 

o  Interfaclal  property 
L  limiting  current 
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Table  1.  (cont.)  Effect  of  Supporting  Ionic  Species  on  the  Effective  Transference 

Numbers  of  Reacting  Ions  for  Various  Concentration  Levels 
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Tab le  I . 

Coipparlson  of  Viscosity  and  Density  Data  with 
Those  of  Elsenberg  al^^ 

Conditions:  1.  -  =  25.0^C. 


K^Fe (CN)p 

K^Fe(CN)K 

NaOH 

d 

Dev. 

P 

P 

Dev. 

(molar) 

(molar) 

(molar) 

(lit.) 

( calc . ) 

{^} 

(lit.) 

( calc) 

.00998 

.01011 

1.958 

1.368 

.1379 

+0.8 

1.0807 

1.0807 

+0.0 

.01224 

.01228 

2  .011 

1.382 

1.400 

+1.3 

1.0833' 

1.0834 

+0.0 

.04^96 

. 05052 

1 . 92  3 

1.394 

1.415 

+1.5 

1.0949 

1.0952 

+0 . 0 

.0994 

.1008 

2  .007 

1.494 

1.505 

+0.7 

1,1174 

1.1169 

-0.0 

.1994 

.2030 

2  .004 

1.642 

1.626 

-1.0 

1.153C 

1.1560 

+0.2 

Table  II. 

16 

Comparison  of  Viscosity  Data  with  Those  of  Lange 

Conditions:  1.  Pure  sodium  or  potassium  hydroxide 
2.  ^  =  25.0®C. 


KOH 

(molar) 

(lit,) 

1 

(calc . ) 

Dev. 

if) 

NaOH 

(molar) 

M  1 
(lit.) 

(calc . ) 

Dev. 

1.00 

1.009 

1.007 

-0.2 

1.00 

1.104 

1.093 

-1.0 

0.50 

0.951 

0.949 

-0.2 

0.50 

0.991 

0.982 

-0.9 

0.25 

0.922 

0.921 

-0,1 

0.25 

0.944 

0.935 

-0.8 

0.125 

0.905 

0.907 

-t-0.2 

0.125  1 

0.921 

0.913 

-0.7 

Table  III. 

Comparison  of  Density  Data  with  Those  of 
Hitchcock  ^  al^ ' 

Conditions:  1.  Pure  sodium  and  potassium  hydroxide 
2.  T  =  25.0®C. 


KOH 

_(molar) 

,  p  , 

(lit.) 

( calc . ) 

Dev. 

(?) 

NaOH 
(molar ) 

,  p  , 
(lit.) 

P 

(calc . ) 

Dev. 

(?) 

4.00 

1,1726 

1.1754 

+0.2 

4.00 

1.1491 

1 . 1500 

+0.1 

3.00 

1,1314 

1.1351 

+0,3 

3.00 

1.1145 

1.1154 

+0.1 

2  .00 

1.0886 

1.0919 

+0 » 3 

2  .00 

1.0776 

1.0803 

+0.3 
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Pig.  1  Alteration  of  Concentration  Profiles 


by  Ionic  Migration 
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Figure  3.  Cross-sectional  view  of  the  rotating  spindle. 
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Figure  4.  Schematic  dlagrauri  of  the  rotating  cell  apparatus 
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PROJECT  IV. 

Mass  Transfer  in  Concentrated  Binary  Electrolytes 
Introduction 

The  flux  of  minor  components  in  a  dilute  electrolytic 
solution  may  be  expressed  as^ 

=  -  z^u^F  C4  V  <&  -  V  +  c<v_  .  (1) 

This  equation  clearly  accounts  for  three  mechanisms  of 
mass  transfer:  the  motion  of  a  charged  species  in  an  elec¬ 
tric  field  (so-calj.ed  migration),  molecular  diffusion  due  to 
a  concentration  gradient,  and  convection  due  to  the  bulk 
motion  of  the  medium.  Although  this  equation  has  been  applied 
fruitfully  to  many  electrochemical  problems,  its  validity  is 
restricted  to  dilute  solutions. 

The  problem  is  most  conveniently  considered  in  two 
parts.  First  one  seeks  transport  equations,  such  as  equation 
(1),  which  correctly  relate  driving  forces  and  fluxes  and 
thereby  Introduce  certain  transport  properties,  such  as 
and  u^ .  The  transport  properties  should  be  measurable  state 
properties  v^hlch  should  be  independent  of  the  driving  forces 
such  as  V<I>  and  Vc^^. 

The  second  part  of  the  problem  is  to  characterize  the 
variation  of  the  transport  properties  with  temperature, 
pressure,  and  composition.  This  usually  takes  either  the 
form  of  a  theoretical  linkage  to  more  fundamental  concepts 
of  the  structure  of  matter  or  the  form  of  an  empirical  corre¬ 
lation  of  available  measurements. 

Transport  Equations 

Equation  (l)  is  Inadequate  even  in  ternary  solutions 
of  non-electrolytes,  since  in  these  there  are  two 


■’Icv'  oJ  on  U.  i  ]^-  'Jn  -  L  I\' . 
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indopcnidont  concentration  <.;radlent3  and  tlu;  diri'na.i  ^n  flux 
of  each  species  can  be  affected  by  both  concentratl  m  gradients 

As  pointed  out  by  Guggenheim  *  »  ^  the  electrostatic 

potential  In  equation  (1)  does  not  have  a  clearly  defined 
physical  significance.  The  mobility- u,^  and  the  diffusion 
coefficient  are  no  longer  related  by  the  Nernst-Elnsteln 
equation  at  high  concentrations.  Few  Independent  measure¬ 
ments  of  single  Ionic  mobilities  and  diffusion  coefficients 
are  available. 

In  order  to  avoid  these  difficulties,  equation  (1)  can 


be  replaced  by 


•i  ^  '"1  = 


This  equation  is  similar  to  the  Stefan -Maxwell  equation'’ 
and  Is  equivalent  to  one  developed  by  Onsager  .  The 
Stefan-Maxv/ell  equations  apply  to  diffusion  In  dilute  gas 
mixtures  and  express  the  driving  force  as  a  mole  fraction 
gradient  or  a  gradient  of  partial  pressure  Instead  of  the 
gradient  of  the  electrochemical  potential.  The  reciprocal 
of  the  can  be  regarded  as  friction  coefficients  simi¬ 

lar  to  those  used  by  Lalty®^^  and  Kiemm"^*^'^^  to  describe 

■JO 

transport  in  ionic  solutions  and  melts.  Burgers  has  also 

used  this  concept  to  treat  the  conductivity  of  ionized 

gases,  and  Lightfoot,  et  ^  ,  have  applied  aquation  (2) 

14  1 S 

to  liquid  solutions.  See  also  Lamm  and  Truesdell'  . 


Equation  (14;,  page  245,  in  reference  7. 
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The  term  -c^vp,^  in  equation  (2)  can  be  regarded  as 
a  driving  force  per  unit  volume  acting  on  species  i  and 
causing  it  to  move  with  respect  to  the  surrounding  fluid. 

The  force  per  unit  volume  exerted  by  species  j  on  specieo  1 
as  a  result  of  their  relative  motion  Is  expressed  as 


that  is,  proportionp.l  to  the  difference  in  velocity  of 
the  two  species.  By  Newton's  third  law  of  motion 
Equation  (2)  thus  expresses  the  balance  between  the 
driving  force  and  the  total  drag  force  exerted  by  the  other 
species , 

The  number  of  Independent  equations  of  the  form  (2) 
is  one  less  than  the  number  of  species.  When  equation  (2) 
is  summed  over  1,  the  3 eft  side  is  zero  by  the  Glbbo-Duhem 
relation,  and  the  right  side  is  zero  by  Newton's  third  law 
of  motion. 

In  contrast  to  equation  (l),  equation  (2)  is  not  a 
direct  expression  for  the  fluxes.  However,  the  material 
balance  for  a  species  reads 

dc^ 

(3) 

and  requires  an  expression  for  the  flux.  In  particular 
cases  we  invert  equation  (2)  to  obtain  these  expressions. 

For  example,  for  *  binary  electrolyte  composed  of  anions, 
cations,  and  solvent,  equation  (2)  yieDds  two  Independent 
equations  which  may  be  rearranged"^^  to  read 
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=  c  V  =  - 


V  ^  it* 


V  •&- 


i  t* 


vRT  °  ^  F  z  ■*■ 


^  "  ■  'R^  °o  ^o  ■’■  P  (  z  z  )  "^ 

X  ^  •» 


The  molar  average  velocity 


,*.± 

Cm 


has  been  chosen  as  the  reference  velocity,  and  the  current 
density  is 

1  =  F  (z^c^X^.  +  z.c_:3L^). 

The  parameters  tj,  and  t*  in  equations  (4)  are 
transport  properties  related  to  the  of  equation  (2). 
The  transference  numbers  (with  respect  to  the  molar 
average  velocity)  are 

c  c^  z.^. 

t*  =  1  -  t*  =  —  +  — -  -F-  — -7C.-  .  i 

The  diffusion  coefficient  of  the  salt,  based  on  a  thermo¬ 
dynamic  driving  force,  is 


^  ~  - = — ie: - r — jss — 


Z^«C^^-Z  iCr 
+  0+  -  O- 


The  diffusion  driving  force  used  in  equations  (4)  is  the 
gradient  of  the  chemical  potential  of  the  electrolyte 
in  the  solution.  This  chemical  potential  is  readily 
measurable,  and  no  reference  to  single  ionic  activity 
coefficients  is  necessary. 
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Now  we  want  to  define  a  "potential  the  solution” 
for  use  as  a  driving  force  for  the  current.  For  this 
purpose  the  potential  of  a  suitable  reference  electrode, 
at  the  point  In  question,  with  respect  to  a  fixed  reference 
electrode  Is  used.  Application  of  thermodynamic  principles 
to  such  a  reference  electrode  yields  In  this  case 

s_  V  M._  +  V  +  Sq  V  Pq  =  -  nF  V  <!>,  (?) 

where  the  electrode  reaction  Is  given  by 

z  z 


s  M  +  s.M.  + 

-  -  +  +  o  o 

16 


ne 


Equation  (7)  can  be  rearranged  to  give 


1  =  -  KVt  - 


F 


!±!± 

n 


z  ^ 

■f  o+ 

z  -z  /)r  ■ 
-t  o+  -  o- 


s  z  c 
o  +  + 


n  c 


o  J 


where  the  conductivity  k  Is  another  transport  property 
and  Is  related  to  the  of*  equation  (2)  by 

z^z_F^c^  (c_^  + 

K  -  c  +  c.  +  c  /Qr 

-  o+  +0-  o  +- 

Note  that  the  potential  defined  here  as  the  potential  of 
a  reference  electrode  is  considerably  different  from 
the  electrostatic  potential  used  in  equation  (l).  This 
definition  avoids  the  ruestlonable  concepts  regarding 
potentials  in  the  solution  brought  up  by  Guggenheim. 

The  transport  parameters  usually  considered  In  a 
binary  electrolyte  are  the  conductivity,  the  diffusion 
coefficient,  and  the  transference  number,  and  these 
are  related  to  ,  and  by  equations  (5),  (6), 

and  (9).  Actually,  the  transference  number  with  respect 


(8) 


(9) 
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to  the  solvent  velocity  I3  usually  measured  and  reported. 
This  is 


1  -  t°  = 


z  S' 

-t  Cf 

z  -  z”  'S' 
+0+  -  o- 


The  diffusion  coefficient  of  the  salt  which  is  usually 
measured  and  reported  is  related  to  the  thermodynamic 
diffusion  coefficient  by 


D  1  + 


d  iH  ■ 
c5  iH  in 


y-)  . 

mj 


Tne  difference  between  D  and  ^  is  due  to  the  choice  of 
the  driving  force  and  to  the  choice  of  the  reference 
velocity. 


(10) 


(11) 


Concentration  Dependence  of  Transport  Properties 

The  theory  up  to  this  point  is  macroscopic  In  nature. 

In  order  to  gain  further  insight  into  the  nature  of  the 
it  is  necessary  to  resort  to  microscopic  theories 
or  empirical  Interpretation  based  on  experimental  obser¬ 
vations.  From  the  latter  it  is  known  that  over  moderate 
ranges  of  concentration  the  product  of  the  diffusion  coeffi- 

g 

cient  and  the  viscosity  n  is  raasonably  constant  .  On 
this  basis,  and  |i*^  may  be  assumed  constant  for 

ion-3olvent  Interactions,  where  electrostatic  forces  are 
small.  Such  an  assumption  will  not  account  for  the  inter¬ 
actions  between  ions.  The  Debye -Huckel-Onsager  theory  of 
conductance  and  the  experimental  confirmation  of  this 
theory  show  that  is  proportional  to  •fo  in  dilute 
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solutions.  As  a  first  approximation  we  shall  assume  that 
this  dependence  persists  at  higher  concentrations.  Thus 
we  write 


=•  ICT/A  , 


=  kT/B  , 


)  (12) 


where  A,  B,  and  0  are  assumed  to  be  constants,  independent 
of  concentration. 

Substitution  into  equation  (9)  yields 
•A»CCm 


-  z^z_eF 


KH 


=  G  + 


AB 

v+A  ^  v.t 


(13) 


For  A,  B,  and  Q  to  be  constant  thus  requires  tha>^  a  plot 
of  •J vcc,j,/k|i  against  should  be  a  straight  line. 

Such  plots  are  shown  in  figures  1  and  2  for  aqueous  KOH 
solutions  and  in  figures  3  and  4  for  aqueous  KCl  solutions, 
both  at  25®C,*  The  value  of  G  can  be  determined  from  the 
Intercept,  and  the  slope  gives  the  value  of  AB/(v^A+v_B) . 

Substitution  of  equations  (12)  into  equation  (lO)  gives 


«  ^  V  B 

1-0  _  1  4.^  _  ** 

'+  -  1  -  +  v'JS  > 


(14) 


which  suggests  that  the  transference  numbers  with  respect 

to  the  solvent  velocity  should  be  constant.  For  KOH, 

o  ?  9 

t^  =  0.23  is  approximately  constant  between  1  and  17  molal 

concentration.  The  transference  number  for  KCl  is  near 


# 

The  values  of  viscosity  \rere  taken  from  references  17,  18, 
and  19  for  KCl  and  from  references  20,  21,  and  22  for  KOH. 

The  values  of  conductivity  were  taken  from  references  23,  2-*, 
25,  and  26  for  KCl  and  from  references  20,  21,  27,  and  28  for 
KOH. 
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^  0.489,  although  a  value  as  low  as  0.4857  at  3  molar 
has  been  measured  by  the  Hlttorf  method^^*^^. 

The  assumption  that  A,  B,  and  G  are  constant,  together 
with  the  slopes  and  Intercepts  from  figures  l”-4  and  the 
above  values  of  the  transference  numbers,  enables  us  to 
calculate  the  following  values  of  A,  B,  and  G; 

I  M  _  T  n  n  /\8  n  _  n  rtB 


A  X  10® 

B  X  10® 

G  X  10 

cm 

cm 

cm 

24.06 

23.02 

39.71 

26.89 

8.63 

9.63 

This  determination  of  A,  B,  and  Q  requires  values 
for  the  viscosity,  conductivity,  transference  number,  and 
density.  Then  It  Is  possible  to  predict  the  diffusion 
coefficient.  Equation  (6)  becomes 

.  I 

^  v^A  +  v_b 


(15) 


Reliable  data  on  the  diffusion  coefficients  of  KOH  are 
not  available.  However,  for  KCl  precise  diffusion  coeffi¬ 
cient  data  are  avallabie^^ Th-  s  values  of 
and  can  be  calculated  dlrectl;;  from  experimental  data 
without  the  assumption  that  A,  B,  and  0  are  Independent 
of  concentration.  The  results  are  plotted  In  figure  5. 

The  corresponding  values  of  A,  B,  and  G  as  functions  of 
concentration  are  shown  In  table  I. 
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Table  I. 


Transport  Parameters  for  KCl 


c 

A  Tc  10'^ 

B  X  10® 

G  X  10' 

mole 

liter 

crn 

cm 

cm 

0.0001 

23.53 

22  .65 

58.84 

0.0005 

23.53 

22.63 

72.72 

0.001 

23.53 

22.63 

76.54 

0.005 

23.48 

22.59 

79.17 

C.Ol 

23.45 

22  .56 

78.60 

0.05 

23.25 

22.34 

80.08 

0.1 

23.11 

22.18 

70.69 

0.5 

23.03 

22.01 

53.87 

1.0 

23.53 

22.48 

46.29 

2  .0 

24.41 

23.28 

38.83 

3.9 

25.21 

23.99 

34.47 

An  inspection  of  table  I  shows  that  A  and  B  do  not 
vary  from  the  predicted  values  by  more  than  5^.  This  is 
satisfactory  for  most  calculations.  0  shows  variations 
of  up  to  lOOJ^.  This  could  no  doubt  be  improved  upon  by 
a  more  sophisticated  relation  than  the  last  of  equations 
(12).  However,  varies  by  more  than  three  orders  of 
magnitude,  and  thus  the  proposed  relation  offers  con¬ 
siderable  improvement.  It  should  be  noted  that  over 


a  large  part  of  the  range  of  figures  1-4 


0  « 


A  B 

v_^A  +  v_3 


so  that  these  plots  are  fairly  insensitive  to  variations 


in  0. 


Conclusions 

It  is  felt  that  the  equations  presented  for  the 
macroscopic  description  of  ionic  mass  transport  offer  a 


Section  D.,  Project  IV. 


Page  141. 


rational  framework  for  the  correlation  of  the  transport 
properties  of  electrolytic  solutions.  The  relations  pro¬ 
posed  for  the  description  of  the  concentration  dependence 
of  the  ’s  are  rather  simple  and  can  no  doubt  be  improved 
upon  through  a  more  sophisticated  consideration  of  the 
structure  of  liquids  and  molecular  interactions-  They  do 
allow  the  diffusion  coefficient  to  be  predicted  with 
sufficient  accuracy  for  solutions  where  measurements  are 
not  available. 


o  to> 
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NOMENCLATURE 


subsidiary  transport  properties  defined  so  as  to 
reduce  the  concentration  dependence. 


c  -  concentration  of  binary  electrolyte  =  ~ 

-  concentration  of  species  1  (mole/cm^). 

c?  *=  S  c^,  the  total  concentration  (Including  the  solvent). 
D  -  molecular  diffusion  coefficient  of  binary  electrolyte. 

-  diffusion  coefficient  of  species  1  (cm^/sec). 

(0  -  molecuj.ar  diffusion  coefficient  of  binary  electro¬ 
lyte  based  on  activity  gradients  (cm2 /sec). 

’’diffusion  coefficient”  describing  the  Interaction 
of  species  1  and  J  . 

e  -  magnitude  of  the  electronic  charge  (coulomb). 

e“  -  symbol  for  an  electron. 

f^  -  mean  molar  activity  coefficient. 

P  -  Paradav's  constant  ( coulomb /equivalent ) . 

1.  -  current  density  (amp/cra2). 

k  -  Boltzmann's  constant  (erg/’K) . 

m  -  molality  (moles  per  1000  grams  of  solvent). 

-  a  symbol  for  the  chemical  formula  of  species  1. 
n  -  number  of  electrons  transferred. 

^  -  flux  of  species  1  (mole/cm^-sec) . 

R  -  universal  gas  constant  ( joule/bole-®K) . 

Sj^  -  stoichiometric  coefficient  of  species  1. 

t?  -  transference  numbev  of  species  1  with  respect  to 
molar  average  velocity. 

t?  -  transference  number  of  species  1  with  respect  to 
the  solvent  velocity. 
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T  -  absolute  temperature  (®K) . 

O 

u^  -  mobility  of  species  1  {cm^-molo /Joule-sec) . 
y.  -  fluid  velocity  (cm/sec). 

y.  -  velocity  of  species  1  (an  average  velocity  for  the 

species,  not  the  velocity  of  individual  molecules). 

-  molar  average  velocity  (cm/sec). 

z.  -  valence  or  charge  number  of  species  1  (negative 

for  anions ) . 

y  -  mean  molal  activity  coefficient. 

K  -  conductivity  (mho/cm). 

|i  -  viscosity  (poise). 

-  chemical  potential  of  binary  electrolyte  ( joule/mole) . 

p,.  -  the  electroc^amlcal^  potential  of  species  1 
( Joule/mole) . 

v_^  -  number  of  cations  per  molecule  of  electrolyte. 
v_  -  number  of  anions  per  molecule  of  electrolyte. 

V  =  +  v_. 

$  -  electrostatic  potential  (volts). 
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CONDUCTIVITY  CORRELATION 
FOR  KOH  AT  25°C 

Low  conceniration  range 


^/v  Co 

Experimental  points  based  on  conductivity  and 
Interpolated  values  of  viscosity  and  density.  FIGURE  I 
The  curve  is  the  best  least  squares  straight 
line  through  the  data.  Slope  =  6.646  x  10“°cm. 

Intercept  =  9.627  x  10-8cin. 
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CONDUCTIVITY  CORRELATION 
FOR  KOH  AT  25°C 

High  concentration  range 


Experimental  points  based  on  conductivity  and 
interpolated  values  of  viscosity  and  density. 

The  curve  is  the  beat  least  squares  straight  FIGURE  2 
line  through  the  data.  Slope  -  6.646  x  10"°cn». 

Intercept  «*  9.627  x  10“t^cm. 
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CONDUCTIVITY  CORRELATION 
FOR  KCI  AT  25"C 


Low  concentration  range 


^/v  Cq 
/CQ 

Experimental  points  booed  on  concructivi ty  and 
Interpolatea  valuea  or  viscosity  and  density.  rir^iiDE  o 
The  curve  Is  the  ijost  least  )auo.r*t»s  straight  rlLfUKt  J 


line 
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CONDUCTIVITY  CORRELATION 
FOR  KCI  AT  25°C 

High  concentration  range 


Co 

Experimental  points  based  on  conductivity  and 
Interpolated  value::  of  viscosity  and  density. 

The  curve  la  the  best  least  squares  straight 
line  through  the  data.  Slope  =  11.766  x  10"8cm. 
Intercept  =  39.710  x  10-8cni. 


FIGURE  4 
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Diffusion  coeflicienf,  cm^/sec 

Concentration  of  KCl,  moles/liter 


f/GURE  5 
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PROJECT  V. 

Diffusion  Coefficient  and  Solubility  of 
Oo  in  Aqueous  KOH  Solutions 


Introduction 

The  dynamic  behavior  of  gas  electrodes  is  strongly  depen¬ 
dent  on  the  solubility  and  diffusion  coefficient  of  the  reacting 
gas.  These  properties  which  greauly  influence  the  transport 
process  should  be  known  with  a  reasonable  degree  of  accuracy 
over  a  wide  range  of  temperature  and  electrolyte  concentrations. 
Of  particular  Interest  is  the  case  of  oxygen  in  aqueous  potassium 

hydroxide  solutions.  A  search  of  the  literature  produced  few 

1-4 

data  concerning  the  solubility  of  oxygen  in  basic  media  .  No 

comprehensive  data  on  the  diffusion  coefficient  oi  oxygen  in 

aqueous  KOH  solutions  was  found. 

The  solubility  of  oxygen  in  distilled  water  has  been 

investigated  by  many  workers  at  25® C,  and  it  is  now  generally 

_ 

agreed  that  the  value  is  1.26  x  10  tl^  mole  per  liter. 

~  S  2  6  “  9 

Values  of  1.8  to  2.6  x  10  cm  /sec  have  been  reported  for 

the  diffusion  coefficient  In  distilled  water.  Other  workers 
have  recently  investigated  the  solubility  of  oxygen  in  potas¬ 
sium  hydroxide  solutlons^^^ .  However,  one  finds  significant 
inconsls  -encies  when  comparing  the  results  of  the  various  ex¬ 
perimenters,  especially  xn  regard  to  ^.oncentrated  caustic 
solutions . 

In  this  laboratory  two  methods  were  utilized  in  making  the 
solubility  measurements.  In  the  first  a  Van  Slyke  manometer 
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apparatus  v/as  used  .^n  v/nlcii  dissolved  gas  vv'as  boiled  off  Prom  a 
presaturated  solution  and  the  pressure  of  the  gas  measured  at 
a  known  volume-  The  second  method  Involved  bringing  a  known 
amount  of  gas  into  contact  vfith  degassed  solution  and  measuring 
the  amount  absorbed. 

The  experimental  methods  chosen  for  determining  the  diffu¬ 
sion  coefficient  involved  the  measurement  of  transport  rates  to 
an  electrode  at  which  oxygen  was  consumed.  Since  it  was  felt 
that  the  equation  describing  transport  to  the  electrode  should 
be  unambiguous  and  free  of  empirical  constants,  a  rotating  disk 
electrode  and  a  stagnant  diffusion  cell  were  employed,  rather 
than  the  dropping  mc-'rcury  electrode. 

0^  solubility  measurements 

Experimental ;  The  Van  Slyke  manometer  was  used  for  the 
major  portion  of  the  solubility  values  reported.  The  apparatus 
is  shown  In  Figure  1  The  experimental  procedure  consisted 
of  admitting  a  kno\'/n  amount  of  liquid  through  stopcick  A  into 
chamber  E  which  was  initially  filled  with  mercury.  After  all 
of  the  liquid  was  in  the  chamber,  stopcock  A  was  ^.losed  and  the 
mercury  level  dropped  to  the  50  ml  level.  The  sample,  now  under 
vacuum,  was  agitated  by  a  magnetic  stirrer  in  order  to  facili¬ 
tate  the  boiling  off  of  the  dissolved  gases.  After  a  suitable 
period  of  agitation,  the  liquid  level  was  raised  to  the  2  ml 
mark  which  confined  the  gas  to  a  precalibrated  voluie.  The 
pressure  of  the  evolved  gas  was  determined  by  the  height  of  the 
mercury  in  the  manometer  column  C.  Tne  gas  was  ejected  and  the 
liquid  returned  to  the  50  ml  mark.  After  rebelling,  a  second 


Section  D.j,  Project  V. 


Page  154 . 


pressure  reading  was  taken.  The  difference  between  the  two 
presbares  represented  the  pressure  of  the  dissolvea  gas  at  the 
calibrated  volume. 

The  absorption  apparatus  and  experimental  procedure  are 

12 

described  in  detail  by  Kabatake  and  Hildebrand  and  will  not 
be  discussed  here. 

The  solutions  for  the  Van  Slyke  experiments  were  prepared 
by  bubbling  99.6^  pure  oxygen  through  the  sample  for  15  minutes. 
The  sample  was  in  a  constant  temperature  bath  controlled  to 
±0.2®C  during  this  period.  The  gas  was  presaturated  with  water 
vapor  by  first  bubbling  the  oxygen  through  a  solution  of  the 
sacne  temperature  and  composition  as  the  sample.  After  a  2-5 
minute  period  during  which  the  bubbling  was  stopped  and  an 
oxygen  blanket  was  maintained,  the  sample  was  withdrawn  and 
analyzed  for  oxygen  content. 

For  the  absorption  experiments  two  methods  were  tried  in 
preparing  the  degassed  solution.  In  the  first,  KOH  solutions 
were  boiled  for  24  hours  under  vacuum  while  boiling  and  freezing 
the  liquid  sample.  A  single  stage  vacuum  distillation  was  used 
in  the  other  procedure.  Distilled  water  v:as  frozen  over  pel¬ 
lets  under  vacuum  and  the  resulting  solution  boiled  and  cooled 
to  0®C  under  vacuum  to  iiisure  a  gas-free  sample. 

The  principal  factors  affecting  the  solubility  rtsuj.ts 
utilizing  the  Van  Slyke  apparatus  were:  a)  measurement  of  the 
sample  volume;  b)  measurement  of  the  pressure  and  volume  of  the 
gas  boiled  off;  and  c)  the  degree  of  saturat'^on  of  the  sample 
introduced  into  the  manometer.  Sample  volumes  were  measured 
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with  a  maximum  probable  error  of  ±1.0^,  the  measurement  of 
the  gas  volume  and  pressure  Introduced  a  maximum  probable 
error  of  0.5^  with  a  maximum  precision  of  3  x  1C  mole  gas 
per  liter. 

Supersaturation  was  believed  to  bo  a  source  of  error  since 
the  samples  v/ere  prepared  by  bubbling  gas  through  the  liquid. 
Experiments  were  conducted  In  wn^ch  samples  were  taken  from 
saturated  liquid  and  the  solubility  determined  as  a  function 
of  time  after  the  bubbling  was  stopped.  After  two  to  three 
minutes,  the  solution  had  developed  its  equilibrium  state  and 
no  change  In  the  solubility  could  be  detected  within  the  limits 
of  sensitivity  of  the  Van  Slyk  apparatus.  It  was  also  sug¬ 
gested  that  the  hydrostatic  head  could  cause  a  supersaturation 

13 

condition;  however,  Livingston,  Morgan,  and  Richardson 
showed  that  with  a  ten  centimeter  high  flask  the  height  of  the 
column  would  have  little  effect  on  the  solubility.  This  is 
especially  true  when  agitation  occurs  in  the  liquid  which  was 
the  case  In  these  experiments . 

Livingston  £t  aQ  also  showed  that  Henry's  law  is  obeyed 
for  oxygen  dissolved  In  water.  In  order  to  compare  our  results 
with  those  by  other  experimenters,  v/ho  present  their  data  at 
a  partial  pressure  of  0^  of  760  mm  Hg,  we  assumed  that  Henry's 
law  applies  to  solutions  of  KOH. 

Kabataka  and  Hildebrand  stated  that  the  largest  source  of 
error  In  their  experiments  was  the  degassing  process.  Substan¬ 
tial  Improvement  was  made  by  our  addition  of  the  vacuum  distil¬ 
lation  to  the  degassing  process. 
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The  error  In  al  sorption  determinations  was  estimated  to  be 
leas  than  1%}  however,  as  with  the  Van  Slyke,  this  error  in¬ 
creases  as  one  Increases  the  KOH  concentration.  The  solubility 
drops  by  a  factor  of  100  as  t:he  KOH  concentration  Increases 
from  0  to  10  molar,  thus  approaching  the  limits  of  the  apparatus. 
Results  and  discussion 

The  results  of  the  Van  Slyke  experiments  are  presented  in 
Figure  2 .  This  plot  shows  the  variation  of  oxygen  solubility 
as  a  function  of  potassium  hydroxide  concentration  at  a  total 
pressi  .'•e  of  one  atmosphere.  At  25”C  and  a  partial  pressure  of 
760  mm  Hg  oxygen  this  variation  is  represented  by  the  equation 

Log  S  =  -  0.1746  C  +  Log  1.26  x  10"^ 
where  S  =  the  solubility  of  oxygen  in  gram-moles  per  liter  and 
C  =  the  concentration  of  KOH  in  gram-moles  per  liter. 

Tne  linear  variation  of  the  logarithm  of  oxygen  concentra¬ 
tion  with  KOH  concentration  is  typical  of  the  ’’salting  out”  of 
a  non-electrolyte  in  an  electrolyte  solution.  The  vapor  pres¬ 
sure  of  the  solvent  changes  by  15  mm  Hg  at  25®C  over  the  range 
of  KOH  concentration  investigated;  therefore  the  partial  pres¬ 
sure  of  the  non -electrolyte  is  very  nearly  constant.  At  0®C 
the  difference  in  partial  is  negligible  and  at  60® C  the  change 
is  in  the  order  of  100  mm  Hg. 

The  theoretical  explanation  of  the  "salting  out"  (or  in 
some  cases  "salting  in")  of  a  non-electrolyte  by  an  electro- 
lyte  has  been  attempted  by  several  workers'^  '  ,  In  all 

cases  assumptions  must  be  made  which  cannot  be  expected  to  be 
valid  over  a  wide  range  of  salt  concentrations.  With  the 
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excectlon  of  tne  theory  of  McDevlt  and  Long,  the  assumption 
is  made  that  the  on-electrolyte  affects  the  dielectric  con¬ 
stant  of  the  s  elution.  This  variation  is  considered  to  be  the 
principal  effect  of  the  non-electrolyte.  By  assuming  the 
dielectric  constant  is  a  linear  function  of  the  non-electrolyte 
concentration,  a  Icgarltlimlc  "salting  out"  behavior  is  pre¬ 
dicted.  The  theory  of  McDevit  and  Long  is  based  on  the  concept 
of  an  "Internal  pressure"  of  the  solution. 

In  none  of  the  available  theories  are  the  parameters  suf¬ 
ficiently  well-knovm  to  allow  verification.  The  temperature 
variation  of  the  parameters  is  even  less  precisely  known. 
However,  for  most  cases,  experimental  measurements  show  that 
the  slope  of  the  log  concentratlor.  non-electrilyte  versus 
salt  concentration  plot  decreases  with  increasing  temperature. 

After  the  experiments  with  the  Van  Slyke  were  completed 
it  was  noted  that  there  were  large  differences  between  the 

4 

results  of  this  work  and  those  presented  by  Posplsll  and  Luzny  . 
Although  their  value  agreed  with  ours  for  distilled  water, 
the  discrepancy  was  in  the  order  of  200^  in  4  molar  KOH  solu¬ 
tion.  Results  obtained  with  the  Hildebrand  apparatus  are 
plotted  on  Figure  3  for  25®C.  Included  in  the  figure  are 
points  from  other  investigators  as  well  as  results  from  the 
Van  Slyke  experiment  corrected  to  a  partial  pres  sure  of  Og  of 
760  mm  Hg. 

10 

Walker  e^  ^  used  a  gas  chromatographic  method  for  their 
measurements  in  which  the  precision  is  said  to  be  of  the  order 
of  ±l^i  however,  because  of  the  very  low  oxygen  solubility  in 
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concentrated  KOH  solutions,  they  say  these  results  are  probably 
only  accurate  to  ±15-205^  for  KOH  concentration  higher  than 
7  molar. 

1"* 

Knaster  and  Apel’baum  have  recently  investigated  the 
solubility  of  oxygen  in  KOH  solutions  at  2^®C,  45® C,  and  75® C. 

An  apparatus  similar  to  the  Van  Slyke  was  used  for  their  experi¬ 
ments.  As  in  the  pre<='''nt  work,  the  method  of  saturation  was 
bubbling  oxygen  through  KOH  solutions.  Solubility  values  for 
4,  7,  and  10  molar  at  25®C  were  Interpolated  from  their  data 
and  sho;m  in  Figure  3. 

4 

Posplsil  and  Luzny  estimate  their  experimental  error  to 

be  less  than  0.9$^;  however,  their  results  do  not  follow  the 

trend  seen  by  other  investigators. 

Diffusion  coefficient  measurements; 

Theory  of  experimental  methods:  The  rotating  disk  elec- 

18 

trode  has  been  considered  by  Levlch  in  detail.  Gregory  and 
19 

Rlddlford  have  pointed  out  an  error  introduced  by  an  approxl- 

20 

matlon  in  the  Levich  solution.  Newman  has  provided  a  correc¬ 
tion  to  the  Levich  statement  which  is  more  exact  than  that 


given  by  Gregory  and  Rlddlford.  The  Levich  expression  for  the 

2 

limiting  current  density  in  Amp  per  cm  is 


=  0.6205  Sc 


-2/3 


where  n  is  the  number  of  electrons  transferred  per  molecule, 

F  is  thp  Faraday,  is  the  bulk  concentration  in  moles 

X 

per  cm  ,  w  is  the  angular  velocity  in  radians  per  second,  v  is 
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the  kinematic  viscosity  -  and  Sc  is  the  Schmidt  number  (^) . 
The  Newman  correction  which  Is  valid  for  the  range  400  <  Sc  < 
1400  is 


The  case  of  diffusion  through  a  stagnant  liquid  to  an  elec¬ 
trode  has  been  the  subject  of  investigations  by  Stackelberg^^, 
Laltlnen  and  Kolthoff  ^  and  Lingane^^.  The  absence  of  natural 
convection  is  essential  to  this  method.  Even  though  this  con¬ 
dition  may  be  difficult  to  realize,  the  low  current  densities 
used  in  stagnant  diffusion  cells  make  them  attractive  in.  the 
cases  where  activation  polarization  prevents  determination  of 
limiting  current  densities  in  the  forced  convection  regime. 
Diffusion  in  a  pore  to  a  flat  electrode  appears  to  provide  the 
best  opportunity  to  eliminate  natural  convection. 

Although  chronopotentlometry  has  been  demonstrated  to  be 
a  powerful  analytical  tool,  it  was  felt  that  for  work  with  the 
oxygen  electrode  constant  potential  operation  offered  special 
advantages.  In  chrcnopotentlometrlc  determinations  any  cur¬ 
rent  which  is  not  due  to  diffusion  of  the  species  under  inves¬ 
tigation  will  increase  the  transition  time  and  thus  the  appar¬ 
ent  diffusion  coefflcler t  will  be  incorrect.  In  the  case  of 
constant  potential  operation  such  effects  are  not  as  Impor¬ 
tant  as  long  as  the  concentration  of  the  diffusing  species 
eventually  become  zero  at  the  electrode. 

The  equation 


_ 0.6205  Sc~^/^ _ 

+  0.298  Sc’V^  +  0.14514  Sc‘^/^ 
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solved  for  the  boundary  conditions 

I.  C  =  ^bulk  ^  t  =  0 

II-  C  =  0  at  X  =  0  for  t  >  0 

III.  C  =  ^bulk.  X  00  for  t  >  0 
yields  the  solution 

1  =  nF  D  Id  (2) 

™|x=0  ^  Vt 

Thus  If  diffusion  In  the  pore  Is  unidirectional  and  the 
electrode  Is  operating  at  limiting  current,  the  current  should 
drop  as  t  .  The  correction  for  the  finite  length  of  the 
pore  Is  found  to  be  negligible  for  time  periods  of  several 
hours . 


If  a  first  order  rate  equation  Is  assumed  such  that 

1  =  K  C  ,  ^  ,  e^ 

electrode 

where  T]  Is  the  overpotential  and  the  cell  Is  to  be  run  at 
constant  overpotential,  boundary  condition  II  becomes 


NF  D  dC 


K  e  '  K  e 


_ 


at  x=0  for  t>0. 


A  typical  solution  of  equation  (l)  with  boundary  condition 
II  Is  plotted  In  Figure  4.  It  can  be  seen  that  the  current 
converges  to  that  predicted  equation  (2).  Therefore  even  If 
the  electrode  Is  not  at  llmi^..ng  current  when  the  potential 
is  first  applied.  It  does  converge  on  equation  (2)  regardless 
of  the  rate  constant  or  overpotential  applied. 


Experimental 


Disk  electrodes  were  constructed  as  indicated  In  Figure 
5.  The  disks  were  rotated  on  a  double  ball  bearing  shaft  by 


!  . .  n„i‘  ■  i||i~  . . . 
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a  1/4  horsepower  D.C.  motor.  The  speed  was  varied  from  900  rpm 
to  3600  rpm  and  measured  with  a  stroboscopic  tachometer. 

Current  was  provided  by  a  Lambda  Electronic  Model  28  constant 
current  supply  and  measured  by  means  of  the  voltage  drop  across 
a  calibrated  precision  resisror.  Electrode  potentials  were 
measured  relative  to  a  saturated  mercury /taercury  oxide  re¬ 
ference  electrode  and  recorded  by  a  Bausch  and  Lomb  V.C.M-5 
recorder. 


Platinum,  silver,  nickel,  amalgamated  silver  and  gold 
plated  copper  were  used  as  working  electrodes.  The  counter 
electrode  was  a  platinum  screen.  The  diameter  of  the  rotating 
disk  was  measured  by  a  comparator  to  an  accuracy  of  ±0.5^. 

The  cell  used  in  the  stagnant  diffusion  experiments  is 
shown  in  Figure  6.  The  working  electrode  and  counter  elec¬ 
trode  were  silver.  The  diameter  of  the  diffusion  pore  was 
measured  with  mlcrometered  wire  to  an  accuracy  of  ±1.5^. 
Current  and  potential  measurements  were  made  in  the  same  man¬ 
ner  as  the  disk  experiments  with  an  accuracy  of  ±0.55^. 

All  solutions  were  prepared  from  distilled  water  and 
reagent  grade  potassium  hydroxide.  99. 6^^  pure  oxygen  was 


used  saturating  the  solution.  Viscosity  of  the  solutions  at 
25® C  was  taken  from  the  data  of  Mcllhenny*^  and  from  the 

p5 

data  of  the  Solvay  Technical  Service’  for  60®C. 


The  electrodes  in  all  cases  were  pre treated  by  cathodic 
evolution  q±  hydrogen.  Experiments  v.  re  carried  out  in  a 
constant  temperature  bath  controlled  to  0.2®C.  Saturation  was 
achieved  by  sparging  oxygen  into  the  solution  through  a  fritted 
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glass  tube  for  15  minutes  prior  to  making  a  determination.  In 
some  cases  the  bulk  concentration  of  oxygen  was  determined  by 
a  Van  Slyke  analysis  to  Insure  that  saturation  was  being 
developed . 

In  the  rotating  disk  experiments,  the  limiting  current 
was  determined  at  3  to  4  rotational  speeds  with  both  silver  and 
plainum  electrodes.  These  limiting  currents  v;ere  plotted 
versus  the  square  root  of  the  rotational  speed  and  the  best 
straignt  line  fitted.  From  these  plots  and  the  reproducibility 
of  duplicate  runs,  the  error  in  the  determination  of  limiting 
current  was  estimated  to  be  of  the  order  of  3^.  The  diffusion 
coefficients  were  calculated  on  the  basis  of  the  Newman  correc¬ 
tion  to  the  Levlch  solution  for  the  transport  equation. 

The  current  versus  time  behavior  of  the  stagnant  diffusion 
cell  was  recorded  for  periods  of  15  to  30  minutes-  A  constant 
potential  between  1.0  and  1.2  volts  was  applied  to  the  cell 
during  each  experiment.  Before  applying  the  potential,  the 
liquid  in  the  pore  was  withdrawn  into  a  syringe  through  a 
hypodermic  needle,  thus  drawing  fresh  solution  from  the  bulk 
into  pore.  Five  to  ten  minutes  was  allowed  for  the  solution  to 
become  quiescent. 

In  general  the  experimental  curves  did  not  have  the  same 
shape  as  those  calculated  assuming  a  first  order  rate  equation, 
but  did  converge  to  a  straight  line  of  slope  -1/2.  The  diffu¬ 
sion  coefficient  was  calculated  using  equation  (2).  From  the 
scatter  of  the  points  on  the  current  versus  time  plot  the 
probable  error  in  the  diffusion  coefficients  was  estimated  to 
be  ±8^  in  the  high  concentration  range. 
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solutions,  the  activation  overpotentlal  was  low  and  the 
limiting  current  plateau  was  well  defined.  As  the  KOH  concen¬ 
tration  was  Increased  the  activation  overpotential  Increased 
and  the  limiting  current  plateau  became  less  distinct. 

When  platinum  was  the  metal  used  at  a  constant  current 

and  rotational  speed,  the  overpotential  was  found  to  vary  with 

2 

time.  As  was  noted  by  Palous  and  Buret  ,  the  potential 
depended  not  only  on  time  but  also  on  the  previous  current 
setting.  If  the  current  had  been  decreases  the  overpotential 
was  considerably  lower  than  if  it  had  been  Increased.  As  the 
concentration  of  KOH  was  Increased,  the  variation  of  potential 
with  time  Increased  until  in  2.0  molar  KOH  it  was  nearly  im¬ 
possible  to  determine  a  reproducible  current  versus  a  poten¬ 
tial  curve.  Therefore,  at  higher  than  2  molar  KOH  concentra¬ 
tion,  investigations  were  made  using  the  stagnant  diffusion 
cell . 

Walker  et  al^*^  have  recently  measured  the  diffusion  coeffi¬ 
cient  for  Og  in  KOH  solutions  at  25** C  by  the  polarographic 
technique.  The  calculated  dlffuslvlty  values  are  10-20;^ 

higher  than  those  obtained  in  this  laboratory.  The  modified 

27 

form  of  the  Ilkovlc  equation  which  employs  an  empirical  con¬ 
stant  was  used  to  evaluate  the  diffusion  coefficient  from  the 
polarographic  data.  Because  of  the  uncertainty  in  the  con¬ 
stant,  as  reported  by  Walker  et  al,  no  critl ',al  evaluation 
can  be  made  of  their  results. 
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Results  and  Discussion 

m 

The  diffusion  coefficient  for  25® C  and  60® C  Is  plotted  as 
a  function  of  potassium  hydroxide  concentration  in  Figures  1 ,  8. 

Because  of  the  rapid  increase  of  the  oxygen  diffusion  coeffi¬ 
cient  with  decreasing  KOH  concentration  In  dilute  solutions,  it 
Is  difficult  to  make  an  accurate  extrapolation  to  determine  the 
value  In  water.  From  Figure  7  It  may  be  estimated  to  be  be- 
tween  1.9  and  2.0  x  10  cm  /sec  at  25®C,  and  from  Figure  8  to 

be  between  4.6  and  4.8  x  10~^  cm^/sec  at  60® 0.  The  diffusivlty 

-5 

at  25® C  compares  well  with  values  between  1.9  x  10  and 
2.0  X  10~^  cm^/sec  reported  by  Hlmmelblau^^  and  2.12  x  10”^ 

9 

2  8 
cm  /sec  reported  by  Jordan  e^  al  . 

These  estimates  correspond  to  a  temperature  coefficient  of  ^ 

per  degree  for  the  diffusion  coefficient  In  water.  The 

7 

data  summarized  by  Millington  Indicates  a  temperature  coeffi¬ 
cient  between  2  and  3.5^  per  degree  In  the  temperature  range 
16®C  to  25®C. 

The  product  of  the  diffusion  coefficient  and  the  viscosity 

Is  plotted  in  Figure  9.  As  Is  noted  In  this  figure,  the  ae- 

crease  in  viscosity  with  Increasing  temperature  does  not  account 

for  the  temperature  coefficient  of  the  diffusion  coefficient. 

The  rapid  Increase  of  the  product  Dp,  In  dilute  solutions 

0 

was  also  reported  by  Jordan,  Ackerman  and  Berger  .  It  may  be 

I 

explained  by  assuming  a  cage-llke  structure  for  water.  As  the 
solution  becomes  more  concentrated  the  cages  become  filled 
and  diffusion  becomes  more  hindered.  ** 

With  silver  rotating  disk  electrodes  In  dilute  KOH 
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Solubility  of  Oxygen  in  Aqueous  KOH 
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Figure  6. 


Diffusion  Coefficient  of  Oxygen 
in  Aqueous  KOH  at  25  °C 
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Concentration  of  KOH,  moles/liter 

Figure  7. 


Diffusion  Coefficient  of  Oxygen 
in  Aqueous  KOH  at  60  °C 
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Concentration  of  KOH,  moles /liter 

Figure  8. 
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PROJECT  VI. 

Ionic  Mass  Transport  and  Current 
Distribution  in  Channels  at  High  Velocity  Plow 

Background  and  objectives 

From  a  number  of  studies  made  on  forced  convective  mass 
12  3 

transfer  *  *  It  has  become  clear  that  the  limiting  current 
technique,  l.e.,  measurement  of  cathodic  limiting  current  under 
the  assumption  of  diffusion-controlled  mass  transfer,  is 
eminently  suitable  for  obtaining  precise  mass  transfer  coeffi¬ 
cients  and  thus  establishing  generalized  mass  transfer  corre¬ 
lations.  In  a  number  of  cases  the  method  has  provided  a  check 
on  such  correlations  as  can  be  derived  by  theoretical  arguments, 

3 

in  the  case  of  particular  geometries  .  In  turbulent  mass  trans¬ 
fer  it  has  served  to  obtain  a  measure  of  the  validity  of  certain 
assumptions)  e.g.,  regarding  the  dependence  of  eddy  viscosity ,c.q. , 
diffuslvlty,  on  radial  coordinates,  in  an  otherwise  Independent 
development  of  analogy  between  momentum  -nd  mass  transfer^. 

Since  mass  transfer  coefficients  are  of  Immediate  practical 
Importance  in  chemical  and  electrochemical  engineering,  the 
limiting  current  technique  has  also  been  applied  to  natural 
convective  mass  transfer,  in  geometries  that,  softer,  defy 
analytic  solutions.  For  example,  the  correlation 

NUj^  =  0.13  (Grj^  Sc)V2 

can  be  mentioned,  which  was  found  to  adequately  represent  the 
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emplrloe.l  data  for  mass  transfer  by  natural  convection  be- 

4 

tween  two  horizontal  plates  . 

5 

In  the  work  undei'taken  by  Hickman  ,  a  first  attempt  was 
made  to  emplrlc.ally  disentangle  the  Influence  of  forced  and 
free  convection  In  superimposed  forced  and  free  convective 
flow.  The  case  studied,  mass  transfer  at  horizontal  electrodes, 
represents  a  situation  of  considerable  practical  Interest, 
where  analytic  development  has  had  fow  useful  results^.  The 
data  obtained  In  the  laminar  region  and  which  suggest  a 
separation  between  distinct  convective  regimes: 

Gr  «  >  920  Re ,  free  convective  mass  transfer 

d'^L 

Gr  2  <  920  Re^  forced  convective  mass  transfer 

d  L 

are  very  encouraging  as  to  the  continued  usefulness  of  limiting 
current  techniques  In  such  complicated  situations. 

At  the  same  time,  the  refinement  Introduced  by  using 
sectioned  electrodes  seems  particularly  advantageous  to  the 
empirical  solution  of  these  correlation  problems,  since  It  is 
necessary  to  distinguish  current  density  Inequalities  along 
the  electrode  which  are  caused  by  not-yet-developed  mass 
transfer  profiles,  from  those  current  density  Increases  which 
are  caused  by  the  onset  of  natural  -convective  Instabilities. 

With  the  limited  system  capacities  at  his  disposal  (5  QPM 
total  flow  rate,  14  amp  total  current),  Hickman  was  not  able 
to  pursue  his  Investigations  through  the  forced  convection 
transition  region  and  Into  the  turbulent  regime.  Apart  from 
the  capacity  factor,  there  are  also  a  few  complications  which 


1 


ir|;|jl' 


Section  D.,  Project  VI. 


Page  178. 


arise  from  the  nature  of  the  flow  situation  and  go  relatively 
unnoticed  in  the  laminar  range : 

(1)  Especially  beyond  the  laminar  regime  it  Is  necessary 
to  depend  on  available  semi -empirical  expressions  for  pure 
forced  convective  mass  transfer.  If  one  wants  to  "sift  out" 

any  free  convection  influences.  This  in  turn  makes  It  advisable 
to  employ  simple  geometries.  Accordingly,  the  presently 
available  wide  electrode  channel  (  cross  section  15  by  1  cm) 
is  assumed  to  approximate  an  Infinite  parallel-plates  geometry. 

(2)  Fnirthermore ,  to  be  able  to  distinguish  the  effect 

of  free  convection  on  the  initiation  of  turbulence  In  the  chan¬ 
nel,  It  Is  necessary  to  exclude  external  disturbances  of  the 
flow  as  much  as  possible  (specifically,  vibrations  of  the 
circuit  Induced  by  the  circulation  pump)  and  to  control  flow 
rates  closely- 

(3)  As  expected,  and  quantitatively  confirmed  by  Hickman, 
any  forced  convection  Introduces  a  non -homogeneous  boundary 
condition,  since  the  leading  cathode  section  reaches  limiting 
current  later  than  the  trailing  sections.  This  effect  is 
especially  noticeable  at  high  mass  transfer  rates  (high  He, 
high  Gr)  and  Is  somewhat  leveled  off  by  natural  convection. 

In  the  turbulent  regime  It  is  pronounced  but  Limited  to  the 
first  entrance  sections  of  the  channel. 

The  Inequality  can  be  overcome  by  using  local  cathodic 
overpotentials  to  Indicate  when  local  (section)  current  den¬ 
sities  reach  limiting  values.  Accordingly,  a  multiple 

5 

measurement  system  developed  by  Hickman  has  been  added. 


Section  D.,  Project  VI. 


Page  179. 


scanning  the  section  current  densities  and  corresponding 
section  overpoten tJ.als .  This  refinement  is  most  profitable 
when  one  takes  an  Interest  in  the  development  of  concentra¬ 
tion.  profiles  in  free  convective  (and  forced  convective) 
mass  transfer,  and  consequently  in  the  mechanics  of  the  develop¬ 
ment  of  natural  convection  flow  patterns. 

The  apparent  presence  of  reproducible  convection  patterns 
termed  "roll  cells"  seems  to  make  estimates  of  free  convection 
flow  pattern  dimensions  possible. 

Disregarding,  for  the  present,  any  theoretical  approach 
to  the  second  part  of  the  above  objectives,  i--  can  be  said  that 
continuation  of  experiments  on  combined  free  and  forced  convec¬ 
tion  poses  the  following  requirements  for  the  apparatus: 

(1)  The  system  should  allow  high  enough  flow  rates  to 
cover,  with  the  existing  electrode  cells,  the  range  Re^<10,000, 

(2)  The  power  source  for  the  system  should  be  able  to 
supply  correspondingly  high  currents. 

(3)  The  flow  rate  measurement  should  have  a  uniform  rela¬ 
tive  accuracy  over  the  range  of  flow  rates  covered. 

(4)  The  flow  regulation  should  be  precise  over  as  wide  as 
range  as  possible. 

(5)  The  velocity  profile  in  the  flow  cell  should  be  well 
defined  before  mass  transfer  starts. 

(6)  Freedom  from  external  vibrations  should  be  assured  as 
much  as  possible. 

(7)  Sectional  current  densities  euid  reference  electrode 


t 


potentials  should  be  scanned  sequentially,  to  allow  construction 
of  local  limiting  current  curves. 
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Moreover,  from  practical  considerations,  a  number  of 
additional  requirements  can  be  expressed. 

(8)  The  temperature  of  the  electrolyte  should  be  kept 
constant  to  within  1®C. 

(9)  It  should  be  possible  to  apply  current  to  the  elec¬ 
trodes  at  a  pre-set  rate  of  increase.. 

(10)  The  resistances  used  to  measure  the  section  currents 
should  be  adjustable  within  such  ranges  as  to  make  the  poten¬ 
tial  variation  between  the  sections  less  than  the  assumed 
criterion  of  2  mV. 

(11)  With  a  view  to  possible  related  studies,  employing 
a  modified  arrangement,  the  parts  of  the  flow  circuit  should 
be  easily  replaceable  or  interchangeable. 

(12)  An  adequate  protection  should  be  provided  for  the 
operator  in  case  of  a  catastrophic  leak  or  other  accident 
involving  the  circulation  system. 

Design  and  construction  of  equipment. 

Supporting  structure  and  protection: 

A.  Design  considerations 

It  was  felt  that  an  expansion  of  the  circulation 
circuit  made  it  necessary  to  provide  an  adequate  support  for 
its  elements,  especially  the  heavy  electrode  flow  cells  and 
the  necessary  valves  and  flowmeters.  The  weight  of  any  kind 
of  piping  and  the  liquid  hold-up  could  also  be  expected  to  be 
appreciable,  since  the  circuit  dimensions  are  conditioned, 
in  part  by  the  flow  cell  length  (approximately  6  ft),  in  part 
by  the  moderate  head  of  the  available  pump  (14  p.s.l.).  The 
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supporting  framework  would  have  to  stand  free  from  w?lls,  to 
give  maximum  accessibility.  It  could  also  serve  as  a  c'‘ntinuing 
structure,  such  that  it  would  be  easy  to  attach  shielding  on 
all  sides  to  protect  the  operator  of  the  system  against  acci¬ 
dental  leaks  of  the  corrosive  electrolyte. 

B.  Execution 

A  basic  supporting  structure  was  erected  lonslstlng 
of  8  Unistrut  Double -Reinforced  Columns,  mounted  between  floor 
and  celling  and  forming  a  4  by  2  rectangle  with  equal  spacing. 

A  separate  support  was  provided  for  mounting  the  heavy  circu.''a- 
tion  pump.  In  such  a  way  that  propagation  of  pump  vibrations 
through  the  frame  would  be  avoided.  All  heavy  parws  of  the 
circulation  system  are  supported  by  single  or  double  Unlstrut 
channels  mounted  onto  this  basic  structure.  Immediately  abo\?e 
the  floor,  covering  the  entire  area  of  the  rectangular  structure, 
a  dropping  tray  was  mounted  through  which  spilled  electrolyte 
is  drained.  On  the  long  sides  of  the  rectangle,  at  working 
height,  stainless  steel  frames  were  mounted.  In  each  of  which 
four  suspended  plexiglass  windows  can  slide.  The  windows  are 
39  In  high  by  30  In  wide  and  fabricated  from  1/4  In  thick  Luclte. 
All  area  under  and  above  the  windows  can  be  covered  with  rectan¬ 
gular  shields  of  stainless  steel  sheet.  Fixed  plexiglass  win¬ 
dows  mounted  onto  the  short  sides  of  the  structure  give  a  view 
of  the  Inside  over  the  whole  height.  All  valves  can  be 
operated  from  outside  the  shielding.  Connections  to  measuring 
equipment  outside  are  protected  against  corrosion  and  stray 
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Electrolyte  circulation  system; 

A.  Design  considerations 

Design  calculations  were  based  on  the  characteris¬ 
tics  of  the  available  Aurora  centrifugal  pump  (capacity  70  GPM, 
total  head  14  PSI).  Assuming  a  total  circv.it  lenguh  of  approxi¬ 
mately  10  m  and  employing  a  friction  factor  of  0.01,  the  Fanning 
friction  relation 

_  32  f  L  e  q"' 

dS  ■ 

gives,  for  water,  at  maximum  capacity  D  =  3.5  cm.  The  minimum 
holdup  by  the  circulation  system  in  that  case  is  approximately 
10  liter.  The  real  holdup  was  expected  to  be  greater,  because 
regulation  of  the  flow  rate  requires  a  by-pass  of  the  electrode - 
cell  and  because  the  anode  and  the  cathode  compartment  of  the 
flow-cell  have  to  be  fed  separately.  The  bl-sectlon  of  anode 
and  cathode  feed  line  had  to  be  designed  early  enough  to  in¬ 
clude  a  flowmeter  in  one  of  them  (presently  the  cathode  feed 
line).  This  was  thought  preferable  over  measuring  the  total 
electrode  flow  rate,  since  Imbalances  and  disturbances  could 
easily  be  expected  during  an  experiment.  To  avoid  the  cost  and 
complications  of  a  second  flowmeter  in  the  alternate  feed  line, 
outlets  for  differential  pressure  measurement  between  cathode 
and  anode  were  designed  at  the  entrance  of  the  electrode  flow 
channel.  The  design  of  the  circuit  tubing  would  have  to  be 
such  that  any  of  the  two  feed  lines  co  .jLd  be  connected  to  any 
of  the  two  electrode  compartments  in  vjhatever  electrode  orien¬ 
tation  with  respect  to  the  vertical  would  be  desirable.  This 
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could  be  achieved  either  by  employing  flexible  tubing,  or  by 
a  rigid  construction  employing  interchangeable  elements.  The 
more  costly  solution  of  constructing  a  four-way  combination 
valve  of  resistant  material  was  not  considered.  Finally,  the 
design  would  have  to  provide  drain  lines  and  a  storage  possi¬ 
bility  for  the  electrolyte  in  such  a  way  that  handling  the 
large  volumes  cf  corrosive  solution  would  not  present  great 
difficulties  or  risks.  Since  in  previous  experiments  tempera¬ 
ture  increases  of  several  degrees  centigrade  due  to  pump  energy 
dissipation  were  observed,  a  modest  cooling  facility  would  have 
to  be  included  somewhere  in  the  circuit,  preferably  in  the  main 
reservoir. 

B.  Execution 

It  was  decided  to  use  1  1/2  in  I.D.  pyrex  Double- 
Tough  Glass  Tubing  for  the  main  circuit  piping.  In  this  way 
sagging  connections  were  avoided  and  continuous  support  of 
long  lines  was  unnecessary.  Also,  flanged  connections  could  be 
employed.  The  glass  tubing  was  available  in  unit  elements 
of  various  length,  down  to  3  1/2  in.  Elements  are  connected 
by  Corning  type  flanges  and  sealed  by  Teflon  0-rings.  Connec¬ 
tions  to  metal  flanges  are  sealed  with  neoprene  gaskets.  The 
inlet  and  outlet  of  the  electrode  flow-cells,  as  well  as  all 
valves  and  flowmeters  in  the  system,  were  provided  with  1  1/2  in 
I.D.  stainless  steel  flanges  of  150  lbs.  standard  ASME  speci¬ 
fications  (a  few  connections  were  1  in  I.D.). 

Regulation  of  flow  rates  through  the  return  line,  the 
cathode,  and  the  anode  was  provided  by  inserting  in  the  circuit 
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three  Powell  globe  valves  of  1  1^/2  in  I.D.  stainless  steel, 

150  lbs  ASME  specifications. 

Plow  rate  measurements  ore  made  by  a  magnetic  flowmeter, 
which  offers  the  advantage  of  minimal  pressure  loss  and  high 
accuracy  over  a  wide  range  of  flow  rates,  through  adjustment 
of  the  full-scale  range  between  5  and  60  QPM.  The  transmitter 
signal  is  converted  to  a  current  and  continuously  recorded. 

For  extremely  low  flow  rates  rotameter  measurement  can  be  pro¬ 
vided  by  a  circuit  extension. 

Propagation  of  piamp  vibrations  through  the  circuit  piping 
is  prevented  as  much  as  possible  by  flexible  connections  to 
pump  inlet  and  outlet. 

The  drain  lines  of  the  circulation  system  consist  of  two 
1  in  I.D.  branches  at  the  lowest  points  of  the  circuit  (on  the 
suction  and  on  the  pressure  side  of  the  circulation  pump) . 

They  lead,  through  1  in  I.D.  Edwards  globe  valves  of  stainless 
steel,  to  1/2  in  polyethylene  tubing  and  a  small  Cole -Parmer 
polyethylene  pump  which  can  pump  the  electrolyte  into  one  of 
two  Nalgene  storage  bottles  of  15  gallon  capacity,  located 
above  the  system.  These  storage  reservoirs  can  empty  into  the 
main  reservoir,  likewise  a  Nalgene  vessel  of  15  gallon  capacity. 
The  main  reservoir  can  be  cooled  by  means  of  copper  coils 
through  which,  depending  on  the  need,  tap  water  or  refrigerated 
water  is  circulated. 

Figure  1.  gives  a  schematic  representation  of  the  elect’^o- 
lyte  circulation  system,  in  the  form  of  a  block  diagram. 
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Electrode  flow-cell. 


A»  Design  considerations 

The  electrolyte  enters  the  electrode  channel 
coming  out  of  a  tube  of  circular  cross  section  and  3.75  cm  I.D. 
It  has  to  conform  Itself  to  a  rectangular  channel  of  15  by 
1  cm  cross  section.  The  distribution  device  used  for  the  wide 
electrode  channel  proved  to  be  unsatisfactory  at  Re  numbers 
beyond  the  laminar  range.  It  was  thought  necessary  to  have  a 
completely  new  inlet  section  constructed,  which,  like  the  old 

5 

one  ,  could  be  flanged  onto  the  main  flow-cell.  In  view  of 
the  lack  of  experience  with  entrance  behavior  at  high  flow 
rates,  it  seemed  advisable  to  provide  an  excess  of  flow 
straightening  devices,  with  the  possibility  to  remove  one  or 
more  of  them  if  the  resulting  pressure  drop  would  prove  to 
diminish  the  flow  range  capacity  of  the  cell  too  much.  It  was 
not  considered  necessary  to  revise  the  outlet-section  of  the 
cell,  except  for  an  enlargement  of  its  capacity. 

As  mentioned  already  in  the  previous  section,  outlets  for 
differential  pressure  measurements  had  to  be  made  in  the  en¬ 
trance  of  the  cathode  and  the  anode  channel.  The  pressure  drop 
along  both  channels  was  considered  of  Interest  for  the  inves¬ 
tigation  and  therefore  similar  outlets  were  designed  at  the 
end  of  the  two  channels. 

B.  Execution 

A  new  inlet  section  was  constructed  out  of 
Luclte.  Each  circular  inlet-channel  branches  symmetrically 
into  2,  then  4,  then  8,  finally  16  circular  channels,  with 
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calming  sectlona  of  proportional  length  inserted  between  the 
bi-sectlon  points.  The  total  cross  section  of  the  channels 
stays  constant  throughout.  The  last,  sixteenfold,  stage  emp¬ 
ties  into  a  short  length  of  rectangular  channel  of  dimensions 
15  by  1  cm.  Here  the  fluid  expands  slightly  in  the  width  dimen 
slon.  The  channel  is  then  ended  with  a  15  mesh  stainless  steel 
wire  screen,  which  can  be  taken  out.  Also,  the  Lucite  elements 
in  which  the  8-  aiid  the  16-fold  division  of  the  flow  tedces 
place,  are  not  glued  to  the  first  and  last  elements  but  held 
in  place  and  sealed  by  neoprene  O-rings  under  pressure  from 
the  inlet-section  flanges.  They,  too,  can  therefore  be  removed 
to  alleviate  the  pressure  drop. 

Between  the  Inlet  section  and  the  main  electrode  channel 
an  additional  calming  section,  of  15  cm  length,  can  be  in¬ 
serted  by  means  of  flanges.  It  is,  in  turn,  ended  with  a 
removable  stainless  steel  wire  screen  of  15  mesh. 

Pressure  reading  outlets  were  drilled  in  the  short  sides 
of  the  cathode  and  anode  channel,  at  7  cm  and  5  cm,  respec¬ 
tively,  both  downstream  from  the  inlet-section  and  upstream 
from  the  outlet-section.  The  holes  were  1/4  in  diameter 
throughout  the  wall  thickness,  except  for  the  last  1/16  in 
length  opening  into  the  channel  which  had  1/16  in  diameter. 

The  actual  connection  to  the  mercury  U-tubes  was  a  1/8  in 
diameter  branch.  Just  above  this  point.  This  connection 
could  be  severed  by  lowering  a  plug  into  the  main,  1/4  in 
diameter  hole.  The  plug  was  machined  to  fit  the  1/16  in 
diameter,  flush  with  the  channel  wall. 
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To  avoid  great  lengths  of  tubing  necessary  for  differen¬ 
tial  pressure  measurement  between  inlet  and  outlet  of  the  same 
channel,  it  was  decided  to  measure  the  absolute  pressure  at 
these  points,  by  means  of  an  open  mercury  tube.  The  differ¬ 
ent."'  al  pressure,  hov;ever,  is  measured  between  the  two  charjiels 
at  the  inlet  and  outlet  of  the  flow-cell. 

Finally,  the  cathode  holder  was  provided  with  a  groove 
and  a  neoprene  0-rlng  for  seal,  to  avoid  any  degree  of  contamina 
tlon  of  the  electrolyte  by  grease. 

Power  supply: 

A.  Considerations 

A  new  power  supply  had  to  be  designed,  capable 
of  supplying  currents  of  up  to  100  amp  as  a  current  ramp. 

Because  of  the  necessity  of  trying  out  several  increase  rates 
in  mass  transfer  situations  where  the  limiting  current  curve 
has  no  conspicuous  plateau,  the  current  source  should  be  able 
to  reach  any  maximum  in  a  time  that  could  be  pre-set.  A  maxi¬ 
mum  voltage  should  not  be  exceeded  during  this  process. 

B.  Execution 

The  system  designed  to  meet  these  requirements 
is  shown  in  the  block  diagram.  Figure  2  .  A  Kintel  Operational 
Amplifier  connected  as  an  analog  integrator  is  used  to  generate 
the  rsimp  function.  Current  maximum  can  be  set  at  any  va3ue 
between  0  and  100  amp;  time  ranges  are  1,  2,  5,  10  and  20 
minutes.  Current  stab'  .ioy  is  ,1%  of  range;  linearity  is  .1^. 

A  maximum  voltage  of  18  V  is  not  exceeded. 

The  transistor  bank  in  the  power  supply  Is  cooled  by 


water. 
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Electric  measurement  system: 


A.  Considerations 


At  higher  total  currents  applied  to  the  cathode, 
it  was  thought  probably  that  the  variable  resistances  (combined 
in  the  current  shunt  box)  which  serve  to  measure  sectional 
currents,  might  have  too  little  variability  to  allow  equal 
IR-drops  between  the  cathode  section  contactors,  and  the  com¬ 
mon  busbar,  wMch  is  connected  directly  to  the  power  source. 
This  would  be  caused  by  the  fixed  resistance  presented  by  the 
leads  and  12 -pin  connections  in  the  cathode  holder,  as  well  as 
the  slide-wire  resistance  and  .several  contact  resistances  in 
the  shunt  box.  A  calibration  of  the  total  resi^tarc-^,.  nvolved 


showed  that,  indeed,  assuming  uniform  current  density  over  the 
cathode,  deviations  from  equipotentlality  of  the  central  sec¬ 
tions  would  amount  to  more  than  2  mV  (the  ass’  ed  criterion) 
at  total  currents  above  approximately  30  amp.  The  situation 
would  seem  more  serious,  however,  even  at  10  amp,  with  respect 
to  the  equlpotentlal  condition  of  the  long  side-sections  (carrying 
each  4/15  of  the  total  current,  regardless  of  current  density 
distribution  along  the  cathode).  Deviations  of  more  than  20 
mV  were  to  be  expected,  with  respect  to  the  central  sections. 

It  seemed  advisable,  as  a  provisional  solution  for  low-current 
measurements,  to  drastically  diminish  the  fixed  resistance  of 
the  connections  to  the  side -sections. 

To  enable  a  continuous  measurement  of  (l)  total  current; 

(2)  section  currents;  (3)  local  overpotentials;  (4)  electro¬ 
lyte  flow  rate,  an  Integrated  electric  system  had  to  be 
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designed,  which  could  be  compactly  housed  In  one  or  two  cabi¬ 
nets,  within  the  operator’s  reach.  S  ich  a  system  could  be 
envisaged  as,  either  leading  to  simultaneous  recording  of 
all  data  mentioned  by  means  of  single  and  multi -point  recorders 
or  as  a  rapid  sequential  scanning  of  all  signal -points  by  a 
scannei'-digltal  voltmeter  combination,  the  output  from  which 
could  be  printed  and  inspected,  or  directly  punched  Into 
IBM  cards  to  facilitate  possible  computer  treatment  of  the 
data. 

B.  Execution 

The  current  shunt  box  was  modified  by  taking  out 
the  variable  resistance  wire  corresponding  to  the  (combined) 
side -sections  and  replacing  it  by  two  fixed  copper  wires  of 
1  raft  resistance,  eachlng  leading  directly  to  one  of  the  slde- 
sectlons.  The  total  contactor-torbusbar  resistance  of  the 
side-section  connections  Is  now  1.7  mci,  which  is  expected  to 
make  the  current  shunt  box  usable  up  to  total  currents  of 
approximately  30  amp.  A  decrease  of  total  resistance  has  also 
been  achieved  for  the  other  leads  by  eliminating  the  12 -point 
plug  connection.  The  current  shunt  box  has  been  fastened  to 
the  cathode  holder. 

An  Integrated  system  leading  to  simultaneous  recording 
of  most  of  the  data  Involved  has  been  designed  and  constructed 
(see  block  diagram.  Figure  3.).  Such  a  system  can  without 
effort  be  simplified  to  suit  a  rapld-scannlng-dlgital  volt¬ 
meter-printer  combination. 

The  present  system  combines: 
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(1)  a  pajiel  of  voltage  dividers  which  make  it 
possible  to  scale  each  sectional  current  output  in  such  a  way 
that  the  recorded  signal  is  proportional  to  the  current  den¬ 
sity  on  the  sectlonj 

(2)  a  stepping  switch,  actuated  by  the  10-point 
recorder,  which  selects  the  sectional  current  density  output 
signals  for  sequential  recording) 

(3)  a  microvolb-ampllfler,  with  ranges  between 
50  and  2000  m-V,  v;hich  serves  to  amplify  low  current  section 
current  density  signals) 

(4)  a  selector-switch  which  det3rmlnes  which 

10  of  11  available  current  density  signals  will  be  recorded;. 

(5)  a  ten-point  Leeds  and  Northrup  x-y  recorder, 
which  reglstrates  section  current  density  against  total 
current) 

(6)  a  programming  pin  board  which  combines  the 
leads  originating  in  the  9  reference  electrodes  with  branches 
from  the  11  leads  which  contact  the  variable  resistance  wires 
at  the  outlet  of  the  current  shunt  box; 

(7)  ein  attenuator  to  bring  any  signal  from  a 
chosen  combination  in  (6)  down  to  the  level  required  by  the 
recorder; 

(8)  an  x-y-y  recorder  which  reglstrates  the 
chosen  overpotential  signals  against  total  current; 

(9)  an  x-y  recorder  to  reglstrate  the  flow  rate 
signal  from  the  magnetic  flowmeter  to  current  converter 
against  the  total  current; 
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(10)  a  shunt  Introduced  Into  the  power  source  cir¬ 
cuit  to  give  an  output  corresponding  to  the  Instantaneous 
current,  the  order  of  magnitude  of  which  is  kept  within  the 
mV-range  by  the  decade -^set ting  of  the  current  maximum. 

Literature  study. 

Forced  convective  mass  transfer;  The  results  obtained  by 

5 

Hickman  for  systems  considered  under  laminar  forced  convective 
control  show  a  remarkable  consistency  with  empirical  correla¬ 
tions  available  for  recwangular  ducts; 

Nu^  •  1.85  ReV2  gc^/^  (d/L)V3  (Norris  and  Streld) 

In  particular  the  current  density  profiles  for  the  cathode - 
over-anode  situation,  l.e..,  involving  indefinite  hydrodynamic 
stability,  are  in  perfect  agreement  with  boundary  layer  theory 
predictions  of  mass  transfer  rates  proportional  to 

A  similar  agreement  has  not  been  claimed  for  data  obtained 
in  the  range  2500<Re<6000,  which  were  supposed  to  be  consistent 
with  either  the  equation  for  fully  developed  mass  transfer  in 
pipes  (originally  developed  empirically  for  heat  transfer 

Nu^  =  0.027  Re°’®  SeV^  (Sleder  and  Tate) 

or  the  equation  for  the  entrance  region; 

Nu^  =  0.276  Re^*^®  (Sc  d/L)V3  (Shaw,  Reiss,  Hanratty®) 

It  is,  however,  doubtful  whether  the  narrow  range  and  the  tran¬ 
sition  character  of  tne  Re  numbers  concerned  can  be  considered 
representative  for  the  turbulent  regime.  It  should,  in  parti¬ 
cular,  be  considered  possible  that  natural  convection  intro¬ 
duces  an  earlier  transition  to  turbulent  flow.  Exactly  in  that 


Section  D.,  Project  VI. 


Page  19P 


case,  the  Influence  of  natural  convection  occurring  In  the 
flat,  rectangular  geometry  employed  would  se^^ra  to  contribute 
to  marked  anisotropy  of  turbulent  motion  such  as  to  make 
comparison  with  results  obtained  in  pure,  turbulent  forced 
convective  mass  transfer  deficient. 

In  the  range  of  Re  >  10,000,  however,  such  a  comparison 
would  be  Justified  and  the  quality  of  the  results  obtained  in 
correlating  laminar  forced  convective  mass  transfer  suggests 
that  the  limiting  current  technique  is  sufflclonbly  accurate 
to  provide  a  check  on  the  validity,  for  mass  transfer,  of  the 
sevei’al  refined  analogy  expressions  developed  for  turbulent 
mass  transfer^^®'^  as  well  as  the  familiar  Chilton-Colburn 


analogy 


St  =  (f/2)  Sc”^/^ 


The  most  Important  aspect  of  this  check  would  concern  the 
variation  of  the  dependence 

(3  iH  St/d  In  f) 

as  a  function  of  the  Sc -number.  Recent  experimf^nts  by  Liglitfoot 
and  Hubbard^®  ccnflrm  the  evidence  already  available  that  the 
Sc— dependence  given  by  the  Chilton-Colburn  expression  agrees 
better  with  the  data  than  that  Implied  in  Delssler’s  analogy^. 
Though  the  geometry  of  the  channel  employed  for  such  tests 
is  of  importance  in  the  laminar  range  of  Re  numbers,  it  appears 
that  in  the  range  indicated  above  the  results  for  mass  trans¬ 
fer  between  parallel  plates  do  no<-  deviat?  mere  than  from 
those  for  round  pipes  '  .  A  theoretical  development  of  the 
heat  transfer  to  friction  anr..lofcT  for  the  case  of  special 
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Interest  here,  transfer  from  one  wall  to  flow  between  parallel 
plates,  .3  aval lab le^^. 


The  sectioned  electrode  technique  enables  measurements  of 
entrance  effects  In  turbu!  mt  transfer.  Mass  transfer  tech¬ 
niques,  in  particular  electrolytic  methods,  require  less  compli¬ 
cations  In  this  regard  than  measurements  of  heat  transfer, 
where  thermal  Insulation  of  short  lengths  of  channel  presents 
a  problem.  Entrance  region  measureme  its  are  of  Interest 
especially  at  very  short  lengths  and  relatively  low  Sc  numbers, 
where  the  assumption  that  e/v,  l.e.,  the  rate  of  eddy  dlffu- 
slvlty  to  kinematic  viscosity,  depends  only  on  the  dimension¬ 
less  coordinates,  y+and  loses  credibility. 

Influence  of  free  convection;  instability  and  transition; 

5 

The  correlation  developed  by  Hickman  to  distinguish  betv/een 
regimes  of  laminar  forced-convective  mass  transfer  and  free- 
convective  mass  transfer*  employs  the  modified  Grashof  number: 

Gr  p  =  pgoAcd^L/v^  . 
d^I-  ® 

Though  this  dlmensj.onless  combination  represents  a  convenient 

way  of  summarizing  the  basic  finding.  It  dees  not  contribute  to 

Its  explanation  since  the  actual  value  of  the  mass  transfer 

boundary  l^er  Is  not  employed. 

Also,  Insofar  as  available  Instability  criteria  cam  be 

used  for  comparison.  It  should  be  noted  that  the  criterion 

12 

derived  by  Velte  ,  Ra  >  5300,  refers  only  to  square  cross 
section  and  utilizes  exclusively  the  equlvalert  diameter. 

The  application  of  Instability  criteria  and  the  analysis 
of  transverse  flow  patterns  with  regard  to  the  results  obtained 
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in  laminar  flow  through  a  rectangular  duct,  would  require 
first  a  careful  investigation  of  the  transition  region  to  tur¬ 
bulent  flow  under  conditions  of  hydrodynamic  stability,  i.e.,  i 

in  the  cathode -over-anode  situation.  Precise  control  over  flow 
rates  and  absence  of  any  gross  flow  disturbances  a.re  essential 
conditions. 

Influence  of  free  convection:  transverse  flow  patterns; 

Convection  patterns  termed  “roll  cells”  are  apparently  present 

in  laminar  flow  influenced  by  free  convection.  Though  their 

existence  was  not  demonstrated  conclusively,  it  has  been  made 

plausible  by  comparison  with  similar  phenomena  observed  in  heat  * 

15 

transfer  to  a  Horizontal  layer  of  air  in  shear  flow  . 

For  any  further  analysis  t  is  necessary  to  have  quantita¬ 
tive  and  reproducible  data  on  tne  characteristic  dimensions  of 
these  convection  patterns  and  their  point  of  Initiation.  The 
most  promising  approach  to  an  explanation  would  be  an  analysis 
of  similarity  with  the  Taylor -Goer tier  vortices  observed  ir 
laminar  flow  on  highly  concave  walls.  These  "neutral"  vor¬ 
tices  are  likewise  the  result  of  a  three-dimensional  distur- 
14 

bance 

Such  an  analysis,  however,  presents  extremely  complicated 
mathematical  problems  and  would  only  be  warranted,  as  well  as 
considerably  helped,  by  substantial  data  on  the  dimensions  of 
the  phenomena  observed. 

Since  the  localization  of  roll  cells  is  still  highly 

4 

dependent  on  the  measurement  of  deposition  profiles,  it  is  im¬ 
perative  that  .Ibrational  disturbances  be  avoided. 
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Immediate  objectives. 

After  test  runs  have  been  completed  and  the  regulation 
characteristics  of  the  electrolyte  circulation  system  have  been 
established,  the  investigation  will  be  conducted  along  the 
following  lines; 

A.  Employing  the  wide  els'" t- ode -channel  and  a  copper 
deposition  electrode  reaction. 

1.  Collection  and  correlation  of  data  for  lajninar 
pure  forced  convective  mass  transfer,  employing  the  cathode- 
over-anode  position. 

2.  Establishment  of  the  laminar- to- t.rbulent  transi¬ 
tion  Re  range  in  pure  forced  convective  flow. 

3.  Collection  and  provisional  correlation  of  data 
for  turbulent  forced  convective  mass  transfer: 

a.  below  Re  =  10,000 

b.  above  Re  =  10,000. 

4.  Evaluation,  over  the  total  Re  range  covered,  of 
the  bounda.  idltions  existing  at  individual  sections  and 
their  influence  o'  generalized  correlations. 

5.  Considerao.ion  of  modified  cathode  segmentation 
(1st  part) . 

6.  Collection  and  correlation  of  sectional  limiting 
current  data  for  combined  free  and  forced  convective  miss  trans¬ 
fer,  employing  the  anode -over-cathode  position  and  a  wide  range 
of  bulk  concentration. 

7.  Collection  of  deposition  profile  data  indicating 
onset  of  free  convective  flow  patterns  and  correlation  with 


. . .  . 
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section  data  obtained  for  the  initiation  of  free  convsct j  ve 
mass  transfer  regime. 

8.  Comparison  of  data  with  respect  to  the  Initiation 
of  full-scale  turbulence  for  pure  forced  convective  and  mixed 
convective  mass  transfer  situations. 

9.  Evaluation,  over  the  Re  range  covered,  of  the 
boundary  conditions  existing  at  individual  sections  in  mixed 
convective  mass  transfer. 

10.  Estimation  of  expected  current  density  distribution 
at  other  orientations  with  respect  to  the  vertical. 

11.  Consideration  of  modified  cathode  segmentation 
(2nd  part) . 

B.  Employing  the  narrow  electrode  channel  and  a  ferrl- 
ferrocyanide  reduction  reaction  on  nickel  cathode. 

1.  Collection  of  section  current  density  vs  overpoten- 
tlal  data  in  the  laminar  forced  convective  mass  transfer  regime. 

2 .  Evaluation  of  boundary  conditions  existing  at  in¬ 
dividual  sections,  and  the  reliability  of  overpotential  as 

an  indication. 

3.  Collection  aind  correlation  of  turbulent  convective 
mass  transfer  data  above  Re  =  10,000  with  existing  analogies, 
using  extreme  variations  of  electrolyte  Sc  number. 
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PROJECT  VII. 

Transport  Properties  of  Concentrated 
Electrolytic  Solutions 

In  order  to  make  accurate  mass  transfer  cal'^ulatlons  in 
electrolytic  systems,  it  is  necessary  to  have  correct  values 
of  the  transport  properties  of  the  electrolytic  solution.  Cases 
of  practical  Interest  usually  require  a  knowledge  of  the  diffu¬ 
sion  coefficients  of  Individual  species  as  well  as  conductance 
In  a  multicomponent  solution  over  a  range  of  temperatures  and 
concentrations.  Unfortunately,  few  such  data  are  available  In 
the  literature.  Most  existing  measurements  have  been  made  on 
binary  electrolytes  In  dilute  solution  at  25®C.  Because  of  these 
circumstances  a  program  has  been  undertaken  to  measure  and  to 
co'^relate  the  transport  properties  of  concentrated  electrolytic 
solutions  over  a  wide  range  of  conditions. 

Correlation  of  transport  properties 

The  first  step  In  systematizing  the  information  on  such 
solutions  was  to  construct  a  rational  mathematical  framework 
of  flux  equations  wherein  all  the  transport  properties  are 
unambiguously  defined  and  ascribed  with  a  specific  significance. 
Newman,  Bennlon,  and  Tobias^  have  pointed  out  the  shortcomings 
of  the  classical  dilute  solution  flux  equation  with  regard  to 
concentrated  and  multicomponent  solutions.  They  proceeded  to 
eliminate  the  inherent  difficulties  oy  using  the  generalized 
multicomponent  diffusion  equations: 
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CiVu^  =  RT 


Q>  1  C  4 

(ij  ■  ii* 
VlJ  ^ 


where  the  transport  properties  D^j  each  represent  the  inter¬ 
action  between  a  pair  of  species  in  the  solution.  It  is  seen 
that  this  relationship  can  apply  to  an  arbitrary  number  of 
components.  Also  the  formulation  avoids  a  frame  of  reference 
problem  usually  encountered  in  measuring  properties  of  con¬ 
centrated  solutions.  It  Is  felt  that  because  of  the  straight¬ 
forward  significance  of  the  they  should  lend  themselves 

more  readily  to  theoretical  interpretation  and  correlate  more 
easily  than  the  properties  defined  by  the  classical  flux  equa¬ 
tions.  This  is  the  problem  which  is  presently  under  considera¬ 
tion. 


Before  dealing  with  the  complicated  problem  of  multicom¬ 
ponent  systems,  it  is  necessary  to  investigate  first  the 
behavior  of  the  interaction  coefficients  in  binary  solu¬ 
tions  of  various  electrolytes  in  order  to  determine  their 
dependence  on  the  various  physical  parameters  of  the  system. 

All  available  data  on  densities,  viscosities,  activity  coeffi¬ 
cients,  conductances,  dlffuslvltles,  and  transference  numbers 
for  binary  systems  have  been  collected  from  the  literature. 

Care  has  been  taken  to  critically  evaluate  the  data  where  possi¬ 
ble.  It  has  been  found  that  a  fair  amount  of  data  are  avall- 
abj.e  for  most  of  these  properties,  at  least  at  25®C.  However, 
for  dlffuslvlty  and  transference  number  the  data  are  very  sparse 
for  concentrated  solutions.  It  should  also  be  noted  that 
diffusion  data  for  multicomponent  systems  are  practically  non¬ 
existent  . 
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When  all  relevant  data  are  available  for  a  particular  bi¬ 
nary  system,  it  is  possible  to  calculate  the  three  D. ^  coeffi- 
dents  as  was  done  by  Newman  d  d.  for  potassium  chloride.  The 
raw  data  for  their  calculations  are  shown  in  Figure  1.  It 
is  seen  ti.at  conductivity  increases  monatonically  with  concen¬ 
tration,  the  cation  transference  number  decreases  with  concen¬ 
tration,  and  the  measured  diffusion  coefficient  drops  rapidly 
to  a  minimum  and  then  rises  whereas  the  thermodynamic  diffusion 
coefficient  goes  in  the  opposite  direction.  This  general  be¬ 
havior  appears  to  be  characteristic  of  the  alkali  halides. 

The  calculated  functions  appear  in  Figure  5  of  the  previous 


reference , 

In  Project  IV.  some  tentative  assumptions  about  the  depen¬ 
dence  of  on  concentration  and  viscosity  led  to  a  conducti¬ 
vity  correlation  based  on  the  equation 

AB 


^/  vcc 


-  z_j_z_eF 


VC. 


V^i 


=  G  + 


+  v_B 


(2) 


T 


If  A,  B,  and  G  from  their  equations  (12)  are  constant,  a  plot  of 
vcc^/kp  versus  ■/" should  be  linear.  It  was  found  that 
for  KCl  and  KOH  the  suggested  relationship  was  indeed  linear 
over  the  entire  concentration  ranges  at  25®C.  This  relationship 
has  been  investigated  for  a  number  of  other  systems. 

In  Figure  2  are  plotted  the  conductivities  of  ammonium 
nitrate  solutions  of  various  compositions  and  at  various  tem¬ 
peratures.  It  is  seen  that  for  this  salt  the  conductivity 
ai:  any  temperature  goes  through  a  maximum  in  the  high 
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concentration  range.  This  behavior  is  common  for  a  number  of 
electrolytes  including  KOH.  It  is  such  a  maximum  which  often 
determines  the  choice  of  concentration  of  supporting  electro¬ 
lyte  In  a  working  cell.  The  phenomenon  appears  to  be  related 
to  the  large  viscosity  increase  observed  in  concentrated  solu¬ 
tions  of  many  electrolytes. 

In  Figures  3  and  4  these  data  for  NH^NO^  are  plotted  in 
the  manner  suggested  by  equation  (2).  A  linear  relationship 
appears  to  hold  well  over  most  of  the  concentration  range  but 
is  seen  in  Figure  4  to  fall  at  high  concentrations.  Also  the 
data  for  different  temperatures  fall  on  different  lines  which 
indicates  that  the  temperature  dependence  of  has  not  been 
accounted  for  completely. 

In  Figures  5,  6,  and  7  tiiree  more  examples  are  presented. 
The  conductivity  of  CuSO^  when  plotted  In  Figure  6  is  seen  to 
deviate  from  linearity  In  the  concentrated  region  as  do  the 
data  for  LINO^  in  Figure  7.  In  Figure  8,  however,  it  is  seen 
that  LlNOj  In  ethanol  follows  the  linear  relationship  quite 
well. 

It  has  been  found  from  the  analysis  of  the  conductivity 
data  for  a  large  number  of  systems  that  the  forms  assumed  in 
equations  (12)  of  Nevnnan  et  al« account  quite  well  for  the  be¬ 
havior  of  in  dilute  binary  electrolytes  but  that  they 
sometimes  fal^.  in  the  concentrated  region.  The  reason  for 
success  or  failure  in  a  particular  case  Is  not  obvious.  How¬ 
ever,  it  dees  appear  that  these  first  approximations  do  take 
into  account  the  proper  effects  and  do  successfully  reduce 
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the  experimental  conductivity  data  to  a  well  behaved  form. 

As  suggested  in  the  original  presentation  a  more  sophisticated 
approach  to  the  structure  and  interactions  in  the  solutions 
should  make  improvements  in  the  correlation  possible.  This 
problem  is  being  given  further  consideration. 

In  Figure  9  are  shown  the  cation  transference  number, 
measured  diffusion  coefficient,  and  thermodynamic  diffusion 
coefficient  of  copper  sulfate.  With  these  data  and  the  con¬ 
ductivity  shown  in  Figure  5  it  is  possible  to  calculate  the 
®1J  copper  sulfate.  The  results  are  shown  in  Figure  10. 

Similarly,  the  diffusion  coefficient  dat  \  for  ammonium 
nitrate  are  plotted  in  Figure  11.  Unforti... lately,  no  measure¬ 
ments  of  the  transference  number  of  this  salt  have  been  mr  '''=. 

By  assuming  that  =  0.487  it  was  possible  to  compute  the 
D^j’s  for  this  substance.  The  results  are  shown  in  Figure  12. 

By  comparing  Figures  10  and  12  with  Figure  5  of  Newman 
et  al.,  one  sees  that  the  0^^  of  these  three  systems  behave 
qualitatively  in  the  ssime  way.  However,  a  quantitative  corre¬ 
lation  is  not  yet  obvious  for  these  functions.  In  particular, 

the  behavior  of  D  for  the  three  salts  differs  considerably 

+  — 

in  the  high  concentration  range,  and  the  slopes  are  different 
even  in  the  dilute  range.  Again  it  appears  that  successful 
correlation  will  require  consideration  of  liquid  structure  and 
molecular  interactions  in  conjunction  with  Investigation  of 
the  behavior  of  many  other  binary  electrolytic  solutions. 

During  the  collection  from  the  literature  of  the  data 
which  are  necessary  for  these  calculations,  it  has  become 
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obvious  how  scarce  reliable  values  of  diffusion  coefficients 
for  concentrated  solutions  are.  Even  where  measurements  have 
been  made,  one  must  be  careful  in  ur.Ing  the  reported  values 
because  in  many  cases  they  are  no^  the  Intrinsic  properties 
which  are  defined  by  a  differential  flux  relationship.  Some¬ 
times  the  reported  diffusion  coefficient  is  an  integral  value, 
that  is,  one  which  represents  some  average  over  the  concentra¬ 
tion  difference  employed  in  the  measurement  rafter  than  a  true 
differential  value.  As  indicated  before,  it  is  sometimes  not 
clear  what  the  frame  of  reference  is  with  respect  to  vihich  the 
diffusion  coefficient  was  measured.  And  in  some  cases  the 
experiment  was  simply  not  capable  of  very  high  accuracy.  Dis¬ 
crepancies  between  different  Investigators  make  this  fact 
obvious.  Because  of  ^he  need  for  much  more  extensive  accurate 
data,  both  for  practical  application  and  as  a  basis  for  a 
successful  theory  or  correlation  for  concentrated  and  multi - 
component  solutions,  a  program  of  diffusion  coefficient 
measurements  is  being  undertaken.  In  conjunction  with  this, 

*  program  of  transference  number  measurement  is  also  being 
planned. 

Experimental  methods  for  determining  diffusivlty 

In  deciding  what  method  whould  be  chosen  for  the  diffusion 
studies,  we  have  considered  all  the  techniques  which  are 
presently  available.  Since  a  number  of  thorough  reviews  on  the 
experimental  methods  for  measurelng  diffusion  have  appeared 

2  3  4  5 

in  the  last  twenty  years  '  '  '  ,  no  detailed  discussion  seems 
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p  Justified  here.*  We  Si.jll  simply  outline  the  several  methods 
which  appear  to  be  reasonably  accurate  and  note  the  advantages 
and  disadvantages  which  are  particularly  relevant. 

Optical  methods:  Most  of  the  popular  optical  methods  are 
based  on  the  solution  of  Pick’s  law  for  free  diffusion  .  That 
is,  diffusion  taJces  place  from  an  initially  ■^harp  boundary  In 
an  effectively  infinite  column  of  solution.  Most  methods, 
such  as  the  Gouy  method  ,  use  Interferometry  to  observe  the 
gradient  of  index  of  refraction.  Such  a  measurement  can  be  made 
quite  accurately.  Unfortunately,  the  complete  analysis  of 
these  experiments  requires  the  assumption  of  a  constant  diffu¬ 
sion  coefficient.  Therefore,  it  is  necessary  to  work  with  a 
very  small  concentration  difference  which  tends  to  limit  pre¬ 
cision.  Furthermore,  the  analysis  of  the  measurements  depends 
rather  critically  on  the  initial  condition.  One  muat  be  very 
careful  to  maike  the  initial  concentration  boundary  as  sharp  as 
possible;  it  is  necessary  in  any  case  to  make  a  zero  time 
correction. 

These  methods  do  have  the  advantage  of  affording  measure¬ 
ments  in  a  relatively  short  time,  although  the  analysis  of  the 
fringe  patterns  appears  to  be  rather  tedluous.  In  spite  of  the 
various  difficulties  optical -methods  still  offer  the  possibility 
of  highly  accurate  measurements,  and  a  large  portion  of  re¬ 
liable  data,  particularly  in  concentrations  solutions,  has  been 
ODtalned  by  an  optical  technique. 


A  particular  complete  discussion  of  the  measurement  of  dlffu- 
sivlty  has  been  given  recently  by  Geddes®. 
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Diaphragm  cell;  A  method  which  Is  not  as  accurate  as  most 

of  the  optical  methods  but  which  is  easier  to  use  for  quick, 

9 

approximate  measurements  Is  the  diaphragm  cell  .  This  cell 
employs  psuedo-steady  state  diffusion  through  a  porous  dia¬ 
phragm  which  separates  two  large  volumes  of  solution.  The 
method  is  not  absolute,  and  each  diaphragm  must  be  calibrated 
with  some  substance  of  known  diffusivlty.  This  cell  measures 
an  integral  diffusion  coefficient,  but  an  analysis  has  been 
developed  which  makes  it  possible  to  calculate  the  differential 
coefficient  from  the  measured  values^^.  Because  chemical 
analysis  is  possible,  this  type  of  cell  offers  a  means  for 
studying  multicomponent  diffusion. 

It  seems  that  the  primary  alsadvantage  of  this  method, 
aside  from  its  limited  accuracy,  is  that  one  may  make  measure¬ 
ments  without  being  aware  of  a  sizable  systematic  error.  For 
example,  the  effective  diffusion  path  length  which  is  deter¬ 
mined  by  calibration  may  change  from  solution  to  solution  be¬ 
cause  of  a  change  in  the  hydrodynamics  at  the  faces  of  the 
d laphragm . 

Conduc time trie  method;  A  method  has  been  developed  by 
Harned^^  which  has  made  it  possible  to  measure  diffusivities 
of  electrolytes  at  extremely  low  concentrations  where  the  values 
can  be  compared  with  those  calculated  from  the  Onsager-Fucss 
limiting  law.  The  excellent  agreement  obtained  has  been 
gratifying  in  verifying  the  theory  as  well  as  indicating 
high  experimental  accuracy.  It  should  be  noted  that  Harned’s 
higher  concentration  measurements  agree  quite  well  with  the 
results  of  the  Gouy  optical  method. 
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*  f  The  essential  features  of  Harned’s  cell  are  that  it  is  of 

I  finite  length  such  that  diffusion  is  of  the  restricted  type  and 
,  I  that  the  decay  of  a  concentration  gradient  with  time  Is  monl- 

1  tored  by  a  conductivity  difference  between  two  points  In  the 

I  cell.  Unfortunately,  this  method  cannot  be  used  as  accurately 
I  in  its  present  form  for  studying  many  electrolytes  in  the  con- 
f  centrated  region  because  of  the  maximum  in  the  conductivity 
I  versus  concentration  curve. 

I  After  consideration  of  these  and  the  other  techniques  for 

f  measuring  diffuslvlty,  it  has  been  decided  that  two  methods 
4  j  which  offer  specific  advantages  are  a  rotating  disk  electrode 

1  and  a  modified  Harned  type  cell  which  Incorporates  an  optical 

*■  1  technique.  The  latter  combines  the  mathematical  advantages  of 

I  the  Harned  cell  with  the  precision  of  an  optical  measurement 

§  in  concentrated  solutions;  the  electrochemical  method  offers  a 

I  means  of  studying  diffusion  in  concentrated  multicomponent  sys- 
I  terns  by  measuring  the  mass  transfer  rate  of  a  salt  to  an  elec- 

I  trode  in  the  presence  of  a  supporting  electrolyte, 

i  Rotating  disk  electrode 

I  The  rotating  disk  is  one  of  very  few  cases  of  forced  con- 

J  vection  where  complete  solution  for  the  velocity  profile  is 
I  possible.  The  important  conclusion  of  the  hydrodynamic  analy- 

I  sis  is  that  the  normal  component  of  the  velocity  is  Independent 

•  i 

I  of  position  on  the  surface  of  the  disk.  Therefore,  if  the  disk 
I  is  made  a  cathode  such  that  mass  transfer  takes  place,  it  is 

I  J 

-I  pDsslble  to  make  limiting  current  measurements  and  with  the  use 
I  of  the  known  convection  term  to  determine  the  diffusion 
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coefficient  of  the  species  reacting  at  the  electrode. 

These  results  have  been  known  for  some  time  for  the  case 

12 

of  a  constant  property  solution  and  low  rates  of  mass  transfer 
An  analysis  has  been  made  recently  in  this  laboratory  to  take 
into  account  the  effect  of  varying  solution  properties  as  well 
as  finite  interfacial  velocity  due  to  the  mass  transfer.  One 
conclusion  of  this  analysis  is  that  non-constant  properties  do 
not  affect  the  uniform  accessibility  of  the  surface.  The  quan¬ 
titative  results  make  more  exact  Interpretation  of  the  experi¬ 
mental  mass  transfer  measurements  possible. 

A  rotating  disk  electrode  and  related  apparatl  have  been 
assembled  and  will  be  used  to  measure  the  diffusion  coefficient 
of  copper  sulfate  in  sulfuric  acid  solutions  at  various  tem¬ 
peratures  and  concentrations.  A  photograph  of  the  apparatus 
is  shown  in  Figure  13.  In  testing  the  apparatus  several  diffi¬ 
culties  have  been  encountered.  It  was  found  that  under  certain 
conditions  the  observed  polarization  curves  were  unstable  and 
poorly  reproducible.  These  results  indicate  that  turbulence 
occurred  even  though  the  Reynolds  number  was  low.  In  these 
experiments  the  disk  was  actually  the  end  of  a  cylinder  as 
shown  on  the  right  in  Figure  14.  It  was  decided  that  the  disk 
should  be  more  streamlined.  Therefore,  the  shape  has  been 
modified  as  shown  on  the  left  in  the  same  figure.  It  is  hoped 
that  this  will  eliminate  the  problem  of  turbulence. 

Another  problem  in  the  analysis  of  the  copper  rotating 
disk  electrode  is  the  electrode  reaction  kinetics.  Because  of 
the  low  exchange  current  density  the  limiting  current  which  is 
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i  measured  depends  on  both  the  concentration  overvoltage  and  the 

activation  overvoltage.  A  complete  equation  for  determining 
#  j  the  diffusion  coefficient  of  copper  sulfate  which  takes  both 

these  effects  into  account  has  been  derived. 

Optical  restricted  diffusion  cell 

I  Some  of  the  advantages  of  the  Harned  type  diffusion  cell 

have  already  been  mentioned.  Another  advantage  is  that  the 
measurement  is  independent  of  the  initial  concentration  distri¬ 
bution  since  actual  determination  involves  an  extrapolation  to 
infinite  time.  Therefore,  the  special  care  which  is  required 
f  in  setting  up  most  optical  methods  is  not  necessary  with  this 

cell . 

As  Indicated  previously,  it  is  difficult  in  many  experi¬ 
ments  to  determine  the  true  differential  diffusion  coefficient 
because  of  the  variation  of  solution  properties  with  concentra¬ 
tion.  This  complication  also  obscures  the  frame  of  reference 
in  many  cases.  In  the  Harned  cell,  however,  it  is  possible  to 
eliminate  these  difficulties. 

Harned’ s  solution  of  Pick’s  second  law  for  the  concentra- 
Gion  profile  in  a  column  of  liquid  of  height  was  made  for  a 
solution  with  constant  properties,  and  he  accordingly  carried 
out  his  measurements  using  small  concentrations  differences, 
liy  measuring  conductivity  to  monitor  concentration  at  z=a/6 
and  z=5a/6,  he  was  able  to  reduce  the  solution  of  the  equation  to 

io  =  o(|)  -  o(f-)  =  2A^  coa  (J)  e"  Dt(ir/a)^ 

+  2A5  003  (It)  e-  25Dt(F/a)^  +  .  .  .  .  {3I 


Section  D.,  Project  VII. 


Page  212. 


It  Is  seen  that  a  series  of  this  form  converges  very  rapidly 
with  time  so  that  at  large  times  a  plot  of  lx  (Ac)  versus  t 
will  be  a  straight  line,  the  slope  of  which  yields  the  diffu¬ 
sion  coefficient. 

Because  of  the  several  advantages  of  this  type  of  measure¬ 
ment,  we  have  investigated  the  possibility  of  modifying  Harned’s 
cell  for  application  to  the  study  of  concentrated  solutions. 

To  Qo  this  it  is  first  necessary  to  show  that  the  diffusion 
coefficient  which  is  obtained  is  a  differential  one. 

Bj  considering  the  complete  flux  equation  for  concentrated 
solutions  derived  from  equation  (1),  which  allows  for  flux  of 
solvent  as  well  as  solute  in  the  restricted  cell,  and  by 
treating  the  concentration  and  all  variable  properties  as  per¬ 
turbation  series  about  their  respective  values  at  infinite  time, 
it  has  been  possible  to  derive  the  equation  which  corresponds 
to  equation  (3)  but  which  is  completely  general  and  complete 
for  this  geometry: 


Ac  =  An 


-d”  (TT/a)  t 


+  A, 


-3D”  (F/a)^t 


(4) 


In  this  equation  xs  a  well-defined  differential  diffusion 
coefficient  corresponding  to  the  uniform  cell  concentration 
at  infinite  time.  The  coefficient  A^  is  a  complicated  combina¬ 
tion  of  concentration  derivatives  of  various  solution  proper¬ 
ties;  this  terra  vanishes  in  the  constant  property  case.  Even 
though  the  series  in  equation  (4)  does  not  converge  as  rapidly 
as  that  in  equation  (4),  it  obviously  will  converge  at  suffi¬ 
ciently  long  times.  These  times  will  be  on  the  order  of  hours 
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}  to  days,  depending  upon  the  magnitude  of  cell  height.  Thus 

It  Is  seen  that  even  In  the  case  of  non-constant  properties 
#  and  concentrated  solutions  a  plot  of  in  (Ac)  versus  t  should 

j  still  yield  a  linear  region  from  which  a  differential  diffusion 
I  coefficient  can  be  obtained. 

I  It  was  pointed  out  that  In  many  cases  conductivity  cannot 

be  used  to  monitor  a  difference  In  concentration  in  a  concentra¬ 
tion  solution.  It  Is  felt  that  it  would  be  well  to  make  use 
I  here  of  the  high  accuracy  of  an  Interfermotric  technique  to 
observe  the  behavior  of  the  solution.  To  be  specific,  it  Is 
i  planned  that  multiple  beam  Interference  fringes  will  be  generated 

In  the  cell.  Then  as  diffusion  takes  place,  the  contour  of 
these  fringes  will  give  a  direct,  precise  quantitative  indica¬ 
tion  of  the  concentration  profile  in  the  cell  as  it  decays  to 
uniformity.  In  this  way  Ac  of  equation  (4)  can  be  determined 
as  a  function  of  time. 

It  Is  hoped  that  a  cell  such  as  the  one  just  described 
can  be  perfected  to  give  highly  accurate  diffusion  coefficients 
for  concentrated  solutions.  If  this  can  be  done,  such  a  cell 
then  offers  the  possibility  for  some  systems  of  combining 
Harned's  electrodes  with  the  optical  measurement  suggested 
here  to  study  multicomponent  diffusion. 

In  summary,  plans  have  been  made  to  obtain  more  data  on 
the  transport  properties  of  concentrated  electrolytic  solu¬ 
tions  In  order  to  Investigate  the  general  behavior  of  the 
coefficients  of  the  multicomponent  diffusion  equation. 

In  this  way  It  Is  hoped  that  we  shall  discover  how  the  various 

I 
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physical  parameters  of  the  system  affect  their  behavior  and 
thereby  develop  some  means  of  correlating  the  values  of  the 


s  in  binary  systems. 


Prom  such  information  attempts  will 


then  be  made  to  correlate  and  predict  the  behavior  of  the 


D. ,’s  in  multicomponent  systems. 
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CONCENTRATION  (  MOLES/LITER  ) 


CONCENTRATION  (  MOLES/LITER  ) 

Figure  1.  The  conductivity,  transference  numbers  and  diffusion 
coefficients  of  potassium  chloiide  in  water  at  25®C. 
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Figure  10.  The  multicomponent  diffusion  coefficients  of  copper  sulfate  In 
water  at  25® C  over  the  entire  concentration  range. 
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Figure  13.  Photograph  of  the  rotating  disk  apparatus. 


Figure  l4.  Photograph  of  the  rotating  disk  electrode 
_ (actual  size).  On  the  right  is  the  old  disk. 
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g  E.  CLOSING  REMARKS  AND  RECOMMENDATIONS 

j  ’hne  seven  research  projects  supported  in  full  or  in  part  by 

I  this  contract  were  aimed  at  gaining  quantitative  understanding 
I  of  the  role  of  mass  transport  processes  in  relation  to  charge 
I  transfer  reactions  at  electrodes.  Although  the  models  chosen 
I  for  study  are  illustrative  of  certain  of  the  conditions  pre- 
I  vailing  in  fuel  cells  and  primary  and  secondaiy  batteries, 

I  clearly  the  conclusions  reached  can  only  be  applied  to  a  g^ven 
I  practical  situation  by  Intelligent  examination  of  all  charac- 
I  teristic  system  parameters. 

I  lu  is  evident  to  those  who  performed  this  research,  as  it 

I  will  be  to  the  reader,  that  further  sustained  research  efforts 
I  related  to  mass  transport  are  needed  to  sUpport  a  rational 

3 

I  development  of  future  hardward  devices.  The  following  areas 

■M 

§  deserve  particular  attention;: 

I  1.  The  effect  of  flow  past  electrode  surfaces  and  through 

I  porous  electrodes  on  local  and  overall  overpotential,  and  on 

I  current  distribution,  has  been  characterized  only  for  a  restrlc- 

I  ted  range  of  flow  and  geometric  variables.  In  particular  the 
I  laminar  to  turbulent  transition  regime,  pulsating  flows  in 
I  parallel  and  through  flow,  need  to  be  considered.  Although 

I  transport  in  free  convection  has  been  elucidated  successfully 

I  for  extended  surfaces,  further  efforts  should  be  directed 
I  toward  the  study  of  free  convection  effects  where  the  flow  is 

I  confined  by  cxosely  spaced  solid  boundaries,  such  as  between 

I  battery  plates,  or  in  macropores. 
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2 .  Work  under  this  contract,  tcgether  with  most  recenu 
results  here  and  abroad,  on  the  elucidation  of  the  mechanism 
of  porous  gas  electrodes  clarified  the  most  elementary  aspects 
of  the  role  of  Ionic  and  gas  transport.  However,  so  far  the 
models  studied,  both  In  respect  to  geometry  and  materials 
employed,  and  In  regard  to  the  simplifying  assumptions  made, 
can  only  serve  as  qualitative  guides  for  further  more  detailed 
and  more  exact  studies.  The  effect  of  Interfaclal  tension  as 
It  Is  relate'^  to  electrocaplllarlty  phenomena,  has  not  been  as 
yet  considered,  neither  has  It  been  even  qualitatively  shown 
v/hether  convection  due  to  surface  forces  plays  a  significant 
role  In  porous  gas  electrodes.  Although  it  was  clearly  demon¬ 
strated  that  when  water  Is  a  reactant,  its  supply  must  occur 
at  least  to  a  large  degree  through  the  gas  phase  (and  not  as 
It  was  believed  before,  from  the  bulk  electrolyte  by  diffusion 
Into  the  film-reaction  zone),  this  aspect,  because  of  Its  Im¬ 
mediate  significance  In  fuel  cell  design,  should  be  the  3ubj‘:;ct 
of  more  detailed  scrutiny  Involving  both  simple  experimental 
models  and  actual  hardware  type  porous  elfctrodes. 

3.  The  degree  of  usefulness  of  understanding  of  mass 
transport  processes  depends  on  the  availability  of  data  on 
equilibrium  and  transport  properties.  As  long  as  solubilities, 
viscosities,  conductances,  transference  numbers  and  diffusion 
coefficients  of  cell  reactants  and  products  In  the  concentra¬ 
tion  range's  that  are  typical  of  actual  hardware  applications 
a.'e  missing,  the  design  and  optimization  of  every  piece  of 
equipment  will  remain  a  costly  "Edlsonian"  task.  The 
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determination  of  the  solubility  and  diffusion  coefficient  of 
oxygen  .n  KOH,  carried  out  under  this  contract  was  long  over¬ 
due,  and  far  from  being  enough  to  allow  rational  development 
of  the  many  types  of  systems  of  practical  Interest  involving 
a  variety  of  electrolytes,  reactants  and  products.  Direct 
measurements  of  these  properties  over  representative  concen¬ 
tration  ranges  are  extremely  time-consuming,  and  therefore 
the  development  generalized  methods  for  their  prediction  should 
receive  high  priority. 

The  research  program  initiated  under  this  contract  during 
its  last  year  is  being  continued  under  the  sponsorship  of 
the  Inorganic  Materials  Research  Division  of  the  Lawrence 
Radiation  Laboi  y.  (Projects  VI  and  VII.)  A  new  project 
was  initiated  recently  as  a  logical  outgrowth  of  Project  I, 
aimed  at  the  elucidation  of  the  mechanism  of  electrolytic 
gas  evolution.  As  a  recognition  of  its  support  during  the 
initial  phases  of  these  projects,  the  U.S.  Array  Electronic 
Research  Laboratory  will  be  kept  Informed  about  the  progress 
of  this  work. 
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Basis  for  evaluation  of  work-hours  -  22  days  (8-hr)  per  month. 

The  University  provided  compensation  for  research  done  during  the  academic  year  (9/15-6/15). 

This  figure  is  obtained  by  subtracting  1-1/4  day  vacation  time /full  work  month 
(l.e.,  10  hrs.  out  of  each  176  hrs.). 
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